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SUMMARY 


Presented in this report are the results and methodology used to derive linear 
models from a nonlinear simulation. It is shown that bias errors in partial derivative 
calculations using the method of finite differences can cause large enough errors in ma- 
trix elements such that eigenvalues of the system can be changed. These errors can be 
reduced and in some cases eliminated through positive and negative perturbations about 
the operating point. Linear models developed in this manner can show the same eigen- 
values for different finite-difference perturbations within the linear range about the op- 
erating point. 

It is also shown in this report that positive and negative perturbations in the system 
inputs can cause the linear model to represent better the nonlinear simulation in both 
directions about an operating point than when only a single perturbation is used. 

As an illustration of these positive and negative perturbation techniques, linear 
models of the F100 engine were derived from the nonlinear engine simulation. Compar- 
isons of step-response transients run on a nonlinear simulation with those generated 
using the linear models were very good. It is also shown that, because of the small 
perturbations required for transient comparisons, great care must be used to insure 
that there is as little startup transient in the nonlinear simulation as possible. A large 
startup transient can completely overshadow the effects of a small input. 

Finally, techniques to obtain reduced-order linear models from full-state linear 
models were investigated. The modal- and normal-reduction techniques were applied to 
the full-state linear models of the F100 engine. Results from transients run on the full- 
state and reduced-order models compared quite well. 


INTRODUCTION 

In recent years much work has been done on the application of multivariable (op- 
timal) control theory to the design of controls for many applications; specifically, for 
aircraft propulsion systems. Emphasis on this type of work has resulted from a desire 
to obtain improved engine performance over a wider operating envelope. This desire 
has necessitated the development of increasingly complex engines having more control 
variables and requiring more engine measurements. A great deal of multivariable con- 
trol work has been done on turbofan engines. For example, the Pratt & Whitney F100 
turbofan engine has variable fan and compressor geometry, engine bleeds, gas-generator 



fuel flow, and exhaust-nozzle area. Experience gained from this engine can be used to 
develop multivariable controls for even more complex variable-cycle engines (VCE). In 
addition to the controls already mentioned, VCE may have variable turbine geometry, 
duct-burner fuel flow, and a separate duct nozzle. The proper control of these engines 
will require the measurement of even more engine parameters. 

Work currently done using multivariable (optimal) control theory is given in refer- 
ences 1 to 4. In reference 1 linear quadratic regulator (LQR) theory was used to derive 
optimal control of the F100 engine. The controls developed in reference 1 were evalu- 
ated on a real-time, hybrid simulation of the F100 in reference 2. In reference 3 con- 
trols were designed for the F100 engine using output feedback regulator (OFR) theory. 
Minimum-time acceleration control histories are derived for the F100 engine using non- 
linear programming in reference 4. 

All the optimal control design techniques described in these references require 
linear engine models, since no body of optimal control theory has been developed for 
nonlinear systems. Thus, a good linear representation of the engine system must be 
obtained. Since most simulations developed to describe these physical systems are non- 
linear, linear models that are valid about an operating point must be derived. Also, be- 
cause of the nonlinear characteristics of the engine system, more than one linear model 
is necessary, in general, to describe the nonlinear system. For example, for the F100 
engine six linear models are needed to represent the engine over its wide operating 
range . 

A method used to linearize the nonlinear system about an operating point is the 
M offset derivative 11 method. This technique makes use of the method of finite differen- 
ces to approximate partial derivatives about an operating point. For the F100 studies 
previously described, the linear models were derived in reference 5 from the nonlinear 
F100 digital simulation described in reference 6. 

A problem that arises in using the method of finite differences is that small numer- 
ical errors (bias errors) can occur in the partial-derivative approximations. This re- 
port describes the derivation of linear models from a nonlinear simulation using posi- 
tive and negative perturbation techniques so that bias errors are minimized. These per- 
turbation techniques were applied to the F100 nonlinear simulation, and new linear mod- 
els were derived. These models are different from those derived in reference 5 due to 
the different perturbation technique and due to many changes and updates having been 
added to the nonlinear simulation. 

Once the linear (full state) models of the system have been determined, it may be 
possible to reduce the order of the linear models while still retaining the important dy- 
namic characteristics of the system. This is desirable for two reasons: First, a lower 
order system is easier to handle mathematically; second, controls developed using op- 
timal control techniques (such as LQR) require that all specified states of the system be 
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used to derive the optimal control. However, some of the states are not measurable and 
thus cannot be used for control purposes. Also, if all the states of a large-order sys- 
tem were measurable, the control derived would be too complex to implement. Two re- 
duction techniques were investigated - modal reduction (ref . 1) and normal reduction 
(refs. 7 and 8). Both techniques are applied to the F100 full- state linear models derived 
in this report. Reduction of the original linear models is described in reference 1. 

Once the full- state and reduced-order linear models are generated, they should be 
evaluated as to how well they approximate the dynamics of the nonlinear system about 
the operating points for which they were generated. As an example, this evaluation is 
made for the F100 linear models by comparing transients for the full-state and non- 
linear engine models and the full-state and reduced engine models. In doing this evalu- 
ation the effects of a nonlinear simulation startup transient on small-perturbation input 
response is discussed. 


ENGINE DESCRIPTION 

As an example of the linear-model generation technique, a nonlinear engine simula- 
tion is used. This simulation is of the Pratt & Whitney F100-PW-100 (3) engine schema- 
tically shown in figure 1. The engine is an axial, mixed- flow, augmented, twin-spool, 
low-bypass- ratio turbofan. A single inlet is used for both the fan and engine core air- 
flows. Airflow leaving the fan is separated into two streams: one passes through the 
engine core and the other through the annular fan duct. The three-stage fan is connected 
by a through-shaft to the two-stage, low-pressure turbine. A 10-stage compressor is 
connected by a hollow shaft to the two-stage, high-pressure turbine. The fan has var- 
iable, trailing-edge, inlet guide vanes. The compressor has a variable-inlet guide 
vane followed by two variable stator vanes . Engine airflow bleed is extracted at the 
compressor exit and discharged through the fan duct during starting. Compressor- 
discharge bleed air is also used to cool the high- and low-pressure turbine blades and 
to power the augmentor turbopump. The main combustor consists of an annular dif- 
fuser and a chamber with 16 full nozzles. The engine core and fan duct streams com- 
bine in an augmentor and are discharged through a variable, convergent-divergent noz- 
zle. The augmentor consists of a diffuser section and five concentric manifolds. 


ENGINE SIMULATION 
Nonlinear Model 


The nonlinear digital simulation used in this report is Pratt & Whitney's digital 
simulation CCD 1103-1.0 (ref. 6). The simulation uses overall performance maps to 
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describe steady-state component performance. Performance maps are provided for the 
fan, fan duct, compressor, burner, high-pressure turbine, low-pressure turbine, af- 
terburner, and exhaust nozzle. Dynamic relationships, based on the laws of conserva- 
tion of mass, energy, and momentum are incorporated in the nonlinear simulation. 
These include rotor inertial effects, volume gas dynamics, and heat-transfer effects. 
The dynamic elements are modeled with recursion formulas. The digital simulation in- 
cludes a simulation of the engine control system. The control simulation supplies fuel 
flow, exhaust-nozzle area, compressor variable-vane position, compressor bleed, and 
fan-inlet-guide-vane position inputs to the engine simulation. These control inputs are 
computed as functions of various engine measurements and the pilot-set power level 
angle (PLA) command. 


Nonlinear State-Space Model 

As mentioned earlier, the nonlinear dynamic equations in the simulation include 
rotor inertial effects, volume gas dynamics, and heat-transfer effects. The physical 
quantities represented by these equations are rotor speeds, volume pressures, and 
temperatures. The nonlinear equations are described in nonlinear state- space form; 
that is, they are of the form 


X = f (x,u,t) 


( 1 ) 


and 


Y = g (x, U, t) (2) 

(All symbols are defined in appendix A.) The state variables X cannot change instan- 
taneously with time. The output variables Y, which are important engine parameters 
but not necessarily states, can change instantaneously with time. In the nonlinear en- 
gine simulation, there are 16 of these nonlinear state equations. The 16 states corre- 
sponding to these state equations are 

X^ fan speed, 

X 2 compressor speed, Npj 

X 3 compressor-discharge pressure, P T3 

X^ interturbine pressure, P,p^ ^ 

X^. augmentor pressure, P^j^ 

Xg fan-inside-diameter discharge temperature, T^ 


4 


duct temperature, T p2 5 ^ 

Xg compressor-discharge temperature, T T3 

Xg burner-exit fast- response temperature, T^ 4HI 

X 1Q burner-exit slow-response temperature, T T4LO 

X ^ 4 burner-exit temperature, T p4 

fan-turbine-inlet fast- response temperature, T ^ 4 ^jjj- 
X^g fan-turbine-inlet slow-response temperature, T p4 5Pi q 
X ^ 4 fan-turbine-exit temperature, T,p,- 

X-^g duct-exit temperature, T^g^, 

X 16 exhaust-nozzle-inlet temperature, T^^ 

The five control variables are 

main-burner fuel flow, WFMB 
U 2 exhaust-nozzle-jet area, AN 

U 3 fan-inlet-guide-vane position, CIW 

U 4 compressor variable-vane position, RCVV 

U 5 compressor bleed-flow fraction, BLC 

The afterburner fuel flow was not included in these studies. The output variables are 
engine net thrust, FNMX 
Yg total engine airflow, WFAN 

Yg burner-exit temperature, T p4 

Y 4 fan stall margin, SMAF 

Yg compressor stall margin, SMHC 

Y g empirical fan-exit AP/P parameter, DP 2 & q P 

Y 7 theoretical fan-exit AP/P parameter, DP 2 5gc 

Linear State- Space Models 

To derive optimal controls for the F100 engine, a linear engine model is required 
since no body of optimal control theory has been developed for nonlinear systems. 
Thus, a linear representation of the engine must be obtained from the nonlinear state 
and output equations (eqs. ( 1 ) and ( 2 )).- Also, no single linear model can accurately 
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represent the engine because of its wide operating range and nonlinear characteristics. 
Thus, linear models must be derived at various flight conditions and power settings 
throughout the engine operating envelope. The criteria used to select the various oper- 
ating points is that a combination of the linear models generated at these operating 
points must adequately describe the dynamic operation of the engine over the whole en- 
velope. Figure 2 shows a typical engine operating envelope and the flight conditions 
chosen within the envelope to describe engine operation. Table I gives the reasons for 
choosing these flight conditions. (Table I and fig. 2 are also found in ref. 9.) 

Nonlinear equations (1) and (2) can be linearized to yield the following linear state- 
space equations 

X = A X + B U (3) 


Y = C X + D U (4) 

The mathematics used to derive equations (3) and (4) from (1) and (2) are given in ap- 
pendix B. The matrices A, B, C, and D are constant at an operating point, thus the 
models are linear and time invariant. The linear model is shown in block diagram 
form in figure 3. 

The "offset derivative" method with forced steady-state balance was used to derive 
the linear models from the nonlinear engine simulation at the selected engine operating 
points . The method of finite differences was used to obtain approximations of the sys- 
tem partial derivatives at an operating point (ref. 5). 

A and C matrix generation . - The use of finite differences to approximate partial 
derivatives requires that each state and input variable be perturbed individually, while 
all others are held constant. The A and C matrices are generated by perturbing the 
states. From appendix B 


ax 

A.. = — l - 
3 9X j 


( 5 ) 


and 


8Y i 
C.. = — - 

1J ax. 


( 6 ) 


Next, the linearity of these matrices must be checked. This is done by first generating 
an A matrix with one perturbation size and then a second A matrix A' with a differ- 
ent perturbation size. If both A and A' have the same elements, the perturbation 
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sizes are in the linear range about the operating point. In practice, however, it is dif- 
ficult to compare all the elements of the matrices, especially if the system is large. 
Also, because of iteration tolerances in the simulation, the elements could vary from 
matrix to matrix but still retain valid information of the linearity of the system. A 
more practical method of checking for linearity is to calculate the matrix eigenvalues: 



If X approximates X' within the frequency range of interest, both perturbation sizes 
are within the linear range about the operating point, and both A and A' are valid. 

The choice of which perturbation size to use is arbitrary. The smaller perturbation 
size may cause numerical problems in calculating the partial derivatives, especially if 
it is not much larger than the iteration tolerance size specified in the simulation. The 
larger perturbation size insures easier calculations ofthe partial derivatives, but it may 
approach the linearity limit about the operating point. Further restrictions on the per- 
turbation size for the states are given in appendix C. 

This procedure was performed on the nonlinear F100 simulation. The perturbation 
sizes picked were +0.5 and +1 percent of the steady-state values of states, based on 
work done in reference 5. The operating point selected was at sea-level- static idle, 
which was one of the six operating points selected for control studies in reference 1 . 
This operating point was selected for evaluation because it results in an unchoked noz- 
zle, which can cause computation difficulties. This aggravates the problem of calcu- 
lating finite-difference partial -derivative approximations at this condition. Thus, if 
accurate linear models can be obtained at this operating point, no additional problems 
should be encountered at the other operating points where the nozzle is choked. 

Results for the +0.5 and +1.0 percent perturbations at the sea-level-static-idle 
operating point are given in table 11(a). Note that the eigenvalues are significantly dif- 
ferent for the two perturbation sizes. This is shown graphically in figure 4(a) for the 
five lowest eigenvalues. From this, one real eigenvalue and one complex pair were the 
same for both perturbation sizes. However, the +1 percent perturbation showed two 
real eigenvalues at -3. 97 and -8. 15 radians per second, and the +0. 5 percent perturba- 
tion showed a complex pair. Other perturbation sizes, both larger and smaller, were 
investigated; however, no repeatability was found for the eigenvalues. 

One cause of this problem is that there are small errors (bias) in the finite- 
difference calculations. The bias errors are discussed in detail in appendix C. In the 
F100 nonlinear simulation implicit techniques, which require iteration for solution, are 
used to calculate the partial derivatives . The bias errors occur when implicit techniques 
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are used, and the iterations required to derive the partial derivatives converge within 
the specified tolerance but not at zero. 

A method for reducing the effect of this bias error was investigated. First, the 
states were perturbed in the positive direction from the operating point to generate A 
and C matrices (a + and C + ) . Then, the states were perturbed in the negative direc- 
tion from the operating point to generate new A and C matrices (a_ and C_ ). Fi- 
nally, the A + and A_ matrices and the C + and C_ matrices were averaged to ob- 
tain the A and C matrices. The rationale for doing this is explained in appendix C. 

Eigenvalues from A matrices determined in this manner are shown in table n(b) 
for the ±0.5 and ±1.0 percent perturbations about the sea-level-static, idle operating 
point. Note that there is still a shift in the eigenvalues. However, there is much more 
consistency in the low-frequency range because there is no shifting from complex to 
real eigenvalues or vice versa. The improved accuracy of the averaged A matrix is 
shown graphically for the lower frequency eigenvalues in figure 4(b). From this figure 
the biggest shift is a movement of a real eigenvalue from -12.2 to -14. 1 radians per 
second for the ±0.1 and ±0.5 percent perturbations, respectively. 

From table n (b), the larger errors in the eigenvalue correspondence are at the high 
end of the frequency spectrum where an eigenvalue shifted from -422 to -268 radians 
per second. While this is a large shift, it may be disregarded for most applications 
(especially controls) . Uncertainty about the exact nature of a high-frequency eigenvalue 
is more acceptable than uncertainty in the low-frequency range. 

It is interesting to note here that, as explained in appendix C, the iterative nature 
of the solution of the partial derivatives for an implicit formulation does not guarantee 
an exact eigenvalue match, since there is no guarantee that the bias errors introduced 
by the positive and negative perturbations will cancel when the partial derivatives are 
averaged. If the method of finite differences are used and the partial derivatives are 
explicitly calculated, the eigenvalue correspondence would be exact. Results from the 
±0.5 percent perturbation were used in subsequent analyses. 

B and D matrix generation . - Digital logic built into the nonlinear simulation de- 
rives the B and D matrices from the supplied A and C matrices by forcing a 
steady- state match with the nonlinear simulation results. The mathematics for doing 
this are presented in appendix B. For a spepified control perturbation and given A and 
C matrices, the B and D matrices may be calculated by setting the left side of equa- 
tion (3) to zero. Equations (3) and (4) thus become 

X = -A" 1 BU|^ =0 (8) 

Y = D - C A _1 B u|^ =() (9) 
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As long as the perturbations in the control vector are in the linear range about the op- 
erating point, the B and D matrices are valid. A method of checking the linearity of 
these matrices comes from their definitions. From appendix B 



9X i 


9Uj 


( 10 ) 


9Y i 
D.. = — - 

13 9U. 


J 


( 11 ) 


In the linear range about the operating point the A, B, C, and D matrices are constant. 
Thus, for two perturbation sizes in the linear range, the ratio of the steady- state or 
output vectors should be equal to the ratio of the control- input vectors (eqs. (8) and (9)). 
For the linear models generated in this report, the perturbation sizes for the controls 
were supplied in the nonlinear simulation. They were +3 percent for the fuel flow and 
nozzle area, +5° for the fan- inlet-guide- vane position (CIW), +1° for RCW, and 
0.2 percent for the customer-compressor-bleed-flow fraction (BLC). The reasons for 
the selection of these perturbation sizes are given in reference 5. However, in using 
these perturbation sizes it was found that slightly different steady-state values were ob- 
tained with the nonlinear simulation when perturbing in the positive and negative direc- 
tions; that is, for equal positive and negative perturbations, the resulting state and 
output derivations were slightly different. For example, at the sea- level- static-idle 
operating point, a +3 percent change in fuel flow in the nonlinear simulation caused a 
change of 67.2 rpm in compressor speed; a -3 percent change gave a -73.4-rpm change 
in compressor speed. Thus, 


R. 


67.2 

73.4 


0.916 


( 12 ) 


R 


u 


+3 percent (WFMB) 
-3 percent (WFMB) 


(13) 


Since these differences are considered small compared with the actual compressor 
speed (9198. lrpm), the 3-percent perturbation is within the linear range about the oper- 
ating point. To obtain B and D matrices that are representative at the operating point, 
control perturbations were made in the positive direction to generate B + and D + and 
in the negative direction to generate B_ and D_. As with the A and C matrices, the 
B, and B and the D. and D matrices were averaged. The resultant linear model 
exhibited a 70. 3-rpm change in compressor speed for a 3- percent change in fuel flow at 
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the sea-level-static-idle condition. (Note that (67.2 + 73.4)/2.0 = 70.3 rpm.) 

It should be noted here that the derivation of the B and D matrices is application 
dependent. The averaging technique gives a better overall description of the operating 
point when the inputs are to be perturbed in both positive and negative directions. How- 
ever, for some applications the best model may be one that retains exact information 
about the system in only one direction. An example of this would be using the linear 
models to determine minimum time acceleration controls (ref. 4). However, in all 
cases, the A matrix should be as accurate as possible so that the actual system dy- 
namics are adequately described. 

To summarize, the final method used for generating linear models using the posi- 
tive and negative perturbation technique was as follows: 

(1) Matrices A + and C + were generated by a positive perturbation of the states. 
(For the F100, +0.5 percent of the states was used.) 

(2) Matrices A_ and C_ were generated by a negative perturbation of the states. 
(For the F100, -0.5 percent of the states was used.) 

(3) The matrices were averaged 


- A + + A _ 

A = (14) 

2 


C - 


C + + C _ 


(15) 


(4) Using A and C, matrices B + and D + were calculated by perturbing the con- 
trols in the positive direction. (For the F100, +3 percent WFMB and AN, +5° for 
CIVV, +1° for RCVY, and +0.2 percent for BLC . ) 

(5) Using A and C, matrices B_ and D_ were calculated by perturbing the con- 
trols in the negative direction. (For the F100, -3 percent WFMB and AN, -5° for 
CIW, -1° for RCVV, and -0.2 percent for BLC.) 

(6) The matrices were averaged 


- B , + B _ 

B = (16) 

2 


— D + + D _ 

D = (17) 

With this technique linear models can be derived from a nonlinear simulation. The 
models generated for the F100 engine are presented in appendix D for the six operating 


10 



points defined in reference 1. For each operating point, the A, B, C, and D ma- 
trices, the corresponding eigenvalues and eigenvectors, and the matrix of mode shapes 
are given. 


Linear- Model Validation 

Once the linear models are generated, they should be evaluated as to how well they 
approximate the dynamics of the nonlinear system about the operating points from which 
they were generated. This is accomplished by comparing transients run on both the 
linear and nonlinear simulations. For a given input the linear model should retain the 
dynamics of the nonlinear model for all the states and outputs, and the steady-state er- 
ror between the model responses should be small. 

This comparison was done for the F100 full-state models. A +3-percent fuel-flow 
and nozzle-area step transients were run on both the full-state linear and nonlinear 
simulations. Figures 5 and 6 compare responses of selected states and outputs with 
the fuel-flow and nozzle-area steps, respectively. These transients were run at the 
sea-level-static-idle operating point. Responses of the other states and outputs are 
compared in appendix E. In figures 5 and 6 as well as those presented in appendix E, 
the zero time point is the time at which the step input is applied. The nonlinear simu- 
lation was run for 0.1 second before the step was applied to eliminate any startup tran- 
sient. The importance of this will be discussed later. 

Figure 5(a) shows the response in fan speed to a + 3-percent step in fuel flow. The 
nonlinear simulation results are denoted by dashed lines, while the solid lines denote 
the linear simulation results. The agreement between these two fan-speed responses is 
good. The responses exhibit the same dynamic characteristics and have small steady- 
state differences. The observed steady-state error (7 rpm at 2 sec) in fan speed was at- 
tributed to the averaging done to derive the linear model. The steady-state error would 
have been negligible (due only to convergence tolerances in the nonlinear simulation) if 
the B and D matrices had been derived for only a +3-percent perturbation in the in- 
puts. However, the steady-state error would have been much larger than 7 rpm if tran- 
sients had been compared for a -3- percent perturbation in fuel flow. The compressor- 
speed response to a +3-percent fuel-flow step is shown in figure 5(b). Here also, the 
agreement is good, and the steady-state error is small -about 4 rpm at 2 seconds. Fig- 
ure 5(c) shows the burner-exit total-temperature response to the fuel- flow step. Here 
the agreement is excellent with the steady-state error of only 0.2 K at 2 seconds. The 
response of engine net thrust to a fuel-flow step is shown in figure 5(d). Here also, the 
agreement is very good. Note that at 0.05 second both responses show the same high- 
frequency dynamics with the nonlinear model undershoot somewhat larger than the 
linear-model undershoot. The steady-state error is 0.01 kilonewton at 2 seconds. 
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Figure 6(a) shows the fan-speed response to a +3-percent step in nozzle area for 
the nonlinear and linear simulations. In general, the responses are similar; however, 
the percentage of steady-state error (6 rpm) is somewhat larger than the fuel-flow step. 
Again, this is due to the averaging technique used to generate the B matrix. Note that 
the nonlinear response does not start at zero rpm at zero time. The small error is due 
to a startup transient in the nonlinear simulation. 

A startup transient is the result of a small error in an engine balance at zero time . 
In steady-state calculations an engine balance is forced by the simulation logic. How- 
ever, in transient calculations a different calculation path is generally followed through 
the simulation. Depending on how the tables are read and on how the quantities are cal- 
culated, small errors relative to the steady-state solution may result. Since an engine 
balance is not forced, the error tends to drive the system to a different operating point 
to achieve an engine balance (even if no input is applied). One method used to null out 
the startup transient is to run the simulation for a short time and let the engine balance 
before applying the desired input. This was done on the nonlinear simulation where it 
was run for 0.1 second before the step in fuel flow or nozzle area was applied. Thus, 
in all the figures that show a transient rim on the nonlinear simulation, the zero-time 
point is actually at 0. 1 second. For the linear model transients the zero-time point is 
actually at time zero. In figure 6(a) the offset in fan speed at zero time is thus due to 
the speed drifting slightly lower than steady state. Although the steady- state value is 
subtracted in an attempt to use zero as a reference point, the actual reference point is 
not zero. From figure 6(a) the error was 0.5 rpm, which is considered negligible. It 
should be noted that the same error also occurs in figure 5(a) but the startup transient 
effect is not apparent because of the large change in fan speed to a fuel-flow input. In 
the following comparisons an error at zero time due to a startup transient will be called 
a startup transient offset. 

Figure 6(b) shows the compressor-speed response to a step change in nozzle area. 
The responses agreed well with a steady-state error of 2 rpm at 2 seconds. The re- 
sponse of burner-exit total temperature to a nozzle-area step is shown in figure 6(c). 

At zero time the startup transient offset is 0.1 K. In general, both responses show the 
same dynamics with the linea'r response overshooting slightly higher and earlier than 
the nonlinear response. The error at 2 seconds is 0.2 K. Again, this error is due to 
the averaging technique used to derive the B matrix. Comparison of thrust responses 
from the nonlinear and linear simulations for a +3-percent step in fuel flow is shown in 
figure 6(d). Thrust is an output variable and can change instantaneously. At zero time 
the nonlinear simulation shows an immediate response to nozzle-area change. This re- 
sponse was not due to a startup transient offset. The linear model essentially starts at 
zero thrust. However, the linear model does have the potential of responding instan- 
taneously to a step change in an input (eq. (4)). For an output variable, the D matrix 
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represents the rate of change of that variable to an input. Since equation (4) is alge- 
braic, the output variable can respond instantaneously to an input change. Note that the 
state variables cannot change instantaneously since they are modeled by a differential 
equation (eq. (3)). Thus, at zero time the C X contribution is zero in equation (4). 

For a change in nozzle area the D(l, 2) element, which is multiplied by the change in 
nozzle area, gives the immediate thrust response. However, the D matrix elements 
are formed from a forced, steady-state match using the nonlinear simulation. Since 
the transient comparisons agree quite well (0.005 kN error), the D(l,2) element gives 
the correct sensitivity of the thrust to a nozzle-area change. Also, if the linear model 
had predicted the immediate thrust response, the whole transient would have been 
shifted up, since the system is linear, and the error would have been 0.155 kilonewton. 
Thus, the thrust response observed in the nonlinear simulation at zero time must be 
attributed to some effect not described in the linear- model approximation. From the 
figures shown the linear model does approximate the nonlinear system well at the sea- 
level -static -idle operating point. The same comparisons were made at the other five 
operating points. Although not presented, results of these comparisons showed that the 
linear-model approximations were as good as or better than the sea-level-static-idle 
case. 

These comparisons have shown that the nonlinear simulation has a very small 
startup transient whose effect is considered negligible. Experience has shown, how- 
ever, that a large startup transient can completely overshadow the effects of a small 
input. This can be especially troublesome when trying to run transients to compare 
with linear models. In general, gross transients requiring large inputs are not affected 
by the startup transient. Since the problem of startup transients in nonlinear simula- 
tions is of general interest, it is discussed in appendix F. 


Reduced- Order Models 

The results and methodology used to derive linear models from a nonlinear simula- 
tion have been presented. These are full-state linear models. Depending on the appli- 
cation of these models, it may be desirable to reduce the order of the model. For ex- 
ample, a problem results when trying to use full-state models for designing multivaria- 
ble controls that require full-state feedback. A full-state model could require a very 
complex control as well as a large amount of sensors to measure the states. Some of 
these states, however, may not be measurable. Thus, it is desirable to reduce the or- 
der of the model, keeping the states that are measurable and that contain significant dy- 
namic information for control purposes. 

In general, all model reduction techniques assume the system can be partitioned in- 
to two subsystems - one having fast dynamics (high-frequency eigenvalues) and one 
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having slow dynamics (low-frequency eigenvalues) 



The states having fast dynamics are approximated by their steady-state values by as- 
suming 

X F = 0 (19) 

In essence, this says that the high-frequency eigenvalues can be neglected and the dy- 
namics of the whole system can be adequately represented by the low-frequency eigen- 
values . If the eigenvalues of the system separate well, a good reduced-order model 
will result. Usually, a 10 to 1 separation of eigenvalues is acceptable, but in most 
physical systems the eigenvalues cannot be easily separated, and generally there is 
coupling between the states resulting in complex pairs of eigenvalues. Also, the iden- 
tification of states with specific eigenvalues is difficult. However, in general, for a 
particular application specific states must be included in the reduced-order model. 

Thus, the states that are retained specifically in the reduced-order model are based on 
intuition, engine controls experience, and engine protective requirements. The states 
that are not retained must be reconstructed from the reduced-order model from a linear 
combination of the retained states and outputs. 

Generally, a conflict arises when trying to determine how to reduce the full-state 
model. If all states are measurable, then the dynamic accuracy of the reduced-order 
can be assured by using the modal-reduction technique. In this technique the desired 
eigenvalues are specified, and the full-state model is transformed in such a manner 
that the states with eigenvalues that are not desired are approximated by their steady- 
state values . 

In most cases, however, all the states are not measurable. Therefore, the model 
must be reduced so that implementation of a control is possible. Yet the model must 
retain as much dynamic accuracy as possible . Once the states are specified, the 
normal -reduction technique is used. The full-state model is permuted and then parti- 
tioned into desired and undesired states 



Those states that are not retained (undesired) are approximated by 


( 20 ) 
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In general, problems can arise in using either method, depending upon the circum- 
stances. If the states to be retained are specified and the modal-reduction technique 
desired so that dynamic accuracy is retained, the eigenvalue- state relation must be de- 
termined. This can be difficult. If the states to be retained are not specified, a prob- 
lem similar to the one presented in using the modal technique also arises when using the 
normal-reduction technique. In general, many measurable states are available but not 
all can be used for control purposes. Also, one subset of states may provide a better 
control than another. Thus, an eigenvalue- state relation may have to be determined to 
pick the best subset of states. 

The choice of reduction method is not clear when the states to be retained for con- 
trol purposes are specified. For example, for the F100 engine, the states that must be 
retained are compressor speed, fan speed, burner pressure, and augmentor pressure 
from ref. 1). Thus, a fourth-order, reduced-order model is desired, and the use of 
both reduction techniques is explored. 

Modal reduction . - The modal-reduction technique requires that the eigenvalues of 
the retained states be specified. If the eigenvalue associated with each state can be 
identified, the reduced-order model will be dynamically accurate for the states that are 
retained since the reduced-order model will contain the specified eigenvalues exactly. 
However, it is difficult in practice to identify the eigenvalue associated with each state. 

One method used to determine the eigenvalue-state correspondence is to first form 
the eigenvectors associated with the system eigenvalues. The equation describing this 
is 

A e = Xe 

The eigenvector for an eigenvalue in a 16th-order system contains 16 elements; each 
element gives the weight of the corresponding state on the eigenvector. The eigenvec- 
tors can then be converted into mode shapes (magnitude and angle). Then the elements 
are ratioed to the largest element so that the relative contribution of each state can be 
determined. 

For example, for the sea-level-static -idle operating point the eigenvalues of the A 
matrix are given in table III, and the corresponding mode shapes in table IV. These were 
taken from appendix D. From the mode shape the state associated with each eigenvalue 
is determined. For example, if states 1, 2, 3, and 5 are desired, then for state 1 the 
vectors with the largest contribution of state 1 on each mode-shape magnitude (largest 
first real element) are selected. From table IV the vectors selected are 9, 14, and 15 
(where 14 and 15 are a complex pair). From table III these vectors correspond to a real 


eigenvalue at -2.86 radians per second and a complex pair at -0.744±0.20 radian per 
second. For state 2 the largest contributions of the second element are found in vectors 
14 and 15 from table IV. Note that this state does not contribute strongly to any mode 
shape. For state 3 the vectors 3, 14, and 15 are selected from table IV. (Although 
vector 3 contains a 1.00 for element 3, vector 14 contains a 0.9.) Note that complex 
pairs cannot be broken up. For state 5 vector 2 is selected. Table V (a) shows the state 
eigenvalue correspondence for a modal reduction at this operating point keeping states 
1, 2, 3, and 5. Although there are five different eigenvalues shown, a fourth-order sys- 
tem is desired. Coupling causes this difficulty: The complex pair cannot be broken up. 
Because the complex pair strongly influences three states, it should be chosen. The 
eigenvalue at -268. 16 radians per second is needed for state 5. The choice between the 
remaining two eigenvalues is arbitrary. Thus, two models are possible, and they are 
listed in table V(b). The only way to select one of the two models would be to try each 
and to compare with the full-state model by comparing similar transients . 

Normal reduction . - Another reduction method is to permute the states into two 
groups - those which are to be retained and those which are to be discarded. This is 
the normal -reduction method of references 7 and 8. The advantage of this method is 
its simplicity; the disadvantage is that the resulting model will not, in general, retain 
exact eigenvalues. For the sea-level-static-idle operating point, a normal reduction 
which retained states 1, 2, 3, and 5 was performed and the resulting eigenvalues were 
-0.878, -2.66, -142.8, and -245 radians per second. Comparing these eigenvalues 
with those selected for the modal reduction shows that essentially all three real roots 
were retained but that the complex roots were approximated by another real root. While 
it is true that in the normal-reduction technique some dynamic accuracy is lost for the 
retained states, it is much simpler to use. Thus, as long as the states that must be kept 
are specified, the normal-reduction method seems to be the better technique to apply. 

Fourth-order models were derived using the normal -reduction technique. The 
reduced-order matrices ABR, BBR, CBR, and DBR and eigenvalues are given for 
each of the six operating points in appendix G. 

Reduced-order model va lidation. - Once the reduced-order models are derived, 
they must be evaluated as to how well they approximate the dynamics of the full-state 
model. This is done by comparing transients run on both simulations. For example, 

3 percent steps in fuel flow and nozzle area were input to both linear models. Compar- 
isons of responses of selected states and outputs are shown in figures 7 and 8, respec- 
tively. The operating point is at sea-level-static idle. (A discussion of transient re- 
sponses of other states and outputs is given in appendix H.) 

Figure 7(a) shows the fan-speed responses to a + 3-percent fuel-flow step. The 
reduced-order model rises somewhat faster than the full-state model. This is due to 
the absence of high-frequency information in the model reduction. In general, the 
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transients have the same shape, and the error at 2 seconds is 2 rpm. The agreement 
should be good since fan speed is a retained state. Figure 7(b) shows the compressor- 
speed responses to a step in fuel flow. Compressor speed is a retained state and, 
again, agreement is good. The loss of some high-frequency dynamics model is indica- 
ted by the faster initial rise of the reduced-order linear model. The error at 2 seconds 
is 2 rpm. Figure 7(c) shows the compressor-diseharge-pressure responses. This is 
also a retained state in the reduced-order linear model. Here, again, the reduced-order 
model response rises more quickly, but, in general, the agreement is good. The error 
at 2 seconds is 0.25 kilopascal (kilonewtons per square meter) . Figure 7 (d) shows the 
interturbine-pressure responses. Interturbine pressure is not a retained state, but it 
can be reconstructed since it is linearly related to the retained states and the inputs. 

The immediate response at zero time for the reduced-order model results from the re- 
construction (the reason for this is discussed in ref. 8) . The transients, however, agree 
well, and the error at 2 seconds is 0. 1 kilopascal. Figure 7 (e) shows the augmentor- 
pressure response. This again is a retained state and again the loss of high-frequency 
information in the linear model is evidenced by the faster rise time. The transients 
agree quite well, with the error at 2 seconds being 0. 1 kilopascal. Figure 7(f) shows 
the thrust responses. Thrust is an output variable that can change instantaneously. At 
zero time the immediate response for the reduced-order model is due to the recon- 
struction process. At 0. 1 second the loss of high-frequency information is evident in 
the difference in rise times between the two responses. After 0.1 second the responses 
show the same shape, and the error is 0.005 kilonewton at 2 seconds. 

The loss of dynamic information is evident in the transient responses of the 
reduced-order model (fig. 7). If a state is retained, the reduced-order response rises 
faster than the full state. If it is a reconstructed state variable, there is an immediate 
response at zero time . Finally, if it is an output variable there is an immediate re- 
sponse at zero time and a faster rise time in the reduced-order response . 

Figure 8(a) shows the fan- speed responses to +3 percent in nozzle area. Here, the 
agreement is excellent with a 0.5-rpm error at 2 seconds. Figure 8(b) shows the 
compressor- speed responses. Here, also, the agreement is excellent, although the 
responses do diverge somewhat, with an error of 0.4 rpm at 2 seconds. Figure 8(c) 
shows the compressor-discharge-pressure responses. Like the fan and compressor 
speeds, compressor-discharge pressure is a retained state. However, the loss of 
high-frequency information is much more evident for this state. The full-state model 
shows a much larger undershoot and a much faster rise time than the reduced-order 
model response. After 0.1 second the responses agree quite well, and at 2 seconds the 
error is 0.1 kilopascal. Figure 8(d) shows the interturbine-pressure responses. The 
immediate response for the reduced-order model at time zero is due to the fact that 
this is a reconstructed state. The loss of high-frequency information is evident in the 
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lower undershoot and slower rise time for the reduced-order model response. The er- 
ror at 2 seconds is 0.2 kilopascal. The augmentor-pressure responses are shown in 
figure 8(e). This is a retained state, and the agreement between the two models is ex- 
cellent. Figure 8(f) shows the thrust responses. There is no obvious immediate re- 
sponse at zero time for this output variable, and the agreement is excellent. 

The loss of high-frequency information was not as noticeable for the nozzle-area 
transients as it was for the fuel-flow transients. The better agreement for the nozzle- 
area step is due to the negligible high-frequency sensitivity of the discarded states to 
nozzle area. The only state that seemed affected strongly by the reduction is the 
compressor-discharge pressure. The results probably would have been worse if aug- 
mentor pressure had not been retained. 

The reduced-order techniques used to generate the models discarded some of both 
high- and low-frequency states. Thus, both high- and low-frequency information was 
lost. Errors that show up in the reduced-order models can be attributed to this loss of 
information, especially when there is strong coupling between the retained and elimi- 
nated states. 

Although not shown herein, the full-state models and reduced-order models were 
compared at other operating points for response to transients in fuel flow and nozzle 
area. In all cases the results were as good as or better than those presented for the 
sea-level -static-idle operating point. These results indicate that a fourth-order model 
adequately represents this system. 


CONCLUSIONS 

This report presents the methods used to derive and validate linear models from a 
nonlinear digital simulation. To derive a good linear model at an operating point, good 
partial derivatives must be obtained. When using the finite-difference method to ap- 
proximate partial derivatives, small errors in the calculations can occur. These er- 
rors are bias errors, which must be minimized or eliminated if possible. This can 
accomplished by perturbing the states in both positive and negative directions about an 
operating point and averaging the resulting partials. For simulations in which differen- 
tial equations are solved explicitly, this method results in repeatability in the system 
eigenvalues for different perturbation sizes about the operating point. For simulations 
in which the differential equations are solved implicitly, exact repeatability of system 
eigenvalues for different perturbation sizes is difficult to obtain because of the iterative 
nature of the solutions. This is demonstrated by deriving linear models of the F100 
engine from an implicit, nonlinear simulation of the engine. Comparison of eigenvalues 
obtained by the averaging technique for different perturbation sizes about an operating 
point showed better repeatability of low-frequency eigenvalues, than of high-frequency 
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eigenvalues. For most applications (controls) the low-frequency eigenvalues are of 
greater interest. 

Using the newly calculated A and C matrices and assuming a forced steady- state 
match with the nonlinear simulation, the B and D matrices were calculated. The ac- 
curacy of these matrices was checked by showing that linearity holds for different per- 
turbation sizes in the controls in the nonlinear system. For the models in this report 
the B and D matrices were calculated by perturbing the inputs in both positive and 
negative directions about an operating point and averaging the resulting partial deriva- 
tives. Thus, for a given perturbation size the linear model predicts the average steady- 
state value that the nonlinear simulation would arrive at for the same perturbation size 
in the plus and minus direction about the operating point. 

Once the accuracy of the matrices was assured, the ability of the linear models to 
approximate the nonlinear system at an operating point was investigated. This was done 
by comparing similar transients run with the nonlinear and linear simulations. The in- 
puts were +3-percent steps in fuel flow and nozzle area. These comparisons showed 
that the linear models approximated the nonlinear simulation well and thus the method 
used to derive the linear models was good. The linear models were derived at six op- 
erating points . 

In comparing transients from the nonlinear and linear simulations, a startup tran- 
sient was observed in the nonlinear simulation. A startup transient can occur in a non- 
linear simulation as a result of a small error in the engine balance at zero time. In 
steady-state calculations an engine balance is forced. In transient calculations a differ- 
ent calculation path is made through the simulation and, depending on how quantities are 
calculated, small errors may result. With no other input applied, the errors tend to 
drive the system to a different operating point to achieve a balance. In the F100 non- 
linear simulation this error is considered to be negligible, but experience has shown 
that a large startup transient can completely overshadow the effects of a small input. In 
general, gross transients are not affected by startup transients. A method generally 
used to null out a small startup transient is to let the simulation run for a short time, 
with no input applied, to achieve an engine balance. 

The models derived from a nonlinear simulation are full-state models. For the 
F100 engine simulation the linear models are 16th order. Depending on the application 
of these linear models, it may be desirable to reduce the order of the model. For ex- 
ample, if the 16th-order models are used for multivariable control design studies that 
make use of optimal control theory, full-state feedback is required. For a 16th-order 
model this would result in an expensive and possibly impossible control to implement 
since it may be impossible to measure some of the states. Thus, the full-state models 
were reduced to fourth-order models. The states that were retained were found from 
previous studies to be measurable and capable of being used for engine control. The 
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states are fan and compressor speeds, burner pressure, and augmentor pressure. 

Two reduction methods were considered, modal and normal reduction. The modal- 
reduction technique requires selection of eigenvalues to retain the desired states in the 
reduced-order model. The eigenvalue- state relation is determined by mode-shape 
analysis. If this relation is well defined, the reduced-order models are dynamically 
accurate since the eigenvalues are retained. The problem encountered with this method 
is that it is, in general, difficult to determine the eigenvalue-state relation because of 
coupling between the states . And because of this coupling, the set of eigenvalues needed 
to describe the retained states may not be unique . The normal- reduction technique re- 
quires that the states to be retained be specified. The reduced model resulting from the 
use of this technique does not necessarily retain the eigenvalues associated with the re- 
tained states, and in this sense it is not as dynamically accurate as that generated by 
the modal method. Also, if the states to be retained are not specified, a mode shape 
would be required to determine which states to keep. However the retained states are 
specified, the normal- reduction technique was easy to implement and resulted in a 
unique reduced-order model. The eigenvalues resulting from this reduced-order model 
should be approximately equal to those that were shown to be important in the mode- 
shape analysis done for the modal-reduction technique. 

Using the normal -reduction technique, fourth-order models were generated at the 
six engine operating points. To determine how well the fourth-order models approxi- 
mated the 16th-order models, transients were run on both models for +3 percent steps 
in fuel flow and nozzle area. In general, the transients for the state variables from the 
two models agree quite well. The loss in high-frequency dynamics due to the reduction 
can be observed in the transient response. From the comparisons made at the various 
operating points, the fourth-order models approximate the full-state models well and 
can be used for control analysis. 

Finally, the methods used in this report to generate linear models from a nonlinear 
digital simulation are general. The application presented is for the F100 engine, but 
the method is not restricted to that simulation. The problems of bias errors in finite- 
difference approximations of partial derivatives and startup transients in nonlinear sim- 
ulations, again, are not limited to the F100 simulation. They are universal and must 
be dealt with whenever deriving linear models or comparing transients run with the non- 
linear and linear simulations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 1, 1978, 

505-05. 
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APPENDIX A 


ABR 

AN 

B 

BBR 

BLC 

C 

CBR 

CIW 

D 

DBR 
DP 


2. 5QP 

DP 2.5SC 

d/dt 
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FNMX 
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I 

K 

N. 


H 


N t 


n 

AP 

P T3 


SYMBOLS 

state matrix (nxn) 
reduced-order state matrix 
nozzle jet area, m 2 
control matrix (nxn) 
reduced-order control matrix 
customer-compressor-bleed- flow fraction 
output matrix (fxn) 
reduced-order output matrix 
fan-inlet-guide-vane position, deg 
feed- forward matrix (£xm) 
reduced-order feed-forward matrix 
empirical fan-exit AP/P parameter 

theoretical fan-exit AP/P parameter 

total derivative, sec - '*' 

eigenvector 

engine net thrust, kN 

function 

function 

identity matrix 

imaginary part of eigenvalue, rad/sec 
linear slope 

compressor speed, rpm 
fan speed, rpm 
index value 

change in pressure, kPa 
compressor-discharge pressure, kPa 



P T4.5 

P T7M 

PLA 

R 

RCVV 

SMAF 

SMHC 

s 

T T2.5C 

T T2.5H 



T T4HI 

T T4LO 

T T4. 5HI 

T T4. 5LO 

T T5 

T T6C 

T T7M 

t 

At 

U 

WFAN 

WFMB 

X 

AX 


interturbine pressure, kPa 

augmentor pressure, kPa 

power level angle 
ratio 

compressor variable vane position, deg 
fan stall margin 
compressor stall margin 
LaPlace operator, sec” 1 
duct temperature, K 

fan inside-diameter discharge temperature, K 
compressor-discharge temperature, K 
burner-exit temperature, K 
burner-exit fast-response temperature, K 
burner-exit slow- response temperature, K 
fan-turbine-inlet fast-response temperature, K 
fan-turbine- inlet slow- response temperature, K 
fan-turbine-exit temperature, K 
duct-exit temperature, K 
exhaust-nozzle-inlet temperature, K 
time, sec 

change in time, sec 
control vector (m) 
total engine airflow, kg/sec 
main burner fuel flow, kg/sec 
state vector (n) 
change in state 
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AX change in state derivative 
Y output vector (£) 

e error tolerance 

X eigenvalue, rad/sec 

a real part of eigenvalue, rad/sec 

Superscripts : 

total derivative, sec - ' 1 ' 

— vector quantity 

matrix inverse 
small perturbation 

t 

different quantity 
Subscripts: 
b bias 

c calculated 

D desired 

F fast 

i index 

j index 

n true value 

o operating point 

S slow 

sum summed value 

UD undesired 

u input 

y output 

+ positive perturbation 

negative perturbation 
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APPENDIX B 


LINEARIZATION MATHEMATICS 


The problem is to linearize nonlinear equations 



(Bl) 

(B2) 


where the subscript o indicates some operating point where the equations are balanced 
(ref. 10). Let X Q , U Q , and Y q be the operating point values of the states, controls, 
and outputs, respectively. Assume the system is operated close to nominal conditions 
so that X and U deviate only slightly from X Q and U Q . Then define U and X such 
that 


u = u - U Q 

X=X-X Q 


(B3) 


where U and X are small perturbations about the operating point. 

To linearize equation (Bl), expand in a Taylor series about X and U, and sub- 
stitute (B3) into (Bl) such that 

• 

X + X o =7(X o ,U o ,t) +A(X o ,U o ,t)x+ B(X o ,U o ,t)u + HOT (B4) 

Now, assume all the higher order terms (HOT) are small compared with the other 
terms and subtract (Bl) from (B4): 

X = A(t)X + B(t)U (B5) 

At an operating point the system is assumed to be time invariant; thus, 

X = AX + BU (B6) 

The A and B matrices are Jacobian matrices defined by 
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B ii 
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i = 1, n 
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II 

H- 1 
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9f i 
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II 

!— 1 
3 


(B7) 


Equation (B2) is linearized by again assuming small perturbations about the oper- 
ating point. Let 


Y = Y - Y 0 (B8) 

where Y is a small perturbation about the operating point. Then expand (B2) in a 
Taylor series and substitute (B3) and (B8) into (B2): 

Y - Y q = i (X Q , U Q> t) + C (X Q> U Q , t) X + D (X o , U 0 , t) U + HOT (B9) 

Again, assume that the HOT's are negligible, and subtract equation (B2): 


Y = C(t)X + D(t)U 


where, assuming time to be invariant, 


C(t) 



D(t) = D.. 
' ij 


8y i 

iH 

II 

a Xj 

j = l»n 

0y i 

i = 1 ,£ 

8Uj 

j = l,m 


(BIO) 


(Bll) 


Equations (B6), (B7), (BIO), and (Bll) comprise the linearized model. 

The approximation of partial derivatives of equations (B7) and (Bll) is accom- 
plished by finite differences in the nonlinear simulation. Digital perturbation logic is 
built into the simulation such that once the desired operating point is reached, the states 
are varied one at a time while the others are held constant. This is done to derive the 
A and C matrices. Once these matrices are derived, the logic in the simulation is set 
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up so that there is no steady-state error between the linear and the nonlinear models. 
This is accomplished by forcing 


AX + BU = 0 


Then 


X = -A -1 BU (B12) 

This equation is implemented by making small perturbations in the inputs and calculat- 
ing -A -1 B. Once -A _1 B is known and since A is known, 

B = -a(-A _1 b) (B13) 

For the last matrix (d) , using equation (BIO) yields 

Y = CX + DU 

Rearranging and inserting equation (B12) result in 

Y - C ^-A “ 1 BU ^ = DU (B14) 

where D is calculated by making small changes in the inputs while holding all the 
states constant such that equation (B14) is satisfied. 
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APPENDIX C 
THE PROBLEM OF BIAS 

In using finite differences to calculate partial derivatives, problems arise when 
there are small mathematical errors in the calculations. Since the perturbations made 
are small, changes in the model may also be small, and the resultant mathematical er- 
rors can cause significant errors in the partial derivatives. This is termed a bias 
problem. The discussion of this problem will proceed first with explicit calculations 
which, although not the case with the simulation here, are easy to understand. Then 
the effects of this problem on implicit calculations are considered. Finally, the prob- 
lem of selecting the correct perturbation size is addressed. Some of this material is 
presented in reference 11. 


Explicit Equations 

Explicit equations are those in which solutions are obtained without iteration. For 
example, given 


— =f(X) (Cl) 

dt 

this equation can be solved on the digital computer by using a difference equation such 
as 


X(n + 1) = X(n) + f [x(n)] At (C2) 

This is an open integration method known as Euler's method, wherein X(n + 1) can be 
calculated explicitly from previous values of X(n). When this method is used in the 
simulation, the finite-difference scheme is also explicit and is 


_d 

ax. 


i \ dt / 


f. - f. 

1 X.+d. 1 

3 1 


h. 


V d r x i 


(C3) 


where d. is the perturbation size. In steady state f. I =0 analytically, but mathe- 

’ b 

matically it may not because of small calculation errors . 

If there is a steady-state bias. 


X„ = X, + X„ 
c b n 


(C4) 
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• • • 

where X is the calculated derivative, X, is the bias value, and X„ is the true de- 
c b n 

rivative . 

For a linear system 


AX 

— = constant = K (C5) 

AX 


where K is independent of stepsize or direction of perturbation. Perturbing (C4) by 
AX and dividing by AX gives 


AX X, AX 
c _ b n 

AX AX AX 


(C6) 


Thus, the bias error is a function of the perturbation size and can get large as AX 
gets small. 

The bias problem for explicit calculations can be eliminated in two ways. First, 
the differential equations can be sampled in steady state, and the bias error can be 
calculated. Then, the method of finite differences can be used to determine the partials, 
and the bias can be subtracted. Second, the states can be perturbed in the positive and 
negative directions about the operating point. 

For a positive perturbation 



and for a negative perturbation 


(C7) 



(C8) 


Adding (C7) and (C8) gives 



But, for a linear system 
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Thus, 



2(^i^=2K and ^ = K 
\A X/ AX 


Hence, the bias problem is eliminated for explicit calculations by perturbing about 
an operating point. The actual implementation of this procedure on an explicit simula- 
tion is to make a positive perturbation and calculate an A matrix ^A + ) and then a nega- 
tive perturbation and calculate another A matrix ^A_ Then the matrices A + and A_ 
are averaged. 


Implicit Equations 

Implicit equations are those that require iteration for solution: 



This equation can be written as a difference equation of the form 

X(n + 1) = X(n) + f [x(n + 1)J At (Cll) 

Note that equation (Cll) cannot be solved explicitly because of the dependence of f on 
X(n +1). In the F100 simulation the steady-state part of the simulation uses a modified 
Newton-Raphson technique (SMITE) described in reference 11 to achieve an engine bal- 
ance. This is accomplished by specifying dX/dt = 0 in (Cl) and iterating to achieve 
the steady-state balance specified tolerance. 

Now this iteration capability is used to derive the A and C matrix elements. If 
state 1 is perturbed, a new state vector (x ! ) is formed, and equation (Cl) becomes 

f (x») - X' = 0 (C12) 

Since f (x 1 ) is known, the error variable can be specified by 


f(x') - f(x') c 


(Cl 3) 


= € 
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The iteration proceeds until equation (C13) is satisfied within some specified tolerance. 
In equation (C12) the nonzero derivatives result from the perturbation solutions, and 
they are used to calculate the required partial derivatives. 

The problem of bias with this type of solution results when the answer falls within 
the iteration tolerance, but the error is not zero. Then 


K=\ + K (Cl 4) 

where X^ is the error within the iteration tolerance. 

One obvious solution to this problem is to make the iteration tolerance very tight. 
However, this can have adverse effects, such as the program never achieving a solu- 
tion because of oscillations about the tolerance or the solution converging but requiring 
a large amount of computer time. 

Another possible solution is to calculate the error. This method is discussed in 
reference 7. 

A final method is to make positive and negative perturbations about the operating 
point. However, why this works is not clear. The solution of equation (C12) can con- 
verge within the specified tolerance. There is no guarantee on which side it will con- 
verge. The results in this report and in reference 12 (which discusses the same prob- 
lem) indicate that the positive-negative perturbation method works. A possible explana- 
tion for this (from ref. 12) is that the bias error tends to converge on opposite sides of 
zero for the different perturbation directions, causing the errors to cancel. Mathe- 
matically this is done by 




+ X. 


n 


and 


*c_ = X b_ + *n 

Adding equations (Cl 5) and (Cl 6) yields 

X + X = X. + X, + 2X 
c + c_ b + b_ n 

From experience we know that 

X b + X b -» X^ (almost zero) 


(Cl 5) 


(Cl 6) 


(C17) 
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Thus 


X, 


X„ = X + 
c n 


sum 


(C18) 


Note, however, that the iteration errors could possibly add, causing bias to worsen. 
Once X is calculated, divide by AX to get the partial derivative . 


Perturbation Size 

The perturbation size is restricted. It must be in the linear range about the opera- 
ting point, and it must be larger than the iteration tolerance for implicit solutions. The 
importance of the latter restriction is evident since no meaningful solution could be 
achieved otherwise. The former is important to generate the correct linear model. An 
investigation into the effects of iteration tolerance on perturbation size for implicit so- 
lutions was done in reference 11. 
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APPENDIX D 

F100 FULL-STATE LINEAR MODELS 

As an example of the linear model generation technique, linear models of the F100 
engine were developed at six operating points. The operating-point selection is based 
on engine operational constraints. The method of selection is described in reference 9. 
The table of operating points selected is given in reference 1. The points selected for 
illustration in this report are given in table VI. The A, B, C, and D matrices, the 
eigenvectors, and the matrix of mode shapes are presented for each operating point. 
Note that operating point 1 is the sea-level-static-idle operating point for which the ex- 
amples in this report are given. 


32 



Operating Point 1 (sea- level- static idle) 


THE A MATRIX 

-2.015 1.110 -1.39s 633.5 -1113. -1G.b7 -1.151 -.2920 

-.4462 -4.081 143.4 -252. 5 -124,0 27.70 -.409o -.8621 

1.0C6 5.056 -150. 9 7.055 279.3 -46.83 0.9030 1.311 

0. 1 830 0.1 994 103.6 -817 . C 606. 9 -2. 683 0. b7 16D-C 2- . 4589 

0. 1 C 27 -. 1488 1.797 34. 1j - 281.2 1 .400 - . o 2 4 3 D-0 1 -. 4 1 2 3D-0 1 

0.3548 -.2446 1.556 -2.510 102.2 -17. 6l 0.330o C.5350D-02 

0.2900 -.3917 2.533 -3.C9s 129.1 3.756 -19. s5 C.3250D-02 

-.3928 0.2305 37.58 2.665 -109.1 29. 6C -.3523 -20.00 

8145D-C2-.1792D-Q1-. 2991 -. 1 c27 -2.3s2 0.37J8 -. 1327 D-01 6.105 

-. 1 21 2D-0 3- . 26 55D-0 3- .4431 0-0 2-. 2891D-C2-.3439D-C 10. 5531D-C2- . 2 3 26 D-0 30 . 90 4 JD-0 1 
-4.146 -20.45 142.4 29.05 -1151. 133.0 -3.712 41.20 

0. 6298D-C1 0.2537 -210.7 635. 1 16.92 -2.923 G.6757P-C1 9.411 

0. 2795D-020. 1 1 27 D-0 1-9.36 1 26.23 0.7524 -.1299 0 . 2 9 05 D-0 20 . 4 1 8 3 

0. 7068D-C10.3277D-01-94.00 67.85 295. C -.4450 - . 2d 20U-Q 2 4 . 229 

-. 83 96 D-0 20. 382 3 D-0 1- • 279 1 -, 8c5jD- 01 3.459 -.3o30 1 9.93 0.7746D-01 

-2.597 3.636 -11.94 576.8 -20.16 -34.36 1.374 -.4735 

• 

0. 3121 0.3270 -1.166 2.66 1 2. 66C - . 2C C C 0-0 1 - . 47 65 0.2500D-02 

3.948 3.951 -1.355 0. 22CCD- 02- .5000D-030. 1760D-01 91 40 D-0 1-.5000D-04 

-. 8600D-020.1650D-01 1.537 C.0CC0 0 . 1 000 D-0 2- . 3 93 0 D-C 1 0 . 2 005 0.3000D-02 

-. 897 1 -.9058 1.740 G. 000C 0.0000 0 . 52 86 D-01 0 . 1 7 39 -.2928D-02 

0. 41 COD- 01 0.64 6 7D-0 1-. 1 43 5 0 . 1 06 9D- 0 20 . 5 5 2 j D-0 20 . 2C 4 4 0.20 3 3 0. 3315D-02 

-. 135 0 D-C 2-. 27 50 D-0 1C. 126 0 D-C 10. 3 GOOD- 04- . 2d3 OU-O 1 - . 1 49 5D-01 0.73 1 0 D-0 1 0 . 4 1 50 D -02 
-. 27 5CD-02- .435 0D-0 1C . 1 29C D-0 1 0. 5000D-04 0 . 5000 D-C 3-. 1 625D-01 0. 9420D-010. 1000D-03 
0.3450D-C20.3000D-02-.1200D-01-. 5CCCD-03- . 4 50 0D-C2 0 . 1 585D-01 - .81 60 D-0 1-. 1250D-02 
-50. CO 0.78C5D-02 39.70 1 9C C J-^5 0 . C COO - . ^2 C 2 D-02 0 . 2 1 7 OD-O 2-. 439 2D-02 

-. 6667 — . fc d 6 6 0.5 88 1 4. OOCO - . 5 OCOD-0 1 - .44 C6 D-04 0 . 23 9 0D-04 -. 5 85 9 D-04 

C. 36C0D-C1- .6000D-01-46. 15 0.5CCCD-C3D.30COD-G10.159S -.8595 -.2000D-01 

-10.23 -10.27 39.06 -50. Cl 0 . 2 CS5n-0 1 - . 4 3 95D-02 0 .23 90D-0 1 8 76 5D-02 

-.4544 -.4566 1.736 -2. 001 -1.959 1 760 D-0 3 0 . 1 0 44D-02- . 409 5D-03 

-4.596 -4.614 17.55 -1.420 -1.402 -19.98 C . 86 07D-01 -. 4 098 D-02 

0.345CD- 03 C. 5000D -04 0.131b D-0 1G.0C0C 0.CC00 -,2t39D-02-2C.01 0.0000 

0. 6543 0 .5850 -2.463 - . o5 COD- 02 0 . OOCO 24. sC 26.55 -50.00 

THE B MATRIX 

0.1656 20.66 -13.92 10.28 0.5b37E 05 

0. 3390E-01- .5900 -2.55c -21.0b 0.659 1 E C5 

-. 6665E-C 1-7. 935 5.827 42.42 -.3359E C5 

4033E-C1-9.546 1.538 2.206 -.1925E 05 

0.745 1 E-02-240. 3 1.061 -1.217 0.43UCE 05 

-.8139E-01-<4.59C 1.160 -2. 162 -. 1 33SE C5 

0.3450E-C 10.26C0 2.C34 -3 . 40 2 -.U'783 05 

-.3250E-C 2 1.230 -2.269 -13.20 J^.SO 

0.99b3E-C 1-.3o73 - . 5 28 5 E -0 1 -. 3 4 1 5 G.2333E 05 

0. 149 5E-C2- .73 578-0 2- .7 8 09 E -0 3 -. 5 4 4 2E- C 2 34 5. 1 


30.40 28.65 

-23.9s 

-174.8 

0 . 2 8 o 3 E 

Ob 


0.4394 12.90 

0.2822 

2 . 7 f- 7 

0. 16C7h 

0 5 


0. 1916 E-0 10.5433 

0. 1 257 E 

-01C. 12C5 

755.8 



0.1728 0.3210 

C.6164 

0.3532 

9616. 



0. 10 6 8 E-0 20 ,7 426 

-.1113 

C. 3539 

-655.7 



-.7700E-02-3.550 

-27. 19 

31.1s 

- . 4 155E 

C 5 


THE C MATRIX 






0. 16b2 -.2432 

2.949 

56. 1 C 

311.4 

2.273 

-.7570E-01-. 59925 

C. 1803 E-C 10.2394E- 

-02-. 1562E 

-0 10. 2525E- 

0 1 - . 9 9 9 1 

-. 23C0E-C : 

I-.3240I-02-.2500E 

0.0000 0. COOO 

0.0000 

G.OOCO 

o.ooco 

O.COOO 

O.COOO O.COOO 


-. 75 25 E-C 5 0.7 6 63 E-04-.5000E-0 30. 8C95o- 03- .32 C0E-C1-. 736 7S-03-. 1038 E-0 3 -.9 30 0 1.-0 6 
0. 1020E-C30.5349E-C3-. 1 43 5 E-0 1 -. 7 1 3 C 2-0 3 0 . 28 3 1 E-C 1 -. 42 7C 3-02 0 . S 1 40 E-C 4 C . 3 0 0 0L-0o 
0. 5 1 72 E-040 . 656 0 2-0 5- .3550 E-C * l O. 6 920*3-0 3- . 274 dE-C 1 - . 7 C 6 si -04 0 . 5 1 4 5 £-04 8000 E-0 b 

0. 1950E-C4- . 1 358E-050.89CGE-050. 1 4925-0 3- . 5856 E-0 20 . 458 3E-04 0 . 38 4 1 E-0 4- . 13C0E-G6 

0. 567 OE-C 10. 17 25 £-0 1 -. 2372 Zb 1 03-0* - . 6 050 E-C 20 . 3349 0.3206 -.4900 E-0 2 

0.29 OOE-OitO. 1C50E-03-. 1 1 1C F-C3C. 1 9° it-C-bO. 10003-040. 1 37SE-C3 - . 7425 E-C 3 4 50 OE-C 5 

0.0000 O.OCCO 1.0C0 0.0040 0.0000 G.OOOC O.COOO C.GC03 

0. 9 2C0E-0 6 0.3 4C0E-05-. 35 35 3-0 5- . 5 0CC.-C79.4 59 95-0fc0 . 4400E-05 - . 23 8 0 E-C4 lOOOE-Ot 
-.8b00E-C6-.5499E-C60.3185l-C5C.CC r C - . 5CC0 E-C t -. 3 550 E-C5 0 . 2 1 C5E-C4 C . 35CC t-06 

0. 76 5 0 E-0 60. 29 50 E-0 5-. 262 5E-C 5 C. OOCO -.90G0E-0C0.J65CE-Cs-.21 4s £-04 0. 1000 E-0 6 

C. 2200 E- 06 0.1150 2- 05 -.b7 50 £-06-. 4 99 9 E-0 » - . 5 00 0 E-C7 0 . 8 1 GC E-06 - . 4 5 sO h-05- . 4 99 9E -08 
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THE D MATRIX 

0. 1063E-C1 0.7800 1.747 -2.087 0.2061E05 

0. 3767 E-Q4-. 125QE-C20. 1684 C.2135E-01 126.2 

O.OOCC 0.0000 0.0000 0.0000 0.0000 

0. 3040 E-0 50. 5125E-0 3-. 12953-0 20.68 16E- 03 4. C2J 
3462E-C4- .3842E-C2 0.6C17E-C30.715 9E-02-7.20 2 
0.3485 E-050.346GE-030.5783E-040.7M 5E-04 3.798 
0. 6S55E-0o0.5350E-040.2018E-C3-. 72 0 0E-05 0. 8 109 


EIGENVALUES Of THE A MATRIX 
REAL PART IMAGINARY PART 


-852. 11279 
-268. 16113 
-135.01001 
-56. 317566 
-49.669357 
-48.432312 
-40. 37t 2 6 6 
-14. 1 19783 
-2. 8587856 
-21.319580 
-21.319580 
-19.7^7040 
-1.9813099 
-0.744 1G683 
-0.74410683 
-20.347610 


O.COOOOCOO 
0.00000000 
0.00000000 
0.00000000 
o.ocoooooo 
0.00000000 
C. 00000000 
0.00000000 
0.00000000 
0.63925683 
-0.83925668 
0.00C00000 
0.00000000 
0.20040676 
-0.20040876 
0.00000000 


MODIFIED EIGENVECTOR MATRIX 


-5.0913E-03 
6. 9963E-C4 
4. 2649E-03 

1.C000E 00 
-8.3483E-G2 

4. 2 6 38 E-0 4 

5. 35791-04 
-3. 1 4 77 E-04 

9.9511E-05 
1. 3289E-C5 
-1.7011 r-03 

- 1. 6391E- C2 
-6.5 3 55 E-0 3 

- 1. 1670 L -03 
9.33 6ttE- 07 

-1.9592E-02 


1. 1 117E-02 
4 . 6 97 2E-03 
-6 . 1 82 1 E- 0 1 
4 . 2 98 oE- 0 1 
1 . OCGOL 00 
-1.0936E-02 
-1.3 7 Co E-02 
1 . 985 8 E-0 2 
-1.56742-02 
-2 . Co6 5E-0 3 
7.3 0293-02 
-9.0C-46E-02 
-3.5817 E-02 
-4.3952E-C2 
5.6542E-04 
-2 . b 1 1 OE- 02 


- 1.0607E-C2 
-5.21 65E-C3 
1.0G00E OC 
6 • 6 7 72 E- C 1 
2.3190E-C1 
-8.0417E-C3 
-1. 1C38E-C2 
- 3 . 4 4 75 E - 02 
1.9765E-02 

2 . 5 53o E- C3 
-3. 3C 12 E-0 2 
9. 3643 E-C2 
3 » 6 fc 03 E-0 2 
4.62*81-0 2 
2.05C2E-03 
-1.3C79E-C1 


1.00CCE 00 

- 2. 4653E- 02 
-1.354 IE -0 1 

5.3303E-02 
1.4933E-01 
2. 0466 F -01 

2.C756E-01 
6. 0822E-02 

1. 1723E-01 
-1.5205E-C2 

1. 1 463L-0 1 
1. 7643 E - Cl 

- 1 . 7 4 S 3E-C 1 

2. 362 IE- 0 1 

2. 3 359E-01 

- 1. 1 37 1E-0 1 


-1 .S230E-G2 
-4.7425E-03 
2. 16043-01 
1 .2 95 0 E-0 1 
5.9244L-02 
-8. 32783-02 
-6 .4 36 3E-02 
1. 7573 E- 01 
1 . 05r OE -0 1 
1 . 1 3203^02 
4.7966E-C2 
5 .78C3E-02 
1.98315-02 
3 .979OE-01 
7 . 7737 L-0 1 
1.0000E 00 


3.7075&-C3 
«■ . 6 29d E-C 3 

5.9 131E-C2 
5 . o 5 77 - C 2 
5.677 

-1.O4S0E-C2 
-1 . 54 J*n-C2 
-4. 1 182 E-02 

- 1 . 6 1 37 E-0 1 
-1 .6695E-C2 

- 9. 6 700 l-C 2 
-1 . 4 CO 1 l-C 1 
-4.7 1*oj.-C2 

9 . 4 buG L-0 1 

- 2. o 232 a-C 1 
1. GOGOL CC 


- 1.1 18 62- 0 2 
-3. 2217E-03 
-4.4879E-02 
-4.7 34 12-02 
2.74 1 0E-Q 3 
5.7261E-04 
1 . 6 068 E-0 4 
5. 1 9673-03 
2.91033-01 
3.5792E-02 
-4. 0462E-02 
1 .00003 00 
3.78053-01 
1 .08803-01 
-b . 8 55 3 E-0 5 
-2.3 0503-01 

-2.911 IE-04 
-1.2 327E-04 
1 . 33C7E-02 
2. 28 13 L-0 2 

3 . 4 1 3 23 -0 2 
-3. 19725-03 
-1. 1 54 u e- 0 2 
1 .6 372L-02 
1. 63243-02 
1. 68583-03 
*.43G9c.-0J 

1 . 9 0 4 o E - C 2 
0.3214^-03 
I.OvuG'3 OC 
-9. 7 54 85 -01 
5.7108^-01 


-3.2395E-C3 
6.8762E-C5 
-2. 9384^-04 
5.837 0E-04 
1.0596E-G3 
- 7 . 56 94E-06 
1. 1 035 £-04 
-3.0293 E-05 
-4.6844E-03 
-5. 70 20 E-04 
2 • 5C 72 E-05 
2. 0542 £-01 
7.7036S-02 
8.8351E-03 
-8.3877E-05 

1.00005 00 

d • 93 bOE -02 
-3. 4 627 E- 03 
-1.49768-02 
3.27 37 e- 03 
1.27793-02 
1. 70483-02 
1.7666 E-02 
3 .05383-03 
9.6461 E-03 
-3. G 0613-03 
9 . 992 0E-03 
1. 54738-02 
1 . OOOOE 00 
1. 1 7 7c E- 02 
1 . 88562-02 
-1.781 1 £ -02 


-1 • 6743 E-0 2 
1. 7562E-03 
-4. 4450E-03 
d. S583E-03 
-1.6014E-05 
1 • 5424 E-03 
1 .43383-03 
-1. 3367E-C3 
-1 .2464E-01 
-1.511 6 E-0 2 
-4 . 09 94 E-0 3 
1. OOOOE 00 
3.7434E-01 
2. 7224 E-02 
-1 . 040 IE-03 
6.676 1E-0 1 

1 . OOOOE 00 
2.30303-01 

1.05 07E 00 
5. 1 7 64 E-0 1 
1 . 46423-0 1 
1 . 25 94 E-0 1 
9.691 OE- 02 
2 . 93 9 2E-0 1 
-2. 9361E-01 
7. 69 14E-C1 
-3.3961 £-0 1 
-4.5967E-01 
-5.651 OE-O 1 
-5.54 05 E~0 1 
1.0082E-C1 
-2. 7934E-01 


8.529 5 E-03 
- 1 . 5722E-02 
-8.4 954 E- 03 
-4.41273-02 
-1.C312E-04 
-3.8501 E-03 
-5.0783E-C3 
7 . 84 56 E- 03 
1. OOOOE 00 
1 . 1 882 E- 0 1 
2.0 840 E-0 1 
-1.7101 E-0 1 
-6 .31 70E-C2 
1 . 30 63 £-02 
4.891 2 E- 03 
-7.3419 E-02 

1. OOOOE 00 
1 . 795 6 E-0 1 
8.2109E-C1 
4 . 3737 E-0 1 
1 .4726E-01 
1.3981 E-0 1 
1.1 675E-0 1 
2.51253-01 
-2. 17 10 E-0 1 
- 5 . 5463 E-C 1 
-2 • 520 2 E-0 1 
-3.0140E-C1 
-4.6 728 E-0 1 
-3.7592E-01 
1 .21983-01 
-2.3 34 2 E-0 1 


3 * 9457 E-02 
-9.2711E-C3 
-1,7911 E-0 1 
-2.8508 E-02 
6 . 8297 E-C2 
1.0230E-01 
1 . 35 9b E-C 1 
1.324 1 E-C1 
3 . 33 94 E-0 1 
3 . 02 1 1 E-02 
2.4861 E-01 
3.7216E-01 
1.121 3 E-0 1 
9.8219 E-C 1 
4. 5685 E-0 1 
1.00003 00 

2 . 8 1 74 E- 03 
1 .9916E-03 
-6 . 4 950 E-C3 
1 . 1779E-02 
3.1015E-C2 

1 . 9 1 05 E-03 
1.67413-02 
-2.8498 E-02 

- 4 . 9524E-02 
-5 . 1 494 E-03 
-2.9331 E-02 
-3.5685E-02 
-1. 1953E-C2 

1 .OOOOE 00 

- 9 . 6650E- 0 1 
5. 2575 E-01 
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THE 

MATRIX OF * 

[ODE SHAPES, 

IN MAG. AND 

A N GL L( D LG . ) 

FORM 





1 

2 

3 

4 

5 

6 

7 

g 

1 

-0*50910-02 

0. 1 1120-01 

-O.IOolD-OI 

-0* 1119 D- 01 

-0 . 3240D- 02 

-0. 1674D-01 

0. 8529 D-02 

0.3946D-01 

2 

0 . 6996D-03 

C.46 97 D-C 2 

-0.5217D-02 

-0. 3222D-C2 

0 . 6376 D- 04 

0. 1756D-02 

-0. 1572D-01 

-0.927 ID-02 

3 

0. 4265D-C2 

-0.6182 

1.000 

-0.4488 D- 0 1 

-0. 2938D-03 

-0.4445D-02 

-0. 8495D-C2 

-0.1791 

4 

1,CD0 

0.4296 

0.6677 

-0. 4734 D- 01 

0 . 5 837 D-0 3 

0.5958D-02 

-0.4413 D-0 1 

-0.285 1D-0 1 

5 

-0. S3U6D-C1 

1.000 

C . 28 1 9 

0. 2741D-02 

0. 1C 6 CD- 0 2 

-0. 16C1D-C4 

-0. 1G31D-C3 

0.6830D-G1 

6 

0. 4264 D-C3 

-0* 1 C93D-0 1 

-0.8 04 2D -02 

0.5728D-C3 

-0 . 75 fc 9D- 05 

0. 1542D-02 

-0. 3850 D-0 2 

0.1023 

7 

0. 53 58 D- 03 

-0. 1371D-01 

-0. 1104 D- 01 

0. 1607 D-C 3 

0. 1 103D-03 

0. 14 34 D- 02 

“0.5078 D-02 

0. 1 36 C 

p 

-0. 3146 D- 03 

0. 1986D-01 

-C.3447D-01 

0.5197D-02 

-0. 3029D-04 

-0. 1337D-C2 

0.7846D-02 

0.1324 

9 

C. 995 ID-04 

-C. 1567D-C1 

0. 1 978D-C1 

0.2910 

-C.4684 D-02 

-0.1246 

1.000 

0.3339 

1C 

0* 1329D-C4 

-C. 205bl>-C2 

C.2554D-02 

0. 3579D-01 

-C. 5702D-03 

-0. 1512D-C1 

C. 1 166 

0.3021D-01 

11 

-0. 170 ID-02 

C. 7303D-C1 

-C. 3301D-01 

-0.4046D-01 

0 . 2507 D -04 

-0.4099D-02 

0.2084 

0.2466 

12 

-0. 1639D-01 

-C. 9C05D-C1 

C. 9364 D-C 1 

1.000 

0. 2054 

I.OCO 

-0.1710 

0.3722 

13 

-0.6535D-02 

-0. 356 2D- 0 1 

C . 36 6 CD- 0 1 

0.3781 

C.7704D-01 

0.3743 

-0.6317D-01 

0.1121 

14 

-0. 1167D-02 

-0.M95D-C1 

0.4630D-01 

0.1088 

0 . 88 35 D-02 

0.27 22 D-C 1 

0. 1 306 D-0 1 

0.9822 

15 

0. 9 337 D- 06 

C. j654D-03 

C.205CD-02 

-3. 6655 D- 04 

-0.8 368 D- 04 

-0. 1040D-02 

0 . 4 6 9 1 D-02 

0.4569 

16 

-0. 1959D-C1 

-C. 251 ID-01 

- 0 . 1 3 06 

-0.5305 

1 .000 

0 . 6876 

-0.7 342 D-0 1 

1.000 


9 

10 

1 1 

12 

13 

14 

15 

1 6 

1 

1 .000 

C. 11C8B-01 

-121.3 

-0.2911D-C3 

0. 6936D-01 

1.00C 

o.ocoo 

0.281 7D-02 

2 

-0. 2465D-01 

C. 463bD-02 

-90.69 

-0. 1233D-03 

-0.3463D-C2 

0.2C65 

7.058 

0. 1992D-C2 

3 

-C. 1354 

0. 1564 

29.69 

0. 1331D-01 

-C. 1498D-01 

0.9429 

6.993 

-0.6495 D-0 2 

U 

C.5330D-01 

0. 9S&3D-C 1 

21.40 

0.2282D-C1 

0.32 74 D-02 

0.4792 

4. 804 

0. 1176D-C1 

c 

C. 1493 

C. 5802D-0 1 

1.2 19 

0.341 3D-01 

0. 12780-01 

0.1468 

-0. 1646 

0.3102D-0 1 

6 

0 . 2047 

0. 4 6 24D-0 1 

-149 . 6 

-C. 3 1 97 D-0 2 

0. 17C5D-01 

0.1331 

-2.988 

0.191 ID-02 

1 

0.2C76 

0. 466 ID-0 1 

- 146.5 

-0.11 t4D-0 1 

0. 1767D-01 

0.1073 

-5.304 

0.167 4D- 0 1 

6 

0. t C62D-C1 

C. 1276 

58. 19 

C. 1 6 37 D-0 1 

0. 3054D-G2 

0.2734 

4.476 

- 0. 285 CD-0 1 

9 

0. 1172 

0. 1364 

101. 6 

0. 1632D-01 

0. 9846D-02 

0.2582 

-171.5 

-0.4952D-01 

10 

-0. 1521D-C1 

C. 1438D-01 

101.2 

C . 1c 86 D-C 2 

-0. 30 06 D- 02 

0.6820 

80. 10 

-0. 5 149D-C2 

1 1 

0. 1146 

C. 7760D-0 1 

109 . 1 

G . 94 3 1 D-C 2 

0. 5992D-C2 

0.2990 

-171.6 

-C . 2933D-C 1 

12 

0. 17 6 4 

0. 107 1 

112.6 

0. 19C5D-01 

0.1547D-01 

0. 3087 

-166.3 

-0.3569D-01 

13 

-0. 1 75 C 

0. 3622D-0 1 

112.1 

0.5 3 2 ID-0 2 

I.OCO 

0,5165 

-174.6 

-0. 1 195D-01 

14 

C. 2362 

0.7247 

-22. 7t 

1 .COO 

0. 1 178D-01 

0.4734 

-169.2 

1.C0C 

15 

0.2336 

C.5801 

63. 65 

-0.8755 

0 . 1886D- Cl 

0.1119 

-5.426 

-0.9665 

16 

-0. 1137 

1.000 

C.OCOG 

0.5711 

-0. 1791D-01 

0.2574 

-174.9 

0.5258 
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Operating Point 2 


THE A MATRIX 


-4.307 0.1764 5.058 *05.0 -725.7 -2.044 1.124 -1.059 

-.4275 -5.600 128.4 -232.7 -428.4 27.01 2.242 -2.511 

1.036 6.217 -164.2 -4. 970 1054. -82.41 -5.456 5.084 

0.5753 - . 6 005D-0 1 132.7 -573.6 113.7 -6.603 -.5635 -2. 118 

0.1228 -.1764 1.191 2C.98 -115.9 2.606 0.1715 -.1161 

0.8345 -.1380D-010.4965D-01-.2b44 36.42 -19.83 -.2078 0.1230D-02 

0. 6669 -.7323D-020. 1071 -.752C 86. b4 0. 1185 -20.42 0. 1295D-02 

-.1022 0.8039 16.21 C.4725 -106.9 29.42 0.5370 -20.01 

-.8544D-02-.1859D-010. 1720 1 205D-0 1 -8 .7 1 4 0.6922 0.2361D-C1 6.626 

-. 126 6D-C3- .2751 D-030. 2548 D- 02- . 1607D-C3- . 1289 0. 1025D-01 0 .35 68 D-0 30 . 98 1 5D-0 1 

-1.276 -6.996 23.37 6.14C -1296. 101.5 6.712 36.37 

28t0D— 0 1-. 4633 -52.94 199.2 -28.16 2.224 0.1464 9.757 

-. 12 6 8D-02- .2 059D-0 1-2.353 8.855 -1.251 0. 9889 D-0 1 0 . 64 85D-0 20 . 4 337 

-.1645 -.2527 -24.18 22.28 142.9 1.604 0.1127 4.479 

1225D-010.3097D-01-.1032 -.1445D-01 5.759 -.4474 19.82 0.1262 

-1.637 1.887 -2.940 115.6 347.2 -27.11 -1.802 -.4661 


I. 035 1.014 -4.126 9.511 9.576 -.6806 -.3843 0.0000 

II. 38 11.34 -4.040 0. 5 0C0D-04-. 1 10 OD-O 1 0. 5C0CD-04- . 437 1 0.5C00D-04 

-. 2260D-0 10.7500D-02 5.618 -. 1 500D-020 . 1 000D-C 1-. 1 500D-02 1.059 0.0000 

-2.916 -2.931 5.890 2 450D-02- . 3 30 0 D-0 1 0 . 540 0 D-02 0 . 1 0 65 0.1000D-03 

0.1076 0.1008 -.4231 0. 5000D- 030 . 1960D-C 10. 5293 0.3272 0.C000 

244 0D- 02- .5 000D-04 0.4185D-0 2-. 1 005 D-0 2- . 2250 D-0 2-. 5000 D-05 0 . 5 1 8 5D-0 1- . 50 00D-05 
-. 30 00 D-0 20* 5C 0 OD-O 4 0 . 185 1D-0 1-. 745 0D- 03 0 .500CD-0 4-. 2 60 0D-03 0. 17 24 0. 100 0D-04 

0. 2070D-02 C.5 000D-04-. Ill CD-010. 1 O0GD-O4O. 9 050 D-0 20 . 000 0 -. 1044 0.0000 

-5 0.0 0 0. 5 OOOD-05 39.47 623 5D- 03- . 5 605D-0 20 . 67 1 0D-03 - . 3 1 2 0D-C20 . 0000 


-. 6667 - . 6 b67 

0.5847 

- • 8 3 1 0 D ■ 

-0 5-. 7 450 D-0 40. 89C CD-05 -.6 245 D-040. 0000 

0.30 10 D-0 10. 1750 D- 

-01-47.87 

0. 1500D 

-020. COOO -.25C0D-02-1 .3 03 0.0000 

-9.636 -9.609 

38.50 

-50.00 

7 29 5D-0 10. 26 75D-C2- . 32 79 D-0 10. 500 0D-05 

-.4283 -.4272 

1.7 1 1 

-2. OOC 

-2.0C3 0. 107 CD-03 15 00 D-0 20. 0000 

-4.425 -4.413 

17.58 

-3.119 

-3.168 -19.78 -.2936D-0 10.5000D-05 

1885 D-030. 0000 

0. 198CD- 

•010. 9970D- 

- 03- .5 OOOD-05 C. 1 06 9D-02 -2 0.00 -.5 000 D-0 6 

0.5852 0.5910 

-2.361 

0.5SC0D- 

-030 • 4450D-0 1 30.73 20.16 -50.00 

THE B MATRIX 




3864E-01-200.5 

-111.0 

0.2690 

-247.0 

0. 1987E-0 1-322 . 0 

-6.044 

-65.59 

-6326 . 

0.1351 818.3 

14.75 

62.83 

0.9318E C 5 

0.4313E-01 79.69 

2. 120 

6. 348 

-3859. 

0. 1 1 95E-03- 1232 - 

8.749 

-1.96 6 

-57.55 

0.4 406 E-0 2—3. 520 

9.238 

-.1921 

13.60 

0.811 5E- 02 155.8 

9.147 

1.001 

-184.0 

-.45612-02^86.26 

-1.478 

10.30 

-205. 9 

0. 453 7E- 01-5.923 

-. 1244 

-.5119 

9997. 

0.681 1E-03-.6385E- 

•01-. 2060E- 

■02-.7579E- 

02 148.0 

6.326 -1001. 

-18.57 

-74.85 

1175E 06 

0.1461 -17.52 

-.4516 

-1.648 

0.1040E 05 

0. 6497E-G2- .7313 

2030E- 

■01-.7318B- 

01' 461.6 

0. 86 76 E-0 1-29.35 

0. 1 166 

-1.211 

4243. 

- . 353 8E-03 4.107 

-.8741 

0.3384 

-7.715 

0. 29 45E-0 1 248.4 

-96.13 

21.04 

-1698. 


THE C MATRIX 

0.4734 -.6987 5.018 94.47 22.43 10.37 0.6890 -.5099 

0. 138 2 E-0 10. 3137 E- 050. 0000 0. 5GC6E- 07- . 977 5E-02- . 534 5E-04 0 .53 15 E-040. 499 9E-08 

0.0000 0.0000 0.0000 C.OCCC 0.0000 0.0000 0.0000 0.0000 

0. 74 18E-040.5101E-05-. 1454E-040. 7 1 3CS- 04- . 1 524E-0 1- . 7396E-04 0 . 7 922 E-04-. 1950E-06 
0. 1 655E-04 0 .1287E-03-. 2578E-02-. 7990E-04 0 . 1667E-0 1-. 1 0 75E-02- . 87 1 3E-0 40 . 2450E-06 
0. 154 6E-04-.2 981 E-0 4 0. 1 01 3E-0 30. 7 54 0£- 04- . 1 599E-0 1 0 . 432 4E-03 0 .2 1 28 E-03-. 3300 E -06 
0. 91 1 5E- C50.4 1 28L-05- . 1 320E-04 C. 1 848E-04- . 4 14 5E-0 20. 4306E-04 0 .76 97E-0 4-. 9000E-07 

0.4777 0.484G -1.908 79503-03-. 16 103-C 1 3.37C 2.034 0.1C01E-04 

0.0000 0.0000 0.0000 0.0000 0.0000 0.4999E-080.4999E-060.4999E-08 

0.0000 0.0000 1.000 0.0000 C.0000 0.0000 0.0000 0.CG00 

0. 3 550 E-Co 0.50 00 E-0 7-. 15 30E-050. 350 0E- 070. 1 00CE-06 0. 9 997E-C8-. 1522 E-04 0. COOO 
-. 375 CE- 050.3 OOOE-06C. 17 09 F.-05-. 9 9 99E-07-. 4 000 E-Ob -.4000 E-07 0. 16 91 E-04-. 4 59 9E-08 
0. 3950 E-06 0.00 00 1570E-0 5-. 3 500E-07 0 . 1850E-05 0. 0000 - . 1 5 77 E-040 . 0 00 0 

0. 122 OE- 06 0.0000 - . 4380E-06 0. 3 500 E-080 . 8 0 00E-06 0 . OC OC - . 4 049 E-050 . 0 00 0 
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THE D MATRIX 

0.4435E-C3-235.0 34.51 

0.6343E-040. 1778 C.7403 

0.0000 0.0000 0.0000 
0. 3327E-G5-.4 1 44E-0 2- . 3 1 04E-0 20. o 22 6E- 04 C . 1156E-C1 
0. 1 469 E- 060. 1278 E- 010. 2367E -030. 2431E-02-. 13 26 E-0 1 


-7.816 -626.0 

0.204 23 -03C.6590E-02 
0.00C0 0.0000 


64 90 E-0 7-. 1307 E-0 10. 1 4463-02- . 3 29 8E-0 30 . 64b 52-02 
0. 2120E-07-.3548E-020.6496E-03C.4 o372-C40. 2525L-C2 


EIGENVALUES OF THE A MATRIX 
REAL PART IMAGINARY PART 


-562.04395 
-248.34012 
-60. 906693 
-50.085373 
-49. 964691 
-43. 535172 
-30.572876 
-30.572376 
-21 . 942566 
-21. 942566 
-18.756088 
-12*951914 
-7. 0543 7 0 S 
-2.8662462 
-0.67332358 
-1. 9776640 


O.OOCOOOOO 
0.00000000 
O.OOOOCOOO 
0.00000000 
0.00000000 
0.00000000 
10.531186 
- 10. 531186 
1.5567312 
- 1 . 55 o 7 3 1 2 
0.00000000 
0.00000000 
0.00000000 
o.ooooocoo 

C.CGCOOOOO 

O.OOCOOOOO 


MODIFIED EIGENVECTOR MATRIX 
-7.7 125E-03 -2.2732E-02 


2. C83U-03 
3. 63056-02 
1.00C0E CO 

- 1. 1296E-01 
6. 6690E-04 
1.2905E-03 

-1.017 4E-03 
3.7914E-04 
5.C154E-03 

-4.4 9 24E-0 3 

-1. 8975E-C2 

-7.5490E-03 

- 1. 122 JE-03 
2. 1 b 51 E- 05 

-1. 39 36E-02 


-6 .6 176E-C3 
7.1 345 £-C 1 
1 . 0000 E 00 
-4.5479E-01 
5.400 IE-03 
1 .0590E-C2 
-1.9 113E-02 
1 . 223 1 E- 02 

1 . 6 C87E-03 
-5. 3133E-02 
1 .2 909E-03 
5.C731E-C4 

2 . 4 36 6 E- 0 2 
-7 .7 308E-C5 
-2. 1 5b 3 E- 02 


-2.2661E-C2 
-1. 6533E-C2 
-1.23221' - Cl 
-2. 1 333E-01 
-3. 1290E-02 
7.76 413-03 

6 . 4 Q 1 9 E-03 

9.6 6 27 E-03 
3. 7172 E- 01 

4 . 5 9 98 E-0 2 
-8.6973i.-02 

1 . 0 C 0 0 E CC 
3. 7953 3-d 
1. 9753 L -Cl 
-3. 760 4 E-03 
-5.9316S-C 1 


1.599 jL-O 3 
-4. 8733E-C3 
-1. 1 29 2 E-0 1 
- 1.0536L-C1 
-*.t795E-02 
4. 2 8 4 oF- 02 
- 1. 6749E-02 
-5. 3182E-02 
6.3076E-02 
6. 4540E-03 
4. 6882E-02 
5.617GE-02 
1 . 8 5 66e-0 2 
-4. 39363-01 
1 . 0000 L CO 
-3. 2998E-01 


-2. 3255Z-C2 
-1.641 CE-C2 
1 .2773E-01 
1 . 1 802E-01 
7.116 IE-02 
-2.UC69E-C2 
-1 . 8 1 35E-0 1 

1 . 7 09 9 E-0 1 
1 .297bE-01 
1. 3923E-C2 
6.2859E-02 
1 .1794E-01 
4.038 IE-02 
-1 .4275E-01 
1 .00001 00 
2.419CE-C1 


- 1 .68162-03 
2.9654E-C3 
1. 95633-03 
1 . 16o0E-03 
-1 .09333-02 
-3. 3 C bl £- C3 
6 .35452-02 
-1.4oC6x.-02 
-2.47 15 E-0 2 
-2.504 5 E- 0 3 
-1. 5532S-C2 
-2. 3763L-02 
-7 . 7923E-C3 
-3.2 384 E-C1 
1 . COCOS oc 
-5.4076S-C2 


&.92133-04 
1.8o7lE-06 
1 . C472E-04 
1 .7506 E -04 
-4. 1612E-04 
-2.26253-04 
- 1 . 4 1 2 1 E -04 
7.50203-05 
-4.374CE-C5 
-5. 9259E-06 
-6.22273-06 
-1.7361E-C2 
-*.74172-0 3 
-1 .99633-03 
? . 3 247E-05 

1.00003 00 


-2.^0853-04 
2. 5127E-0o 
-4. 59 0 9 E-0 5 
-b.33053-05 
1.69033-04 
9.71663-05 
5. 96 56 E-05 
-3. 337 1 E- 05 
-4 . 042 0^-05 
-4.91733-06 

2.7 606E-C6 
7. 513 9 E-03 
2.81 69E-G j 
8.2534 E -04 
-4. 1581^-05 

1 .00003 00 


4 .4432 E— 0 2 
-5.983CL-02 
5. 77863-02 
-2. 97682-02 
5. 94943-02 
-2.7 0^5 E-0 2 
-2. 9 86 7 E-0 2 
3.4121E-C2 
1.0000E 00 
1 . 1986E-01 
1 . 37733-0 1 
-6. 1744L-0 1 
-3. 03 273-01 
-2.30^73-02 
2. 4728E-02 
-1 .67 33 E-0 1 

1 .00CG2 00 
6. 86 18E-01 
9 .55463-0 1 
8.5072E-C1 
7. 4 134E-0 1 
7. 32893-0 1 
7. 1032E-01 
9 . 8767 E-0 1 
1 .2814 E-0 1 
-1 .6517E-02 
4. 66233-02 
1. 1BC03-01 
-1 . 1 7 2b E-0 1 

5. 6 3 Ob E-02 
8. 33C5E-01 
2.8523 E-0 1 


4. 3715E-C3 
-1.76043-02 
3 • 6566 E-0 1 
2.9329 E-C1 
3. 420 6 E- 01 
-9.2 389 E-02 
- 1 . b 59 8 E-0 1 
2 . 9 627 E-02 
4.67333-01 
6 . 0 7 92 E-02 
-5.0 266 E-02 
4. 3144E-02 
2.9483E-C2 
5.2665 E-0 1 
1 .4903 E- 01 
1.00003 00 

3 .066b E-0 2 
1.159C3-01 
1.4967 E-C 1 
1 . 1 0 5 3 E-0 1 
4. 5231 E-C2 
5. 0 088 E-02 
4. 6786 E-02 
1.03C5E-01 
-1. 1502E-C2 
1.0GC0E 00 
-2.C849E-02 
-5.0782 E-02 
-6 . 6253 E-02 
-7. 2463 E-02 
5. C^bOE-02 
-1.20903-02 


4.9862E-C2 
8.20873-02 
2.5335E-C1 

2 . 4 6 64 E-C 1 
2.9895E-01 
-1*7315 E-02 
-2. *1963-03 
-7. 1908E-C2 
-1.0765E 00 
-1 .2222E-01 
-5 . 54 04 E-C 1 
-1.C-159E CO 
-3.6816E-C1 
7.30203-02 
-1.43723-01 
1 . 0 0 00 E CO 

5.9C16E-C2 
2.7171E-03 
2.7633E-G3 
7. 1 283E-C3 
3.5441 3-C 2 
4 . 0620 E-C2 
3 . 804 3 E-C 2 

2 . 7 7 44 E- 02 
1.5496E-02 
-4.7589E-C3 
1.38173-02 
1 .62813-02 
1.C000E 00 
3.23523-03 
4. 1 121E-02 
-7.08263-03 


0 • *+ 9333- 02 
3.95123-02 
8. 149*2-02 
8. 34672-02 

3 . 6 o 52 3- C 1 
2. 1 59CE-01 
3. 3573 t -01 
3.101 5E- 0 1 
2.4b2l£-0l 
2. 1265E-02 
1 .7214Z-01 
2.56903 -0 1 
7. 52043-02 
7 . 06162-01 
9.4764E-0 1 

1 .OOOCE 00 


6.24773-01 
-1. 0125 E- 01 
-9.9954F-02 
-1.9077F-02 
6 . 440 33-0 1 
6.463 85-01 
6. 62333-01 
4.874CE-01 
4. 3 1 S 4 E-0 1 
1 .8576E-02 
3.61323-01 
4. 67 97 E-0 1 
7.9*243-02 
b. 4154E-01 
1. GOODS OC 
6.96443-01 
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THE 

; MATRIX OF 

MODE SHAPES, 

I ti MAG. AND 

ANGLE (DEG. } 

FORM 





1 

2 

3 

4 

5 

6 

7 

8 

1 

-0.7713D-02 

-C. 227JD-01 

-0. 2266D-01 

0.6921 D-03 

- C. 29 09D-03 

0. 4443D-Q1 

0. 3539D-0 1 

-39.99 

2 

0 • 2683D-02 

-0.661SD-C2 

-C. 1653D-01 

0.1867D-05 

0 • 25 13D-05 

-0 . 5983D-0 1 

0. 5936D-0 1 

-57.10 

3 

0*3630 D- 01 

C.7135 

-C.1232 

0. 1047D-03 

-0 • 459 ID- 04 

0.5779D-01 

0.3146 

10.28 

4 

1.000 

1.000 

-C.21 33 

0. 1751D-03 

-0.6331D-04 

-0. 2977D-G 1 

0.2711 

4.915 

5 

-0* 1130 

- C. 4548 

-0.31 29D-0 1 

-0. 4161 D- 03 

0. 1690D-03 

0.5949D-01 

0.3212 

3.847 

6 

0.6669D-03 

0.5400D-02 

C.7764D-02 

- 0 . 228 2D- 03 

0.97 17D-04 

-0.2704D-01 

0.6647 D-01 

-145.6 

7 

0. 129 1 D- 02 

0. 1059D-01 

C.8492D-02 

-0. 1412D-03 

0. 5966D-04 

-0.2987D-01 

0.1174 

-136.0 

8 

*0. 10 17 D-0 2 

-0. 1911 D— 0 1 

C . ^ 66 3D- 02 

0.7502D-04 

-0 . 3337D-04 

0. 3412D-01 

0. 5499D-0 1 

112. 6 

9 

0 . 37 9 1 D- 03 

C. 1228D-01 

0.3717 

-0.4874D-04 

-0. 4042D- 04 

1.000 

0.8298 

111.5 

10 

0.5C15D-04 

0. 1609D-02 

0.460 OD-C 1 

-0. 5926D-05 

-0.4917D-05 

0.1199 

0. 9652D-01 

108.6 

11 

-0.4492D-02 

-0.5314D-01 

-0.8897D-01 

- 0 . 6223D-C5 

0.276 ID-05 

0.1377 

0.3934 

140.2 

12 

-0. 1898D-01 

0. 1291D-02 

1.000 

-0. 1798D-0 1 

0. 7514D-02 

-0.8174 

0.7190 

132. 6 

13 

-0.7549 D- 02 

0. 507 3D-03 

C . 37 95 

-0.6742D-02 

0 . 28 17D-02 

-0.3033 

0.2612 

130.4 

14 

-0. 1122D-02 

0.2437D-C1 

0.1975 

-0. 1S98D-02 

0 . 8253D- 03 

-0.2305D-01 

0.3760 

37.11 

15 

0. 2185 D- 04 

-0*773 1D-0& 

-C.3760D-02 

0.9825D-04 

-0.4158D-04 

0. 2473D- 0 1 

0.1489 

89.94 

16 

-0. 1394D-01 

-0. 2158D-01 

-0.5932 

1.000 

1.000 

-0.1673 

1.000 

0.0000 


9 

10 

11 

12 

13 

14 

15 

16 

1 

0. 1648D-01 

131.1 

-0.1 o 82D- 0 2 

0.649 3D- C 1 

0.6248 

1.000 

0 . 8067D-0 1 

0.5902D-01 

2 

0. 1210D-01 

151.5 

0.2S65D-02 

-0. 3951 D- 01 

-0.1012 

0.6882 

0.1159 

0.2717D-02 

3 

0.1206 

-86.48 

0.1958D-02 

0. 8150 D- 01 

-0.9995D-01 

0.9555 

0.1499 

0.2763D-02 

4 

0. 1119 

-86.76 

0. 1186 D- 02 

0.8847D-01 

-0.1903D-01 

0.8507 

0. 1 105 

0.7128 D- 02 

5 

0.8495D-01 

-98.68 

-0.1094D-01 

0. 3665 

0.6440 

0.7413 

0.4523D-01 

0. 3544D-0 1 

6 

0.3475 D-01 

74.32 

-0.3008D-02 

0.2159 

0.6464 

0.7329 

0. 5 009 D-0 1 

0. 4062D-0 1 

7 

0. 1288 

140.3 

0 • 6354D-G 1 

0.3357 

C. 6623 

0.7103 

0* 4679D-0 1 

0. 3804D-Q 1 

8 

0. 1266 

-62.28 

-0. 1481D-01 

0.3102 

0. 4874 

0.9877 

0.1031 

0.2774D-01 

9 

0. 1020 

-19.06 

-0.24 71D- 01 

0.2462 

0.4319 

0.1281 

-0. 1 150D-01 

0. 1550D-01 

10 

0.1095D-01 

-20.14 

-0.2S04D-C2 

0 . 2 1 27 D-C 1 

0. 1858D-01 

-0. 1652D-01 

1.000 

-0.4759D-02 

11 

0 . 5545D-0 1 

-8.283 

-0.1 553D- C 1 

0.1721 

0.3613 

0. 48 82 D-01 

-0. 2Q85D-0 1 

0.1 382B-G1 

12 

0.9237D-01 

-19.53 

-C.2375D-0 1 

0.2589 

0.4880 

0.1180 

-0. 5078D-0 1 

0. 1828D-01 

13 

0.3143D- 01 

-20.31 

-0 .77 92D-C2 

0. 7520D-01 

0. 7942D-0 1 

-0.1173 

-0. 6625D-01 

1.000 

14 

0.3267 

-153.0 

-0.3238 

0.7062 

0.8415 

0.5681 D-0 1 

-0.7 246 D-0 1 

0.3235D-02 

15 

1.000 

0. 0000 

1.00C 

0.9476 

1.000 

0.8330 

0. 5046D-01 

0.411 2D-0 1 

16 

0.2893 

-98.75 

-0.54 C6D-0 1 

1.000 

0.6964 

0.2852 

-C. 1 209 D-0 1 

-0.7083D-02 
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I 


Operating Point 3 


THE A MATRIX 


-7.118 

1.866 

0.9745 

397. 1 

-563. 2 

-21.13 

-3.848 

-1.481 

-2.835 

-7.777 

115.6 

-233.6 

-70.35 

25.26 

-3.8 12 

-3.762 

8.490 

7.988 

-163.7 

-.2650 

211.8 

-83.04 

11 .38 

7.282 

1.400 

-.3630 

133.3 

-572.1 

23.2b 

-7.678 

1.041 

-2.993 

0. 1203 

- .2073 

1.181 

21.80 

-105. 4 

2.468 

-.37 38 

-.1777 

0. 6972 

1376 

0.3516 

C.OGCO 

6.419 

-18.41 

0.3450 

0.0000 

0. 5281 

-.3273 

0.8162 

1 6 50 D-0 1 

14. 28 

3.792 

-19.25 

0 . 500 0D-05 

-.7211 

0.8261 

13.64 

0.123CD-01 

-17.97 

26.03 

-.9668 

-20.00 

3884D- 

01-.3956D- 

■02 0. 9495D- 

-0 1-. 5394D-01 

-.8606 

0. 3953 

-.54 94 D- 

■01 6.064 


-. 5766D-03-.5765D-040. 140 1D-0 2-. 8220D-03- . 1 275D-0 1 0 . 5852D-02 - . 824 1 D-0 30. 8 98 4D-0 1 
-6.1C5 -5.132 13.22 0.2385 -152.3 55. 7C -8.183 38.14 

-.8092D-C 1 - . 3t>4 1 -34.88 128.7 -2.238 0.7838 -.9891 0-01 9.079 

-. 3594D-C2- .16190-01-1 .550 5.719 - . 954 3D-0 1 0. 3 484 D-0 1 - . 4 34 CD-0 20 . 403 5 

-.1780 -.2163 -15.82 15.31 103.2 0.9455 -.1252 4.154 

0. 3153D-020.138SD-01-.3637D-010.0000 0.4131 -.1592 19.92 0.8839D-01 

-. 8383 0.9513 -.3653 68.67 30.24 -10.86 1.626 -.4284 

1.550 1.679 -6.161 15.06 15.01 -1.229 -.2620 0.0000 

17.82 17.71 -6.320 C. 7500D-020 . 0000 1 200 D-C 1 - . 25 4 0 0.2000D-02 

-.400 0 D-C 2-. 8000D- 01 6.779 1 5C CD- 02 0 . 3 COOD-G 1 0. 46 00D-C1 0.77 35 -.5 000D-03 

-4.679 -4.890 9,353 0. 2 GOOD-O 2- . 3 500D-0 1 0 . 1 500 D-C 1 0 . 7 1 OOD-O 1 0 . COO 0 

0.1745 0.1900 -.7062 0. 1 850D- 02- . 5000D-02 1.237 1.087 C.1850D-02 

C. 5000D-050.6750D-020.0000 -.9950D-Q3O.GOOO 0 . 1 88 5D-02 0 . 22 98 D-0 1 0 . C000 

0.000 0 -.67 5 OD-O 20 . 1 7 75D-0 2-. 2 450D- 03 0 . 5 00 0D-04 0 . 484 5D-02 0 . 454 0 D-0 1 0 . 0000 

-.5000D-050.6b50D-02-. 1295D-02-. 50COD-050.0000 5 380 D-02- . 6 5 1 2D-0 1 - . 1 000D-04 

-5C. 00 - . 4216 D-0 2 39.76 0.0C00 O.OCOO - . 201 7D-02 0 .4 1 02D-02-. 5000D-06 

-.6667 -.6666 C.589C O.OOGC 0.0000 35 85 D-C4 0 . 7293D-0 40 . 0 000 

-.50 00 D-0 40. 54 50 D-0 1-4 8. 02 3 1 50D- 02- . 2b 5 OD-O 1 - . 3 495D-01 - . 55 3 7 -.5000D-04 

-9.804 -9.737 38.94 -50.00 11 20 D-0 1-. 4030 D-02 0.41 10D-0 2- . 3530D-02 

-.4357 -.4327 1.731 -2.001 -2.C01 - . 1 6 1 0D-03 0 . 2 1 95 D-0 3 -. 1 4 1 0 D-03 

-4.482 -4.453 17.80 -3.254 -3.243 -19.75 0 . 1 6 50D-02- . 1 05 5D - 02 

O.OCOO 0.0000 0. 13410-010.0000 0 .5 OOOD- C5 0 . 0 000 -20.00 0.5000D-06 

0.5574 0.5740 -2.230 0. 1 0 C OD-O 30 . 5 0 OOD-O 3 2b. 54 23.94 -50.00 

THE B MATRIX 

-. 9729E-01-2S9.6 -232.8 34.54 -263.5 

2765 E-0 2-3. 460 -94.6 9 -40.11 -3370. 

0.1787 36.58 283.6 13C.2 0.3810E05 

0.47S3E-01-3.310 2o.23 11.58 - 1975. 

0. 27 9 7E-02- 16 83 . 13.05 -3. 948 -57. 15 

0. 1244 E-0 1 47.03 4.230 -2.671 18.85 

0. 582 1 E- 02 7.420 18.04 -6.08^ 1G.70 

-.4372L-02-4.640 -24.06 -12.02 41.55 

0.2657E-01 C.6600E-C1-1 .283 -.6352 2224. 

0.3S9 0E-C30. 102 0 E-0 2-. 19 19 E-0 1-. 9 63 2E-02 32.88 
4.026 -12.75 -2C4.1 -94.11 -.2787E05 

0.8976E-C 1-3. 403 -2.659 -1.222 2522. 

0. 4C12E-G2- .1501 -.1189 -.54a0E-01 112.0 

0.554 IE-01-2.712 -1.083 -1.4S9 1059 . 

1342E-C4- .76C5 -. 1646 0.264b -3.620 

0. 1325E-C 1-44.38 -64.71 17.04 -327.5 


THE C MATRIX 

-.7961 -1.036 5.86G 106.9 117.2 12.38 -1.730 -.8501 

0. 3 3 37 E-0 1 0. 5 96 C E-0 3-. 1 495E- C20 . OOOC - . 28 30E -0 1 -. 69 24E-02 1 5 1 OE-O 20 . GOO 0 

C.0000 0.0000 0.0000 O.OCOO 0.0000 O.OOCO 0.0000 0.0000 

G. 5234 E- 04 0. 3697 E-04-. 9496E-G 40. 1 8 ODE- 05- . 17 37E-C 2-. 4301 £-03- . 9341E-C40. 3550E-07 
0.8548E-040.8350E-04-. 16 1 8 E-0 2-. 2 44 OE- 05 0 . 2 1 32E-G2- . 4fc 64 £-03 C . 1 1 46 E-G 3-. 10 50 E-0 7 
0. 20C8E-C40. 1128E-04-.2831E-0 40.29CCE-05-. 3694E-0 2-. 1 31 5E-03 0. 10 39E-C 30. 4COOE-07 
0. 1 28 2E-050 . 104 6E-0 4-. 27 06E-04C. 9 00 0£-0t- .9298E-0 3-.542 IE-04 0.52 07 E-C4 -. 75 00 E-08 

0.8784 0.9275 -3.44b 1 300E-02 0 . 1 6b 0 E-0 1 6.051 5.333 -. 2500E-02 

4959 E-0 3- . 1 000E-06-.4 999E-08-. 9 9 97 E-08 IC-OOE-Ot-. 5C 00 £-07- . 1 0 94 E-C 30. 0000 
O.OCOO 0.0000 1.000 0.0000 0.0000 0.0000 0.0000 0.0000 

0. 11 50E-C70 .5 900E-06-. 1225E-C60. 3 COOL- 07- . 4999E-C8-. 379CE-06- .6354E-G5-. 4998E-09 
31 50E-G7-.815 0E-06G. 1520E-06C.4 996S-090. 3350E-0t 0 . 4 8 95E-06 0 .77 80 E-0 5 -. 1650E-07 
0. 31 50E-070.1 4 05L-05-. 29452-060. 60C0E-C7-. 1C 00 E-0 7-. 1062E-CS- . 1 3 79E-04- . 5000E-08 
0. 200 OE-O 80. 4 750 E-Ob-. 67 CG E-0 7-. 3 OOOE- 080. COCO - . 2 99 0S-06 - . 33 05 E-05 - . 4 9 98 £-09 
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THE D MATRIX 

-.5910E-02-222 .3 -10.5b -19.33 -313.5 

0.2134E-C3 1.006 1.797 C. 1 C76E-01-2. 71 1 

0.0000 0.0000 0.C00C O.COOO 0.0000 

0. 8332 E- 060. 5758E-G2-. 29 24E-C 20. 6 26 0E-03- . 4040E-0 2 
0. 38 88 E-05 0. 2b3GE- 030. 2863 E- 020*28 C4E-02-. 66G5E-02 
C. 3785E-06C.5352E-02-.5275E-030. 1 6 22E-0 30 . 447 G E-02 
0.3030E-070. 1017E-020.3733L-03C. 1566E-03- .1G25E-02 


EIGENVALUES OF THE A MATRIX 


REAL PART 

IMAGINARY PART 

-570. 10107 

0.00000000 

-197.71367 

0.00000000 

-83. 236343 

0.00000000 

-61. 849380 

0.00000000 

-49.990326 

0.00000000 

-44. 194748 

o.ooococco 

-33. 174606 

0.00000000 

-9.0144463 

3.8102236 

-9. 0 1 44463 

- 3 . 8 1C 2236 

-23.922791 

4.4413338 

-23.922791 

-4.4413338 

-5.2559175 

O.COOCOOCO 

-IS. 180176 

0.00000000 

-0.68158114 

0. COCOCOCO 

-20.93^601 

0.00000000 

- 2 . OC 5784 0 

o.ooooocoo 


MODIFIED EIGENVECTOR MATRIX 


-1. 1233E-02 
4. 7579E-03 
8. 7737E-03 
1. OOOOE CO 
-1. 1567E-0 1 
2. 7634E-04 

3.8 040E-04 
-3. 8844E-04 
9.1 7 30E-05 
1 . 2 1 6 1 E-05 
- 1.C258E-03 
-1. 9684^-02 
-7. 6305E-03 
-1 . 66 OCE- 03 
-7. 8 547 E-Ob 
-1. 4926E-02 


-4.4336E-02 
-8.9331E-03 
7.6875E-01 

1 .CCOQE 00 
-7 . 0 22 3 E- 0 1 
2.6539E-03 
2.9153E-C3 
-2.5517E-02 
1 . 075 4E-02 
1.4027E-03 
-3.313 QE- C2 
6 .545b E- 04 
2.5991E-C4 
3 . 5b6 5 E-0 2 
-4. 8 26 4 E-05 
r 5. 6702 E- 02 


1 . 8 556 E-C 2 
3.57 77 E- 03 
3.7650E-01 

4 . 6 786 E -C 1 
1.0000E 00 
-1 . 05o5E-G2 
-2.0178E-02 

- 1. 567JE-C2 
-1.21 73E-0 1 

- 1.53S0E-C2 
8.9517E-C2 

-2.2830E-C1 

-8.7971E-02 

-1.0104i,-01 

5.9332E-03 

-2.9100E-02 


-1. 4427E-01 
-1* 3081E-01 
-2. 1557E-01 
-1. 8068E-01 
3. 7642E- 0 1 
2. 0787 E -01 
3.6143E-01 

2 . 6523E-C1 
4.C012E-0 1 

3. 3049E-C2 
3. 2490E-01 
4. 3131 E-C 1 
1.2660E-01 
8. 8670E-0 1 
1.000CE 00 
1.064ciE 00 


3.3814E-02 
1 . 088 IE-04 
2.0574E-01 
1 . 0696 E- C 1 
4.451 1E-0 1 
-1 .2290E-C1 
-7 • 5 929E-02 
6 .2 96 7 E-0 2 
1. 1 3 1 2 E-0 1 
1.5965E-02 
-1 . 6 224 E-02 
2. 672 IE-02 
1 .65o3E-02 
7.0633L-C1 
2. 00G2E-0 1 
1 .OCOOE CO 


9. 1 fa 9 9 E - 0 2 
1.2146E-G1 

1 . 9 392 E-C 1 
1.8163E-01 
4. 5244E-01 
-3 .7fc 83 E-C2 
- 1 . 0 378E- 02 
-2.C304E-01 
-7.84443-01 
-8. 5307E-C2 
-4 . 57 30 E-C 1 
-6 . fat b9E-0 1 
-2 . 3842E-C 1 
5.7621E-C1 
-1.7lD2t.-C1 

1 .CCOOF cc 


-3. 1513E-02 
-2.7854E-02 

- 1 . 4 3 1 2E— 0 1 
-2.4209F-01 
-1.5636S- 02 

5. 59 12E-03 
3.1365E-C3 
1 . 6 05o E-02 

3 . 6 384 E-0 1 
4.7567E-02 
-1.0010E-01 
1. OOOOE 00 

3.8 009E-0 1 
2.0469E-01 

-1.3264E-03 
-5.07 1 5 E-C 1 

9 . 4 3 b C E -0 1 
2 • 9 C83E-C 1 
1.00C0E 00 
9.526SS-C1 
9.5324E-01 
o. 1444E-01 
b • C248E- 0 1 

6.9 342t.-0 1 

- 2. 4 293E-02 
-1*3 1 63E-04 
-7.9636E-C2 
-2.2416E-02 
-6.9221E-04 
-8.31 9 IE-02 

o.l 775JJ-C1 

2.5 317E-01 


-1.2759E-04 
1 .151 9E-05 
-2.8009E-05 
- 1 . 0 459 E-05 

7.9 372 E-04 
27 56E-05 
- 1 . 0307E-0b 
- 9 . 53 65E-06 
-1.8C46E-04 
-2. 1 954E-05 
6.4331E-07 
2.8674E-03 
1.0737E-03 

2. 31 48 E-04 
1 .0189E-07 
1 .COOOE 00 

-4. 1088E-03 
4.061 4E-03 
1. 3254E-02 
1.14993-02 
1 . 3368E-02 
-3. 6761E-03 
4.1 344S-02 
-4. 1 262 E-0 3 
-2. 20871.-02 
-2.2596E-03 
-1. 4455 E-02 
-2. 0752E-02 
-6.8635E-03 
-4.1 3893-01 
1 .OCOOE 00 
4.2861 E-02 


-8.8030E-02 
8.9281 E-02 
-B.2266E-02 
-2. 3555E-03 
-1 . 1 4 99 E-0 1 
3.7422E-02 
4. 0475E-02 
-4 .8931E-02 
-9.3122B-01 
-1 . 1184E-01 
-1. 1361E-01 
1 .OOCOE 00 
3.71 78E-0 1 
4.74 14E-02 
-3 .29523-02 
2.9625E-01 

9. 1783E-02 
9. 1 188E-02 
1.6488E-01 
1. 2846E-01 
6 . 24 65E-0 2 
3.6351E-02 
2 • 2648E-02 
7.6402E-02 
-2.5811 E-02 
1. OOOOE 00 
-3.3727E-02 
-6 .5657 E-02 
-8. 5635E-02 
-8. 83793-02 
2 • 4545E-02 
-2.8J12E-0 2 


-3.4469E-02 
-1.1673 E-0 1 
-1. 1282E-01 
-1 . 4015E-01 
-2.0753E-01 
1.1551 E-02 
-8.9416E-03 
6 . 6439 E-02 
1. OOOOE 00 
1.1560 E-0 1 
3.5930E-01 
5.3165 E-0 1 

1.9 228 E-0 1 
-1. 7231 E- 01 
1.37 20 E-02 
-5.2436 E-0 1 

-4.2327E-03 
-7.C472E-03 
-9 . 8844 E-05 
5.32 70 E-04 
2.0661 E-02 
2. 2611E-03 
-4 . 7 28 2 E-02 
2. G185E-02 
3.7801 E-02 
3.9618 E-0 3 
2 . 209 1 E-02 
3.4782 E-02 
1 • 1 714E-C2 
-3 . 8 300 E-0 1 
1. OOOOE 00 

6 . 7 4 97 E- 0 2 


7 . 5 859E-0 1 
-1.2369 E-0 1 
3. 3592 E- 01 
3. 9400E-01 
8. 1150 E- 01 
6 . 8204E-0 1 
7. 4680 E- Cl 

4.D816E-01 
9.9085E-02 
-2.9048E-C3 
8.4279E-02 
1.5179 E-01 
9.8809E-O3 
9.8003E-02 
1 .OOOOE 00 
2.0657E-01 

5.1 845E-02 
-9.2160E-03 
2.3 134E-02 
2. 6319 E-C 2 
3.511 IE-02 
2. 9050E-02 
3. 0759 E-02 
1 .0 J16E-C2 
-5.5420E-03 
1.6512 E-03 
-6.521 0E-03 
-4 . 4 273 £-03 
1. OOOOE 00 
-2.7 052 E-02 
3.374 1 E-C 2 
- 1 . 53 1 7E-02 
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THE 

MATRIX OF MODE SHAPES, 

IN MAG. AND 

ANGLE (DEG. ) 


1 

2 

3 

4 

1 

-0. 1123D-01 

-0.4434D-01 

0.1856 D-0 1 

-0. 3151 D- 01 

2 

0.4758D-02 

-0.8933D-02 

0 . 35 78 D-0 2 

-0.2785D-C1 

3 

0.8774D-02 

0.7687 

0.3765 

-0.1431 

4 

1.000 

1.000 

0.4679 

-0.2421 

5 

-0. 1157 

-0.7023 

1.000 

-0. 1564D-C1 

6 

0 . 27 6 3D- 03 

0. 26 94D-02 

-0.1 0 57 D-0 1 

0. 5591 D- 02 

7 

0 * 3804 D-03 

0.2915D-02 

-0.2018D-01 

0 . 3 1 36 D- 02 

8 

- 0 . 3884D-03 

-0. 255 2D-0 1 

-0. 15 67 D- 01 

0. 1606D-C1 

9 

0 . 9 1 7 3D- 04 

0. 1075D-01 

-0.1218 

0.3838 

10 

0, 1216 D-04 

0. 1 40 3D- 02 

-0.15390-01 

0.4757D-01 

11 

-0. 1026D-02 

-0. 33 1 3D- 0 1 

0.8952D-01 

-0. 1001 

12 

-0, 1968D-01 

0. 6546D-C 3 

-0.2263 

1.000 

13 

- 0 . 7830D- 02 

0. 2599D-03 

-0. 87 97 D- 01 

0.3801 

14 

-0. 1660D-02 

0.356 6 D-0 1 

-0.1018 

0.2047 

15 

-0.7855D-05 

-0.4826 D-04 

0.5933D-02 

-0. 1326D-02 

16 

-0. 1493D-01 

-0.5 67 OD-0 1 

-0.29 1 C D-0 1 

-0.5071 


9 

1C 

11 

12 

1 

55.77 

0. 6 9 24D-0 1 

-24. 80 

0.9485 

2 

178.4 

0. 35 88D-0 1 

-44. 95 

0. 2909 

3 

77.69 

0. 1999 

1.693 

1.000 

4 

69.64 

0.1843 

C.82SC 

0.9527 

5 

20. 12 

0.4486 

-C.4679 

0.9582 

6 

28. 05 

C • 9 094 D-0 1 

-152.1 

0.6144 

7 

19. 17 

0. 5 6 9SD-0 1 

- 142.4 

0.6C25 

8 

14.93 

C. 1503 

117.8 

0.6964 

9 

-31.09 

0. 5604 

126.8 

-0.2429D-C1 

10 

-50.02 

0.61 37D-C 1 

124. 

-0. 1816D-03 

11 

-30. 46 

0. 3C9 4 

137. 1 

-0.7964D-C1 

12 

-25.61 

0.4859 

132.8 

-0.2242D-G1 

13 

-40. 54 

0. 1690 

131. C 

-C.6922D-03 

14 

-38.69 

0. 6449 

5. 815 

-C.8319D-01 

15 

O.OCOC 

0.1863 

85.61 

0.8173 

16 

-34.02 

1.000 

o.coro 

0.2532 


FORM 


5 

6 

7 

8 

-0. 1276D-03 

-0.8803D-01 

- 0. 3447D-0 1 

0.5460 

0. 1152D-04 

0.8928D-01 

-0. 1167 

0.1273 

-0.2801D-04 

-0 . 8227 D-0 1 

-0.1128 

0.2822 

-0. 1046D-04 

- 0 . 2356D- 02 

-0.1401 

0.3065 

0 .79 37 D-03 

-0.1150 

-0.2075 

0.6325 

0. 2276D-04 

0 . 37 4 2D- 01 

0.1155D-01 

0.5042 

-0. 1031D-05 

0.4048D-01 

-0.8942D-02 

0.5867 

- 0 . 9536D-05 

-0.4893D-01 

0. 6644D-01 

0.3743 

-0. 1 8 OSD- 03 

-0.9312 

1.000 

0.2915 

-0.21 95D -04 

-0.1118 

0. 1156 

0.2346D-01 

0.6433D-06 

-0.1136 

0.3593 

0.2373 

0.2867D-02 

I.OCO 

0.5317 

0.3233 

0.10 74D-02 

0.3718 

0.1923 

0. 8979D- 0 1 

0. 2315D-03 

0 .474 ID- 0 1 

-0.1723 

0.6308 

0. 10190-06 

-0. 3295 D-0 1 

0. 1372D-01 

1.000 

1.000 

0.2963 

-0.5244 

0.7668 

13 

14 

15 

16 

-0.41C9D-02 

0.9178D-01 

-0.4 233D-02 

0.5185D-01 

0 . 40 6 1 D- 02 

0.9119D-01 

-0.7047D-02 

-0.921 6D-02 

0.1325D-01 

0.1649 

-0.98b 4 D-04 

0.2313D-01 

0. 1150D-01 

0.1285 

0.5327D-03 

0.2632D-01 

0. 1339D-01 

0.6246D-01 

0. 206 6 D-0 1 

0.351 ID-01 

-0. 36 76 D- 02 

0.3635D-01 

G.2261D-02 

0.2905D-01 

0.4134D-01 

0 . 22 65 D-0 1 

-0.4728D-01 

0. 3076D-01 

-0. 4126D-02 

0.7640D-C1 

0.2018D-01 

0. 1032D-01 

-0.220 9D-01 

-0.2581 D-0 1 

0 . 37 80 D-0 1 

-0.5542D-02 

-0. 2260D-02 

1.000 

0 . 39620-02 

0. 165 ID-02 

-0. 1445D-01 

-0. 3373D-01 

0. 22 09 D-0 1 

-0.652 ID-02 

-C.2075D-01 

-0.6566D-01 

0. 3478D-01 

-0. 4427 D- 02 

-0.6363D-C2 

-0 ,eS64D-01 

0. 1 171D-01 

1 .OOC 

-0.4139 

-0 . 88 38 D-C 1 

-0. 3o30 

-0.2705D-01 

1. OOC 

0. 2455D-C1 

1 . OCO 

0. 3374D-0 1 

0. 4 2 66 D -01 

-0. 2831D-01 

0. 675GD-01 

-0. 1532D-01 
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I 



Operating Point 4 


THE A MATRIX 

-1.171 0.5732D-01 5. 183 416.7 -813.2 -.7133 0.4008 -.3028 

-.1259 -1.502 130.4 -238.2 -424.5 7.532 0.6992 -.7102 

0.2984 1.642 -159.5 -5.241 1034. -23.32 -1.652 1.445 

0.1538 -.3350D-02 129.6 -559.6 118.0 -2.592 -.1830 -.5880 

0.3205D-01-.4310D-01 1.120 21.34 - 110.4 0.6823 0 .4 963 D-0 1 334 1 D-0 1 

0.8100 -.1 984D-010.2541 -1.292 160.5 -19.72 -.2917 0.1410D-02 

0.6266 -.6680D-020.795GD-01-2.265 412.0 0.7880D-01 -20.49 -.4240D-02 

-.1078 0.7927 57.93 1.965 -391.9 30.44 0.5966 -20.01 

-.9830D-02-. 1767 D- 010. 5607 -. 4 675D- 0 1- 3 1 . 9 6 C.7294 0.4744D-01 6.992 

-. 1460D-C3-.2618D-C30.8471D-C2-. 6230D-03-. 4739 0. 1082D-01 0.7028 D-0 30. 1036 

-1.446 -7.291 85.01 25.60 -5015. 113.0 7.993 35.69 

-.3360D-01-.5103 -195.7 734.9 -113.3 2.543 0.1687 10.23 

-. 1493D-02-.2268D-01-8.698 32.66 -5.033 0.1130 0 .75 18 D-0 20 . 4545 

-.1764 -.2732 -89.45 70.34 551.8 1.776 0.1139 4.698 

-. 1324D-0 10.3146D-01-. 3588 -.7485D-G1 21.57 -.4884 19.81 0.1307 

-1.826 1.995 -10.37 438.1 1342. -30.62 -2.177 -.4800 

0.2768 0.26C0 -1.128 2.679 2.579 -.2007 - . 83 55D-0 1 0 . 5000D-04 

3.08 1 3.083 -1. 105 0 . 1 330D- 020 . 1 1 95D-0 1 0 . 5000D-05 - . 1 1 15 0.0000 

-. 7950D-02- .9000D-02 1.484 -.1 35GD-02-. 2 400 D-0 1 0. 2800 D-02 0. 27 16 -.5000D-04 

-. 7597 -.76 25 1.559 0. 1 3 COD- 02- .5 000D-0 30. 1350D-02 0.28 OOD-O 1-. 500 0D-04 

0. 2914D-01 0.2397D-01-. 1170 0. 1400D-030 . 1 600D-C 20 . 1 97 8 0.1224 0.5000D-Q5 

2590D-02-. 1000D-030. 86 75 D-0 20. 0000 - . 7 15 0D-02 0 . 8400D-03 0 . 71 12D-01-. 1000D-04 

0. 1900D-030.2500D-010.8160D-02-. 5000D-050.2850D-01-. 2835D-02 0. 1 1 36 0.0000 

0. 3 195 D-0 2 0. 1000D-0 3-. 1 17 4 D-0 10. 500 CD- 05 0 . 1 9 1 CD-0 1 - . 2255D-02 - . 96 90D-0 1 0 . 1000D-04 

-50. CO 0. 5 C OOD-O 5 39.27 6 6 00D-03- . 59 35 D-0 20 . 21 2 1 D-02- . 1 74 3D-0 1 0. 1 00 0D-05 

-.6666 -.6667 0.5818 -. 87 95D- 05- . 7 9 1 5D-C40 . 2826D-04 - . 2556 D-030. 0000 

0. 3890D-010. 3250 D- 01-47. 72 0. 7 3 50 D-0 20 . 9 550 D-0 1 2040D-01 - 1 . 3 25 O.OCOO 

-9.405 -9.419 38.26 -50.00 - .5 900D-020 .71 00D-03- .4356 D-0 10. 0000 

-.4180 -.4187 1.700 -2.0C0 -2.000 0 . 28 50 D-04 - . 1 952D-020. 0000 

-4.323 -4.336 17.59 -3.192 -3.090 -19.76 - . 44 59 D-0 1 0. 5000 D-05 

3 965 D-0 3- .5000D-050. 2042D-C 1-. 50C0D-06-. 5000 D-0 50 . 5650 D- 03- 20. 00 1000D-05 

0. 5985 0.6006 -2.424 0. 5 GO 0D- 04 - . 4 1 5 OD-0 1 30.85 20.40 -50.00 

THE B MATRIX 

-.411 0 E-G 1-61.57 -30.09 0. 1 855 - 186.0 

0.1934E-01 66.33 -1.522 -17.26 -3519. 

0.1314 -145.4 3.577 16.34 0.5214E 05 

0. 41 C5E-01- 15 .21 0.3764 1 . 7t 1 - 1807. 

0.4120E-03-299.5 1.778 -.5376 -161.4 

6090E-G2 16.11 11.88 -.7b60E-01 310.0 

0.2773E-01 29.35 10.44 0.3425 301.6 

2121 E- 01 46.18 -1.270 9.650 -1060. 

0.1767 3.432 -.1280 -.5200 0.2296E 05 

0. 2580L-C2 0.271 IE-01- . 17 59E-02-. 7722E-02 340.4 
24.47 664.1 -17.71 -81.35 -.2642E 06 

0.5307 5.621 -.4019 -1.9C4 0.2597E05 

0. 238 2 E-C 10.3991 -. 1873E-01-. 64445- 01 10b3. 

0.3357 29.72 0.3654 -1.290 9662. 

1354E-02-4.450 -.7336 0.3465 17.73 

0.777 IE- 01-239.0 -83.75 21. 1C 162.5 

THE C MATRIX 

0.51 18E-G1-. 2216 5.798 111.8 -69.26 3.520 0.2494 -.1784 

0.4174E-C20.1718E-05-.6230E-040.2339E-03-. 1401E-01-. 3 1 C5E-04 0 . 3074 E-0 4- . 4 427E-05 
0.0000 O.COOC O.OCOO O.COOO 0.0000 0.0000 0.0000 0.0000 

0. 15o5£-030 .65C9E-05-.6405E-04C. 3292E-03-. 6 3 20 E-0 1-. 1005E-03 0. 10 21 E-0 3-. 4 650E-06 
0. 20 64 F-0 4 0. 1 44 7 E-0 3- . 1057E-01-. 3 652E- 03 0 . 7 18 7E-0 1-. 1 354E-02 1 143 E-0 30. 3800E-06 
0. 1645E-C4-.3C27E-04 0.3587 E-0 30. 31 1 1 E-C3-. 587 3 E-0 10 . 4691 E-03 0. 2237 E-0 3-. 3050E-06 
0. 9703E-0 50. 462CE-05-. 5000 E- 040.763 OE -04- . 15 04E-0 1 C . 3 944E-04 0. 81 1 0E-04-. 1250E-06 

0.1523 0.1516 -.6119 -. 2000E-03 C. 3^00E-C2 1.074 0.6558 0.0000 

0. 2490E-050.2241E-04-.212bE-C5C.OCCO 0 . 297 1 E-0 4-. 3 533 E-05- . 87 15 E-0 50. 0000 

0.0000 0.0000 1.00C O.OOGC 0.0000 0.0000 0.0000 0.0000 

0. 562 5E-06 0.8 900E-06-. 17 87 E-0 50. 5 0CCE-08 0 . 200 0 E-0 5-. 3 1 5C E-06 - . 1 6 92 E-0 4 0. 4 999E-08 
5140E-C6-.5 2 00E-060. 195SZ-05-. 2 00CE- 07- . 1 55 0E-050 . 2 450E -06 0 . 1 8 5 1 E-040 . 0000 
0.2460E-060.1000E— 07-. 1582E-050 . 0 CCO 0. 1 *8 5E-C 5- . 26 50 E-C6- . 1 5 23E-C4 0 . COOO 

0.2950 E-0 7 0.4999 E-G8-. 436QE-0C0. OOCC 0 . 73 0 05-0fc- . 92 50E-07- .3 963 E-05 0.4998E-09 
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THE D MATRIX 

0.0000 66.03 4.466 

0. 9570E-C40.1 07 0 0. 1735 

0.0000 0.0000 0.C000 


-2.460 -62b. 6 

0. 1800E-03-3. 110 
C.C000 0.0000 


0. 1328E-040.2139E-01-.43132-020.6746E-04-.7927 


0. 89 77 E-0 6-. 9 988 E-0 20. 2439 E-030. 1 96 3E-020 . 1 01 7E-01 


0. 1250E-06 0. 93C5E-02 0. 1 1 S 1 E-C 2-. 3588E-03- . 6599 E-0 1 
0. 2500E -0 60. 2435 E- 0 20. 534 1E-0 30. 342 IE- 04- . 16452-01 


EIGENVALUES OF THE A MATRIX 
REAL PART IMAGINARY PART 


-547.34229 
-242.02245 
-63.958786 
-33. 659134 
-33.659134 
-50. 000046 
-43. 1 4595 C 
-38. 926193 
-21. 324432 
-20. 143570 
-20. 143570 

- 13.043441 
-2.0527363 

- 1. 5882740 
-0.84857059 
-0.681 02890 


0.C0C00000 
0.C00000C0 
0.00000000 
16.220993 
-16.220993 
0.00000000 
0.00000000 
o.ooooooco 
0. 00000000 
1 . 0338373 
-1.0338373 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 


MODIFIED EIGENVECTOR MATRIX 


-2. 2230E-03 
7. 5756E-04 
3. 9353E-02 
1.0000E 00 
-1. 1814E-01 
7.G958E-04 
1. 6463E-03 

- 1. 1070E-03 
4. 5213E-04 
5. 9798E-05 

-5.1 750E-03 
-1.9 5 47 £-0 2 
-7. 77 31E-03 
-6. 9 1 6 IE- 04 
2.3238E-05 

- 1. 4625E-02 


-6 . 616 8 E-0 3 
-1.6480E-03 
7.2 256E- 0 1 
1 . 0 00 0 E 00 
-4 . 6 66 4 E-0 1 
5 . 8 06 1 E- 03 
1.4218E-02 
-2. 047 IE- 02 
1 .3700L-G2 
1.7940E-03 
-5.73G8E-02 

2 . 7 48 1 F- 03 
1 .0854E-03 
2.7272E-C2 
-3 . 2092E-C4 
-2.2 158E-02 


-4.9353E-C3 
-4.4930E-C3 
-2. 1614E-01 
-3. 3354E-C1 
-1. 1 4o6E- 01 
7.9861E-03 
1.75C4E-C2 
1.0254E-02 
4. 1554E-01 
5.2113E-C 2 
-1. 2658E-01 
1.0CC0E 00 
3. 8051 E- 01 
2. 1 527 E-0 1 
-6.9C57E-C3 
-3 • 9b 6 1 E-0 1 


-4. 6184E-03 
-2. 3 999E-03 
6. 3475E-02 
5.7825E-02 
-1.7936E-03 
-2.0655E-02 
-6.907 U-02 
9.71 94E-02 
7.29 b4E-0 2 
7. 6854E-03 
3.7397E-02 
6.8 086E-02 
2. 30 2 IE- 02 
-3.C036E-01 
1.00C0E 00 
-2. C736E-02 


2.4202E-04 
7. 2430E-04 
-2.7162E-02 
-2.5127E-C2 
-2 • 5 885E-02 
7.5549E-03 
4. 6350E-02 
-3 .4 953E-G2 
-2. 1 8 1 3 E- 02 
-2.348 9E- 03 
-9.9504E-C3 
-2 . 029 0 E- C 2 
-6. 9710 E- 03 
-3 . 0 055 E- 0 1 
1.0000E 00 
-6.0929E-02 


-3.C669E-C3 
- 1. 64 69 E- 03 
1 .4036E-C2 
1.2364b-02 
-2. 4534E-02 
-1.84422-03 
-6.01 95?- 0 2 
4.65342-02 
5.4464E-C2 
5.6 144E-C3 
3. 1947E-C2 
5. 1 141E-C2 
1.7G21E-02 
-3. 1C422-C1 
1.0000E 00 
-6.5C86E-02 


3. 16882-03 

- 7 . 5 024E-03 
5.7012E-01 
4.5838^-01 
5.68962-01 

-8.3 14 2E-02 

- 1 .9626E-01 
1 .73912-02 
6.70S9E-01 
8. 6 95 BE- 02 

-1. 1496E-01 
1. 88912-01 
3. 8657E-02 
4.8815E-01 
6.2 566 E- 02 
1.000CE OC 

1 . 1 8t 4 E-0 2 

-5. 1074E-03 
1. 54412-01 
1 .4369E-0 1 
3.676 IE-01 
1. 5645E-0 1 
3.4625E-C1 

2.7 09 9 E -0 1 
1 .29552-01 
1 . 1217E-02 
8. 1 3 5 1 E-0 2 
1. 3743E-01 
4.0 1 07 E-0 2 
4.25992-01 

1.GOOOE 00 
7.93062-01 


1.81 442-02 
1.8o52E-0 2 

1 . 7 232E-0 1 
1. 50S0E-01 
3.4104E-01 
2.3245E-02 
6.65952-02 
-3.2996E-02 
-1.0C28E 00 
-1.1 490E-01 
-5.6957E-01 
-1.22522 00 

- 4 . 503 5 E-0 1 
-1. 1 4 1 3 E-0 1 
-1. 73392-01 

1.00C0E OC 

-1. 5433E-01 
-1 .5730E-03 

- 1.41 70 E-02 
- 2 . 54 b4E-02 
-1.2887E-01 
-1. 18102-01 
-1 .19902-01 
-7 . 649 1 E-02 
-3.7265E-02 

1 .0588E-02 
-3. 2300E-02 
-4 . 5209iL-02 
1.00C0E 00 
-8. 76232-02 

- 1 . 3022 E - 01 
-4.8297 E-02 


-1 . 1993E-07 
7. 0025E-08 
-7. 0356E-06 
-1. 1 374E-05 
-7.0541 E-08 
6 . 8296E-G7 
1. 767t>E-08 
2.8418 E-07 
-3 . 5207E-06 
-4.2675E-07 
-3.57 40E- 08 
2. 34932-05 
8. 8073E-06 
1. 1205E-06 
4 • 6238E-03 
1 .COCOE CO 

1.00 COE 00 
2.131 OE-O 2 
7. 1 875 E- 02 
1.45 98 E-0 1 
7.46 1 6 E-0 1 
7 . 32 54 E- 0 1 
7.3025E-C1 
4.71 54E-0 1 
2. 344QE-01 
-1 .16022-01 
2.064 3E- 0 1 
2. 85812-01 
9 . 4 335E-G 1 
3.8514E-01 
7.7193E-01 
1 . 3458E-01 


9 . 5 7 1 0 E- 03 
-1.5317 E-C2 

2. 3 30 3 E-0 2 
-6 . 8452 E-02 
- 1 . 2389 E-02 
-3.659SE-03 
-7.6 025 E-04 
1.7098E-03 
1.000CE CO 
1 . 1974 E-0 1 
1 . 489CE-01 
-6 • 8257 E-0 1 
-2.530 3 E-0 1 
-1 . 4 91 8E-02 
7.4 150 E-04 
-1 .S590E-C1 

1.0CCCE 00 
6.8 185E-01 
9.4588E-01 
8 . 3z42 E-0 1 
7. 33872-01 
6.8965E-01 
6.9093 E-0 1 
8. 5674 E- 01 
7. 1807E-03 
-2.2515 E-0 2 
-6. 4756E-02 
-1 .7694E-02 
-2 . 6224 E- 02 
-9.7 337 E-02 
7 . 2346 E-0 1 
1 . 1804E-01 


4 . 1 2382- C 5 
- 1 . 6 1 42 E-02 
-3.47592-02 
-1 .00472-01 
-7 . 854 0 E-02 
1.0412E-02 
2.6492E-C2 
-7 . 6388 E-03 
1.0000E 00 
1. 1 823 E-0 1 
2 . 4 1 28E-0 1 
9.6690E-C2 
3. 5564 E-02 
-1.5579E-C2 
-2.6372E-C2 
-2.9190E-0 1 

5.4442E-G1 
4.6742 E-0 1 
6.34 15E-0 1 
5.3156E-01 
3. 8498 E-0 1 
3.6S48E-C1 
3. 6771 E- Cl 
5. 2030 E-0 1 
-2.4968E-02 
1.0000 E 00 
-7.07112-02 
-8. 0364 E-02 
- 1 .0557E-01 
-1 .4903E-G1 
3.8460E-01 
3.6657E-02 
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THE MATRIX OF MODE SHAPES 


I is MAG. AND ANGLE (DEG.) 


FORM 



1 

2 

3 

4 

5 

6 

7 

8 

1 

-0. 2197D-02 

-0. 11S3D-01 

-0.4396D-02 

0.7824D-02 

-164.3 

-0 . 407 2D -05 

-0. 9809D-02 

0.5459D-03 

2 

0.6751D-03 

-0.2612D-02 

-C . 30 52 D- 02 

0. 3223 D- 02 

1 18.5 

-0 . 39 0GD-G9 

0.5419D-02 

0. 1856D-02 

3 

0.5541D-02 

1.000 

-C.2642 

0. 3173D-01 

-31. 38 

0.9352D-06 

-C. 1298D-01 

0.2297 

4 

1.C00 

C. 6 272 

-0.3965 

0.3507D-01 

36.84 

0 • 89 60 D-05 

0. 4 7 b7 D-0 1 

0.2474 

5 

-0.91 21D-C1 

-C.7889 

-0.31C4 

C.9526D-01 

-51. 18 

0. 16 C 2D- 04 

C.6525D-C3 

0.3625 

6 

0. 1431D-C3 

0.4772D-02 

C.5277D-02 

3.2216D-C2 

1 15.5 

0.11 16D-06 

0. 1789D-02 

-0 . 6282D-01 

7 

0.3565D-03 

0. 1121D-01 

0.1 225D-0 1 

0 • 57 06 D-Q 2 

132.0 

-0.3389D-06 

0.2 144D-02 

-0.1352 

8 

-0. 2747D-03 

-0.4913D-01 

0.2056D-01 

0. 2794D-02 

-134.7 

-0. 1442D-06 

-0.7412D-03 

0 . 8884 D-0 1 

9 

0.7902D-04 

0. 2243D-01 

0.3832 

0.3229 

-65. 13 

0.6263 D-05 

-0.5015 

-0.7742D-01 

10 

0. 1031D-04 

0.2S96D-C2 

0.4760D-C1 

0. 3*23D-C 1 

-65. 76 

0.7621D-06 

- 0. 6032D-G 1 

-0.828 1D-C-2 

11 

-0. 8355D-03 

-0. 392QD-0 1 

-0.1222 

0 . 4027D-0 1 

22. 10 

0 • 1 879 D-07 

-0.526 ID-01 

-0. 5331D-01 

12 

-0.2523D-C1 

0.7107D-01 

1.00C 

1.000 

0.0000 

0.51 09D-C3 

1.000 

-0.5024D-01 

13 

-0. 1004D-01 

0.2733 D- 01 

0.3309 

0.3754 

-0.4534 

0.1916D-03 

C. 3723 

-0. 1718D-01 

14 

-0. 2626D-C2 

0. 77 3bD- 0 1 

0.2047 

0.1120 

-32.35 

0 . 4552 D-04 

0 • 2530 D-0 1 

0 . 520 ID- 0 1 

15 

- 0 . 26 3 3D- 04 

-0. 1 7 2 5 D-0 2 

-0.50 1b D- 02 

0. 3714D-02 

-36. 69 

0 • 255 3D- 06 

-0.1567D-02 

1.000 

16 

-0. 1674D-01 

- 0. 9920D-0 1 

-0.2274 

0.6603 

-120. 1 

1.000 

0.1936 

0.9312 


9 

10 

11 

12 

13 

14 

15 

16 

1 

0. 2372 D- 01 

0.54 9 ID-0 2 

147. C 

1.000 

0.7954 

1.000 

-0. 1 131D-C2 

-0.1104 

2 

-0 • 54 86 D-02 

0. 2263D-0 2 

157.3 

0.2221 

0.2161 

-0.2517 D-0 1 

0.30blD-C3 

0. 1150D-02 

3 

0. 7662D-02 

0. 1272 

-78.59 

0. 9807 

0.8949 

0.2459 

-0.1 903D-01 

-0.3157D-01 

4 

0 . 89 1 4 D- 0 1 

0. 1251 

-87.25 

0.8595 

0.7537 

0.3479 

-0.37 42 D-0 1 

-0. 4292 D- 01 

5 

0.3857 

0. 1410 

-115.5 

0.9436 

0.7465 

0.9362 

-0.8514D-0 1 

-0. 1084 

6 

0.1096 

0. 36 9 CD-0 1 

97. 55 

0.2400 

0. 1842 

0.2836 

0.4496D-02 

- 0. 3<t08D-0 1 

7 

0.2583 

0.4728D-01 

112.6 

0.2551 

0.1860 

0.3644 

0.4 158D-01 

-0.4 13 2D-0 1 

8 

0. 1628 

C. 1098 

-57. 37 

0.3216 

0. 2904 

0.9558D-01 

- 0. 2 826D-0 1 

-0. 1285D-01 

9 

0. 2086 

0. 1098 

-33.59 

-0.1854 

-0.17C0 

-0.4129D-01 

-0.3490 D-0 1 

0.497 2D-02 

10 

0. 1830D-C1 

0. 1141D-01 

-34.69 

0.4457 

1.000 

0.1945D-C1 

-C.3561D-02 

-0. 1497D-02 

11 

0.1519 

0.6241D-01 

-28.79 

-0.2246 

-0.2061 

-0.51 1 2 D-0 1 

-0. 2 096 D-0 1 

0.6235D-02 

12 

0.2364 

0.9573D-01 

-37. 29 

-0, 2406 

-0.2483 

-0. 14 18D-01 

-0.3975D-01 

0. 154 ID-02 

13 

0. 6974D-01 

0. J209D-01 

-37. 83 

-0.3580 

-0. 3406 

-0.514 ID- 0 1 

-0.1315D-01 

1 .000 

14 

0.7036 

0.8957 

- 162.7 

-0.2767 

-0.2965 

0.8171 D-0 1 

-0.9180 

-0.2540D-01 

15 

0.7780 

1 .00G 

O.OOCC 

0.2679 

0. 1952 

0.3911 

1.000 

-0.4576D-01 

16 

1.000 

0. 3365 

-116.1 

-0.6832D-C1 -0.7976D-01 

0.2124 D-0 1 

-0.1988 

-0. 1 479D-0 1 
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Operating Point 5 


THE A MATRIX 


-2.289 0.1099 5.133 

4 10.9 

-782.2 

-1. 354 

0.8252 

-.5953 

-. 2272 -2. 935 1 28 . b 

-236.1 

-421.0 

14.51 

1 .235 

-1.384 

0.58C1 3.222 -160.4 

-5.380 

1040. 

-45.28 

-3.187 

2.786 

0.3048 -.2448D-01 13C.4 

-563.9 

115.5 

-4.880 

-.3420 

-1.141 

0. 6C28D-C1-.8437D-01 1.121 

21.42 

-111.5 

1. 331 

0.9735D- 

■0 1- . 656QD-0 1 

0. 8305 -. 163 1 D-0 10 . 1253 

5565 

77.11 

-19.77 

-.24 09 

0. 17 00 D-0 2 

0. o223 1403D-010.635CD- 

01-1. 181 

202.0 

0. 1880 

-20.56 

0. 4400D-02 

-.1182 0.7847 30.02 

C.9SC0 

-203.5 

30.24 

0.6468 

-20.01 

-. 9372D-02-.1714D-01C. 2702 

-.231 OD- 

01-16. 18 

C . 7073 

0 . 5 1 98 D- 

01 6.793 

1380D-03-.2527D-030.4012D- 

02-. 3 1 COD- 

03- . 24 03 

0.104SD- 

■01 0.7829D- 

■030. 1007 

-1.376 -6.967 41.26 

12.55 

-24d9. 

107. 5 

7.575 

36.1 1 

314 ID-01 -.4 74 5 -9 9. IS 

372. 8 

-52.34 

2.283 

0.1996 

9.985 

139 7D-C2- .21 09D-01-4 .409 

16.57 

-2.325 

0.1014 

0.8855D- 

■020.4439 

-.1725 -.2573 -45.32 

39.61 

273.6 

1 .654 

0.1053 

4.589 

1266D-C1 0.3C75D-01-. 1616 

8345D- 

Cl 10.88 

-.4758 

19.81 

0. 1286 

-1.698 1.907 -4.850 

216,3 

67C.2 

-29.04 

-2.058 

-.4760 

0.5626 0.5527 -2.238 

5. 0 56 

5 . C94 

-.3931 

-.1265 

-.5000D-05 

6.153 6.151 -2.134 

0. 50 OC D- 

0 40. 1200D 

-0 1-. 1350D- 

02- . 2 1 30 

0.5000D-04 

-. 1500 D-C 1 - . 1 SCO D-0 1 2.975 

C. 1500D- 

02- .1500D 

-010.0000 

0.5320 

0.0000 


-1.543 -1.543 3.123 0.10GCD-03- 

0. 5892 D-01 0.555 OD-Q 1-. 2324 -.65G0D 

2350D-02-.890CD-020.6670D-020. 50CCD 
-.235CD-02-.2500D-020. 1050 D- 0 1 0. 0 CC G 0 
0. 3150D-C20.7500D-C2-. 1C85D-010. OC-C 0 
-50. CO -.4350D-02 39.37 

-.6667 -.6667 0.5833 

0. 3050 D-C1 0.75 CO D-01 -47.8 C 


-9.603 -9.618 

38.40 

-50 .OC 

0 . 500 0 D 

-.4269 -.4277 

1.7 07 

-2. 001 

-2,000 

-4.415 -4.422 

17. o5 

-3.098 

-3.085 

0. 3900 D- 030. 0000 

0. 20 36 D- 

0 1C, 5 1 5 CD - 

03 C . 0000 

0.5910 0.5950 

-2.379 

0 . oocc 

- . 500 0 D 

THE B MATRIX 

-. 1 302E-C 1-15 9.5 

-5b. 50 

1 .322 

739.0 

-. 14 94 E-0 1-3. 97 5 

-2.387 

-32.03 

-4620. 

0.1657 -31.41 

6.132 

29 .28 

0.6 164E 

0.5106E-0 1-.3OC0 

0.6525 

3. 208 

-2079. 

-. 1 82 6L-C2- 6 C 9 . 9 

3.485 

-.7817 

-231.9 

53C0E-0 2-b. C5 0 

10.76 

-.2800 

-93. Q2 

-.6750E-03 151.3 

10.66 

-.4230 

-1062. 

-.6045E-02-.7868 

-1.230 

9.75C 

- b2 8 . 0 

0.6726E-G1 1.911 

-.1095 

-,47b1 

0 . 1252E 

0. 1263E-020. 1210E 

-0 1-. 149 IE- 

02- . 7 055 3- 

02 185.4 

12. C4 5C.d8 

-14.64 

-68.37 

1 5 1 8 E 

3.2753 36.44 

-.3745 

-1.521 

0 . 1296E 

C.1233E-C1 1.747 

1753Z- 

■G1-.6730E- 

01 575.0 

0.1657 25.49 

0 . 2595 

-1 . 058 

5255 . 

0. 4276E-G3- 1.842 

-.7270 

C. 3429 

40.95 

0. 6 1 5 2 E-C 1-57. 27 

-79.75 

18.90 

-659.0 


1 100D-0 10. 4250D-020.5350D-010.0000 
030. COCO 0.3920 0.2389 -.5000D-04 

C30.0C00 0. 1 100 D-C2 0.64 55 D-01 0.0 000 

COCO 0.5500D-030.101 1 -.5000D-04 

0. 5000D-0 30. 50 00 D-04 - . 1069 0.500 0D-04 

-. o5 COD- 03 0 .C000 C . 5 O00D-05 - . 83 15 D-0 20 . 5000D-05 

86C0D-C5C. 5000D-C60.5000D-07-.1333D-030. 5000D-07 
15CCD-02C.2500D-0 1- . 3000D-02 - 1 .2 54 0.5000 D-0 3 

040. 7645 D-02-. 1995 D- 010. 0000 
0.3335D-03-.8645D-030.5000D-06 
-19.76 4061 D- CIO. OCOO 

-. 500CD-05-2 0. 00 O.COOO 
02 30. 8e 20.15 -50.00 


THE C MATRIX 

C. 9222E-01- .4346 5 .794 1 1 1.7 -70.96 6. 841 0.4809 -. 3399 

C. 776 9 E-C 20. 34 71 E-OS- . 5 36 5 E-0 4 0. 2 03CE-03- . 1 o75E-0 1 - . 6 1 34E-04 0 .6 082 E-0 40 . 4 99 9E-08 
0.0000 O.OGOC 0.0000 0.0000 O.COOO G.COCQ 0.0000 0.0000 

0. 1915 E- 030. 5921 E-0 5- .3160 E-0 4 0. 1 57 CfcJ- 03- . 3 07 6E-0 1- . 90 57E -04 0 . 93 3 9 E-04 37S0E-06 
C. 19 1 2E- C4 0.1 34CE-03-. 50 30 E-C 2-. 1 7 dC E- G 3C . 34 20 E-0 1 1 245E-02 - . 1 05 1 E-C 30 . 500 OE-06 
0. 154 9 E-0 4-. 2 962 E-0 4 0. 177 2E-0 30. 1 53 0S-03- .29b9E-C 1 0. 4555E-03 0 .21 80E-0 3-. 4350 E-06 
0. 9 242E-G50.4 tC5E-C5- . 26 06 E-0 40. 376CE-04-. 7652E-0 20. 382 4 E-04 0 .78 77 E-0 4-. 9500E-07 


0.3063 0.3101 -1.21c 0. 5 1 5 CE-030 . 5 55C E-0 2 2.124 1.281 O.OCOO 

0. 0000 O.CCOG -.36tCE-CbG.4999E-080.4994E-C7C.0000 - . 1 4 56 E-C40 . 500 6E-0 8 

0.0000 C.COOO 1.00C O.OOOC 0.0000 0.0000 0.0000 0.0000 

0. 3600E-CfcC.7 OC 0E-06- . 18 ICE-0 5-. 5 OC OS- 07 0 . 5 CCfcE-07- . 3500E-07- . 1 5 8 0 E-04 0 . 5 006 E-0 8 
4550E- Co- .7*99E-06 0. 196 0 E-0 50. 1 C 0 C E-C7 - . 400 OE-C 60 . 3 50G E-07 0 . 1 7 45 E-0 4- . 4 99 9E-08 
0. 3550E-Ce>0. 1 150L-05-. 1 6 75 E-0 5-. 4 C C GE- 07 0 . 4 994E- 07 0 . 85 OOE-07 14 95 E-0 40. 500 6 E-0 8 
0. 1 250E-C60.3350E-C6-. 4350E-CL-. o5C0E- 07 O.COOO 500CE-08- . 38 1 0 E-0 50. 500 6E-08 


45 



THE D MATRIX 

-. 6056E- 02 125.9 8.720 
0. 8989E-C40. 1341 0.3291 
0.00C0 0.0000 o.oooc 


-4.424 -796.5 
0. 235 CE -04 -.1482 
0.0000 0.0000 


-.8716 E-05- .2 044 E-0 1-.5000S-030.941 5£- 04-.1255E-01 
0. 1905E-0S- .3 901 E-0 30. 2003E-030. 22S3E-020. 50C5E-G 1 
5250 E- 060. 2104 E-020. 1 177E-0 2-. 2 756E- 03- . 4 26 5E-01 
0. 7550E-07 0.6 26 0E-030.5243E-030.5582E-04-.7575E-C2 


EIGENVALUES 
REAL PART 

-551. 70679 
-243.72264 
-63.676819 
-32. 630035 
-32. 630035 
-49.993179 
-43. 530563 
-41. 179535 
-21. 002594 
-21.002594 
-20. C30380 
-12.805978 
-3.4393711 
-0.674 19291 
-1. 9637203 
-1.6019726 


OF THE A MATRIX 
IMAGINARY PART 

0.00000000 
0.00000000 
0.00000000 
1 5.138564 
-15.138564 
0.00000000 
O.OCOOOCOO 
0.00000000 
1.1867085 
-1. 1667085 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
o.ooooooco 
0.00000000 


MODIFIED 
-4.2737E 
1.4729E 
3.8654E 
1.0000E 
-1. 1764E 
6. 87 3 IE 
1. 5789E- 
-1. 1055E' 
4. 2739E- 
5.6505E- 
-4.9419E 
-1. 9486R 
-7.7504E 
-9.4880E' 
2.7583E- 
-1. 4 162E 


EIGENVECTOR 
03 -1.2709E- 
03 -3.5627E- 
7.1884E- 
1.0000E 
-4.6543E- 
5.6185E- 
1.3525E- 
-2.0331E- 
1.3125E- 
1.7183E- 


-02 

00 

■01 

•04 

03 
■03 

04 

05 


MATRI 

02 -9 

03 -8 
01 -2 

00 -3 

01 -9 
03 

02 
02 


02 

03 


-03 -5.5351E-02 -1 

03 1 

04 3 
02 2 
04 -6 
02 -4 


02 

•03 

04 


2. 1663E- 
8 . 5 7 4 8 £- 
2.6150E- 
05 -2.7249E- 
02 -2. 1476E- 


X 

.49233-03 

• 9 745E-C3 

• 0549E-G 1 
. 1 676E-C 1 
. 9 7 53E-02 
• 8418S- 03 

• 5553 E-02 
•0999E-02 
. 1 163E-01 

1097E-C2 
.2357E-C1 
. 0 C OOE 00 
.6061 E-0 1 
.1197E-01 

• 25 14E-03 
. 1420E-01 


5* 2 1 64 E-03 
-1.4072E-02 
5.3J82E-01 
4.29872-01 
5. 3424E-0 1 
-8.33742-02 
-2.0 162E-01 
1.93 16E-0 2 
6.4546E-01 
8.3 6 15 E-0 2 
-9.0202E-02 
1 . 6 28 1 E-0 1 
7.8319E-02 
5.C956E-01 
7.9372E-02 
1.0C0CE 00 


3.29122-02 
3.8955E-02 
2. 1683 E- 01 
1.981 6E-Q1 
3.5593E-01 
1 . 4335E-02 
5.4487 E-02 
-4.6820E-02 
-1.C154E 00 
-1. 15852-01 
-5 • 9081E-01 
-1.1816E 00 
-4. 3304E-0 1 
-7.64 02E-02 
-1.8387 E-0 1 
1.0000E CO 


-2 • 8542E-05 
-7 . 82 62E-08 
8. 94 65E-06 
7.8303E-06 
2 . 96 1 1E-05 
8.58 02E -06 
5.3447E-06 
-4 . 94 57 E-0 6 
-3.0531 E-06 
-3 • 68 98E-07 
3. 97 83E-07 
1 . 4 1 0 1 E-03 
5 • 2885E-04 
1.57 68E-04 
-3.85 89E-06 
1 . 00 OOE 00 


2.2520 E-02 
-3.0951 E-02 
3.6179E-02 
-5.7 556 E-02 
1 . 2733 E-02 
-1.2631 E-02 
- 1 . 2204 E-0 2 
1. 3583E-C2 
1.0000E 00 
1.1 983E-01 
1.3951E-01 
-8. 3629 E-0 1 
-3.1 043E-01 
-2.0409 E-02 
1.0061 E-02 
-2.0327E-01 


9.2817E-03 
-3. 1924E-02 
1.1729E-C2 
-6 • 8026 E-02 
-1.2300E-02 
-4 . 0752E-C3 
-8.6401E-04 
9. 7482 E-C 1 3 
1 , 0000 E 00 
1 . 1 904E-G1 
1 . 9082E-01 
-2.4986E-01 
-9.231 5E-02 
4. 1406 E-05 
Q.9559E-04 
- 1 . 3570E-01 


-1. 5797E-04 
2. 4563E- 04 
-4. 1 9 09E-02 
-3.911 3E- 02 
-4.01 35E-02 
1. 6 537E- 02 
1 . 0675E-02 
-3.1 079E-02 
8. 8536E- 04 
-5. 7 229E-05 
5.4 3 19E-0 3 
1.2066E-03 
6. 9758E-05 
-3. 1495E-01 
I.OOCOE 00 
-9. 3137E-02 


-9.3075E-03 -3. 1619E-03 
-5 . 6 954E-03 -1.0C14E-G3 
5. 7564 E- 02 7.5030E-03 

5 . 2452E- 02 6.5271E-03 

2.88C4E-03 -1.2357E-C2 
-1.1206E-02 -6. 9541 E-0 4 
-1 . 1 297 E-0 1 -1.5745E-03 
1.0 133E-0 1 1.5565E-02 

9.1882E-02 1. 5152E-C2 

9.6481 E-03 1 . 5o88E-C3 

4.941 5E-0 2 8 . 56 3 1 E- 0 3 

8 .5069E-02 1.4132E-C2 

2. 8699E-02 tt.711lE-03 

-2 .7 99 1 E-0 1 - 2 . 93 74 E-0 1 
1.0000E 00 1 . 0000 E CO 

4 . 067 1 E-03 -3. 1310E-02 


2.8587E-02 
-1.2351E-02 
1 .5012E-01 
1. 4329E-01 
3 . 96 59 E-0 1 
1 • 788 1E-0 1 
3. 605SE-C1 
2. 8664E-01 
1 . 5466 E-0 1 
1.325 IE-02 
1 .0107S-01 
1.6423E-C1 
4.7o3CZ-C2 
4.8449E-0 1 
1 . 0 0 0 0 E 00 
8. 3959E-01 


1.0C00E 00 
-3. 7C58E-02 
-6.691 4E-03 
8.5612E-02 
8. 31 58E-01 
6.191 SE-01 
8. 1 24 4 E-0 1 
5.4075E-01 
3. 474 6 E-0 1 
-2.6148E-02 
3. 0534E-01 
4. 1 C 24 E-0 1 
-1.7626E-01 
6.23C9E-01 
9 • 5562E-0 1 
3 . 5954E-01 


1 .2778E-01 
1.5826E-01 
2.0351 E-0 1 
1 . 562 1 E-01 
7.781 3E-02 
8 . 35b 0E-02 
7 . 7666E-02 
1 .4993E-01 
- 1 . 32 02E-0 2 
1 . OOOCE 00 
-2.6681 £-02 
-5 . 4770E-O2 
-7. 1373E-02 
-8 . 303 0E-02 
8. 2782E-02 
-7.7338E-C3 


8 . 0 445 £-02 
-5 . 7909E-C3 
-9 • 5454 E-03 
-9.2 146E-04 
5. 0398E-C2 
5.8 383 E-02 
5.4647E-02 
3.3150 E-02 
2. 4859 E-02 
-7. 7684E-03 
2 • 3090 E-02 
3.001 9E-02 
1.0000E CO 
2. 1156 £-02 
5 . 8 603E-G2 
-3. 202d£-03 


1.00G0E CO 
6.9205E-C1 
9.4362E-01 
8. 3362E-C 1 
7.4002E-C1 
7 . 1359 E-C 1 
6 . 8973E-01 
9. 057 9 E- Cl 
5.94 1 0E-02 
-2.91 34E-02 
-1.3570E-02 
4.1565E-02 
1.4153E-01 
-3.5 185E-02 
7.5 207 £-01 
1.8230E-01 
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THE MATRIX OF MODE SHAPES, IN MAG. AND A NGEE ( DEG. ) FORM 



i 

2 

3 

4 

5 

6 

7 

8 

1 

-0 . 4274 D- 02 

-0. 1271D-01 

-0.9492D-02 

0.2356D-01 

-35.99 

-0. 2854D-04 

0 . 2252 D- 0 1 

0.9282D-02 

2 

0. 1473D-02 

-0. 3563D-02 

-C . 8974D-02 

0. 2929 D- 01 

-64.86 

-0 . 76 26D-07 

-0.3095 D-0 1 

-0.3 1 92D-0 1 

3 

0 . 38 65 D-0 1 

0.7188 

-0.2055 

0.4074 

22. 89 

0.8946D-05 

0.3618 D- 01 

0.1 173D-01 

4 

1.000 

1.000 

-0.318B 

0.3347 

20.25 

0.7830 D— 0 5 

-0. 5756 D-0 1 

-0.6803D-01 

5 

-0.1176 

- 0 . 4b 54 

-C.9975D-01 

0.4539 

11. 33 

0.296 ID- 04 

0. 1273D-01 

-0. 1 23 OD -01 

6 

0. 6873D-03 

0. 561SD-02 

C. 78420-02 

0.6400D-01 

-125.9 

0 • 858CD-05 

-0. 126 3 D-0 1 

-0.4075D-02 

7 

0. 1579D-02 

C. 1352D-01 

C . 1 5 55 D-0 1 

0. 1477 

-119.9 

0 . 5345D- 05 

-0. 1 220D-01 

-0. 8640D-03 

8 

-0. 1 1 05 D- 02 

-0. 203 3D-0 1 

0.110CD-G1 

0. 3581 D- Cl 

112.6 

-0.4946D-C5 

0. 1 356D-01 

0.9748D-02 

9 

0.4274D- 03 

0. 1312D-01 

0.4116 

0.O508 

102.6 

-0.30 53 D-05 

1.000 

1.000 

10 

0. 565 0D- 04 

0. 1718D-C2 

0.51 1 G D-0 1 

0. 10 10 

99. 18 

-0.36 C 0D- 06 

0.1198 

0.1190 

1 1 

-0.49 4 2D- 02 

-0. 5535 D- 01 

-C.1236 

0.4226 

143.7 

0. 39780-06 

0. 1395 

0.1908 

12 

-0. 1949D-01 

0.2166D-C2 

1.0C0 

0.8434 

127.2 

0. 14 10D-02 

-0.8363 

-0.2499 

13 

-0. 7750 D- 02 

0. 36 75D-03 

C.3906 

0. 31 12 

124.7 

0 . 5288D- 0 3 

-0.3104 

-0.9231D-01 

14 

- 0 . 9488D- 03 

0. 26 15D-01 

C.212C 

0.3643 

53.53 

0. 1577D-03 

-0. 2041 D-0 1 

0.4141D-04 

15 

0.2758D-04 

-0.2725D-03 

-0.6251D-C2 

0. 1416 

111.7 

-0. 3659D- C5 

0. 1006D-01 

0.8956D-03 

16 

-0, 1 4 1 6 D-0 1 

-0. 2148 D- 01 

-0.4142 

1.000 

o.ocoo 

1.000 

-0.2033 

-0.1357 


9 

10 

1 1 

12 

13 

14 

15 

16 

1 

0.6582D-02 

136.0 

-C.316cD-02 

3 . 2d 59 D-C 1 

1.0C0 

0. 1278 

0.8045D-01 

1.000 

2 

0. 4 172 D- 02 

1 32. C 

-0.10C1D-02 

-0. 1 2350- C 1 

-0.3706D-01 

0. 1563 

-0. 5791 D-0 2 

0.6920 

3 

0. 5C 35D-01 

-81 . C 6 

C . 7 5 C 3 D-0 2 

0. 15C 1 

-C. 6691D-C2 

0.2C35 

-0.9545D-02 

0.9436 

4 

0 . 4 5 27 D- 0 1 

-8 1.71 

0. 6527 D- 02 

0.1433 

0.856 ID-01 

0. 1562 

-0.9215D-03 

0.8336 

c 

0.2845D-01 

-1 3C. 5 

-0.1^3^0-01 

0 . 3966 

0. 6.310 

0.778 1D-C 1 

0 . 5 04 0D- 0 1 

0.7400 

6 

0. 14 13 D— 0 1 

7 9.12 

-C.0954D -03 

0. 17db 

0.8192 

0 . 83 56 D-0 1 

0.5838D-01 

0.7136 

7 

0. > } 024D-C1 

12S .o 

-C.157f D-02 

0 . 3 b r 6 

0.6124 

0.7767D-01 

0.5465D-01 

0.6897 

8 

0. 7495 D- 01 

-o2 .05 

C . 1 5 j" D- 1 1 

9.2866 

0. 5408 

0. 1499 

0.3315D-G1 

0.9058 

G 

0. b4 97D-G 1 

4 . 4 5 

0. 1515D-01 

1547 

0.3475 

-0. 13 20 D-0 1 

0. 24 86D-0 1 

0. 5941 D-0 1 

1C 

0.cd22D-C2 

-4b. 34 

U . 1 5 c 9 L- 02 

0. 1 J 25 D- n 1 

-0. 2615D-C1 

1.000 

-0. 7768D-Q2 

-0.2913D-01 

1 1 

0 . 1 5 D-C 1 

- 3 o . 7 3 

C . 3j>0 3D“02 

0.1011 

0.3053 

-0 . 2668D-01 

C.2309D-01 

-0.1 357D- 0 1 

12 

O.oCltD-CI 

-44.15 

C . 14 1 3D- 3 1 

0 . 1 o 4 2 

0.4102 

-0. 5 4 77 D- 01 

C. 3002D-01 

0.4156D-01 

n 

0. 20 29 D- 01 

-44 . c6 

C .47 1 1D-C2 

C .47530-0 1 

-C. 17 6 3 

-0.7137D-C1 

1.000 

0.1415 

14 

C. 2^-79 

-17:. 6 

-C.2937 

0.4645 

0.6231 

-0.8303D-01 

0.2116 D-0 1 

-0.3518D-01 

1b 

1 .COC 

C.0CGC 

i.?:c 

I.C^O 

0.8556 

0.8278D-C1 

C.586GD-01 

0.7521 

ie 

C.o592D-G 1 

-1 32.- 

- n . 3 1 3 Vj- 9 1 

0.^3 97 

0 . 35 Q 5 

-0.7734D-02 

-0.3203D-02 

0.1823 
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Operating Point 6 


THE k MATRIX 


-1.015 0.4626 5.427 

430.7 

-660.5 

-5.168 

-1.552 

-.2493 

-.2271 -1.852 102.8 

-227.3 

-70.96 

11.79 

-.7220 

-.6399 

0.5935 1.917 -115.9 

0.3312 

185.9 

-20.86 

1.888 

0.9431 

0. 4220 D- 010. 2219D-01 97.88 

-500.4 

175.0 

-1.445 

0.1167 

-.2981 

0.6233D-01-.4b57D-010.7580 

19.70 

-89.76 

0.5613 

-.1039 

-. 2568D-0 1 

0.3507 -.1233 C.2842 

0.2649 

53.81 

-18.60 

0.5443 

5000D-05 

0. 277d -.3571 0.7406 

0.9950D- 

01 130.8 

4.069 

-18.67 

0.0000 

-.4214 0.5473 83.36 

1325 

-132. 1 

28.20 

-1.336 

-20.00 

-.5023D-020.6630D-02-.7106 

C. 2485 

-1.766 

0.1806 

- * 1 7 94 D- 

■01 5.732 

7502D-04 0. 986 ID-04- . 1062D- 

*0 10.4418D- 

02- • 2569D- 

-0 10. 26 7 3D- 

02-.2690D- 

•C 30. 849 ID-01 

-2.925 -9.041 27.14 

-.7450 

-917.3 

102.8 

-9.271 

42. 16 

0.58C7D-C1-.1009 -227.6 

800.7 

17. 66 

-2.066 

0.1794 

8.816 

0.2580D-02-.4483D-02-10. 12 

35.59 

0.7834 

9180D- 

010.7982D- 

■020.3918 

5120 D-01-. 1301 -103.2 

106.6 

479.4 

0.2203D- 

01-.5740D- 

01 4.001 

-.4983 D-0 20. 1 56 1 D-0 1 0 . 775CD- 

■0 10. 2982 

2.119 

-.1624 

19.97 

0. 5859D-0 1 

-1.360 1.050 8.851 

390. 9 

85.89 

-12. IS 

2. 318 

-.3778 


0.3027 0.3026 -1.126 2.317 2.302 -.1776 -.1923 1000D-04 

3.039 3.043 -1.103 -. 1 0 COD- 05 0 . 2 1 7 1D-0 1- . 50 05D-02 - .67 07D-0 1 0 . 0C 00 

O.OOCO - . 1 0G0D-04 1.186 0.00C0 0 . 2170D-C 10. 9999D-020. 1789 -.5000D-05 

-.7345 -.7165 1.490 0. 2450D- 02- . 2150D-01 C. 5755D-01 0 .2240D-0 10. 5000D-04 

0. 288 5D- 01 0.2 9 05D-0 1-. 106 2 G. 1 c OC-D-C 30 . 550 0D-020 . 1 557 0.19 09 -. 5000D-05 

0.0000 0.0000 -.5000D-05-.50COD-C50.0000 0. 5000D-02 0.52 76 D-010. 0000 

0.0000 -. 5 OOOD-C4 0.0000 -. 5 GOOD- 05- . 5000D-040. 2495D-020. 1270 -.5000D-05 

0. 5000D-050.5000D-040. 3100D-030. 5000D- 050 . 5000D-C4-. 6 995D-02- . 126 1 0.5000D-05 

-50.00 0.1 593D-0 1 39.89 0.0000 0.0000 0. 1250D-02- . 1 500D-05-. 5000D-06 

-.6666 -.6663 0.5910 0. 1 0C0D- 07- .5 000D-070 . 1 66 6D-04 0 . 0000 -. 1000 D-Q 7 

-.3550D-020.5000D-03-48.21 0. 1 250D -020 . 0000 -. 4370D-0 1 - . 88 1 0 C.5000D-04 

-10.71 -10.72 39.63 -50.00 -.1629 0.3749D-02 0. 1788 D-0 10. 5000D-05 

-.4761 -.4766 1.7b1 -2.000 -2.008 0. 1 500D-03 0 . 8 047D-0 30. 5000D-06 

-4.864 -4.872 18.01 -2.534 -2.569 -19.84 0 .3934 D-0 1 0 . 0000 

0.000 0 0 . 2 867D-0 1 0 . 1043D-0 1 0. 0 000 -.3910D-0 1-. 1 500D-02- 19. 99 0. 00 00 

0.5726 0.5175 -2.118 485 0D- 02 0 .0000 22.31 28.64 -50.00 


THE B MATRIX 


0. 20 38 E-0 1-9.805 

-9. 222 

5.907 

-1085. 

0. 374 6E- 01-3.854 

-2.686 

-10.50 

-2153. 

0.9210E-0 1 23.46 

6.906 

23.73 

0.4C57E 05 

0.1 774E- 0 1 1.547 

0.7326 

1.627 

-1690. 

-. 78 6 5 E-03- 16 9.9 

0.7584 

-.6619 

-126.4 

0. 1 730E-02 5.843 

-.3999 

-1 . 594 

806.2 

-. 26 8 3 E-0 1 12.04 

3.075 

-4.668 

-75.32 

8026 E- 01- 24. 82 

-4.828 

-14.57 

-2564. 

0.1188 - . 3885E- 

•0 1 - . 56C9E- 

-01-. 1 965 

0.1783E 05 


0. 1694E-020.2137E-02-.7912E-03-. 2749E-02 264.7 
38.53 -117.2 -34.03 -117.6 -.2212E 06 

0.3485 2.811 0.6698 2.336 C.2261E 05 

0. 1594 E-0 10. 1293 0 . 2956E- 0 1 C. 1 034 1C05. 

0.1864 -5.033 • 0.7565 -.2186E-G1 9765. 

853 0E-02-. 7192 - . 706 1 E-0 1 C. 1 482 -117.2 

0.4443E-C1 20.90 -16.89 13.77 497.9 

THE C MATRIX 

-.6229E-C1-.2344 3.136 79.39 185.9 2.701 -.2590 -.1055 

0.11 32E-G10. 1358E-02-. 2907 E-0 2 -. 6250 E-03 - . 6 156 1 554E-01 - .63 17E-G20 . 4850E-05 

0.0000 0.0000 0.000C 0.0000 0.0000 0.0000 0.0000 0.0000 

D. 156 9E- 0 30. 925 IE-04- . 1 9 35E-03-. 4 4 85E-04- . a 1 32E-0 1-. 1 0 56E-02- .4 1 9 9 E-0 30. 1550E-06 
0. 7778E-G40. 2 000 E-0 3- . 1460E-0 1C. 2 305E- 04 0 . 2438E-C1 - . 1 7 07E-02 0 . 24 66 E-0 3-. 1075E-06 
0. 9565E-040. 9280E-05-. 165GE-04-. 7dOCE-Ci<- . 6 786 E-0 1 1047E-03- . 12 OOE-030. 1999E-07 
0. 3232E-G4 0.4 083E-05-.7 250E-05-. 1 37 0E- 04- . 1 5 80E-O 10.1 37 CE-04 0 . 1 7 98 E-04 0 . 130 0E-07 

0.1170 0.1089 -.4287 1 07 OE- 02- . 1 720E-0 1 0 . b22 6 0.7798 -.4999E-05 

-. 33 5 OE- 05 0.0000 O.OOOC G.3060E-05-. 1350E-C4-. 3 28 CE-04- .6007E-030. 0000 

O.OOCO 0.0000 1.000 o.oeoo O.OOCO 0.0000 0.0000 0.0000 

-. 1030E-06-.3450E-06C.265CE-070. OCOO 0. 1000E-C7-. 1 969E-05- . 39 9 9 E-04 0. 0000 
0. 4300E-0 70.38 50E-06-. 1 150E-07-. 4 8C0E- 07- . 5006E-080 . 1 1 2 9E-05 0 .2353 E-04 -. 5002E-09 
305DE-060.5000E-07C.OOOC 0. 4 999E-080 . 0000 3350 E-05- . 65 39E-040. 0000 

-. 1 120E-0 6-.3500E-070.4 998E-090. 5 0C2E- GSO . 4 995E-C6- . 81 75E-06- . 1542 E-04 0. 1000E-08 
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THE D MATRIX 

-. 60 15E-03-48.41 -.7564 -3. 077 -301.0 

0. 2099E-03-.5145E-010. 1373 0.1794E-G1 2.264 

0.0000 0.0000 C.0000 C.OCOO 0.0000 

0. 433 5E-C5- .6 3 75E-0 2- . 4092E-0 2C. 1249E-020.5119 
0. 450 IE-050. 32 6 9E-0 20. 9 029E-0 3G. 4 112E-02-. 1435 
0. 19 35E-05- .6512E-020.9289E-040. 1 269E-0 3- . 9 6 26E-0 1 
0.7040E-06-. 1519 £-020. 307 IE-030. 5360E-04- . 125 32-01 


EIGENVALUES OF IHE A MATRIX 


REAL PART 

IMAGINARY PART 

-506.73145 

0.00000000 

-139.98358 

0.00000000 

-64. 450974 

0.00000000 

-49.715607 

5.7025499 

-49.715607 

-5.7025499 

-49. 998520 

0.00000000 

-45. 165405 

0.00000000 

-22. 648819 

0.00000000 

-13. 368243 

0.00000000 

-20. 343048 

0.84707332 

-20. 343048 

-0.84707332 

-0.89456981 

0.00000000 

-0.76271462 

0.00000000 

-1. 6238403 

0.00000000 

-19. 166748 

0.00000000 

-1. 9968775 

0 .cocooooo 


MODIFIED EIGENVECTOR MATRIX 


- 2. 1 9 74 E- 03 
6 • 75 07E-04 
5. 54 08E-03 
1.0000E CO 

-9. 1210E-02 

1. 4310E-04 
3.5647E-04 
-2. 7467E-04 
7. 9024E-05 
1.C306E-05 
-8. 3550E- 04 
-2. 5229E-02 

- 1 . 0036E-02 
-2. 6257E-03 
-2. 6327E-05 

- 1. 6740E-02 


-1 .1535E-02 
-2.8 1172-03 
1.0000E 00 
8.272 IE-01 
-7.3886E-01 
4. 7721 £-03 
1.121 IE-02 
-4.910 IE-02 
2 • 2425E-02 
2.8955E-03 
-3.9198E-02 
7. 1069E-02 
2.7833E-02 
7.7346E-02 
-1.7253E-03 
-9.9198E-02 


-4 .3955R-03 
-3.0520E-C3 
-2. 6421E-01 
-3. 98482-01 
-3. 1C38E-C1 
5.27 67 E-C 3 
1.2250E-C2 
2.05592-02 
3.8325E-C1 

4 . 7 c 02E- 02 
-1. 22252-01 
1.0C00E CO 
3. 8C93E-01 

2 . 0 468E-0 1 
-5.0 1572-03 
-2.2742E-01 


-4.071 5E-06 
-3.8995E-10 
9. 3516E-07 
8. 9599E-C6 
1.60172-05 
1.11 62E-07 
-3. 3887E-07 
-1.4419E-07 
6.26792-06 
7. 62 08 E-07 
1. 87862-08 
5. 1092E-04 
1.91642-04 
4.55 17E-0 5 
2 . 5534 E-07 
1.0000E 00 


-9.8094E-03 
5.4193E-03 
- 1 . 298 3 E-0 2 
4 . 7567 E-0 2 
6.5245E-C4 
1.7890E-C3 

2 . 1 43 8 E-C 3 
-7.4125E-04 
-5.0 153E-01 
-6 .03212-02 
-5.26102-02 
1.0000E 00 

3.7 233 E-0 1 
2 . 5302 E- 02 
-1.56722-03 
1. 9361 E- 01 


2. 3719E-02 
-5.4857E-03 
7. 6615E-03 
8.914 4 E-0 2 
3. 8574E-0 1 
1.0957E-01 

2. 5 83 3 E-C 1 
1. 62762-01 
2. 08602-01 
1.830 IE-02 
1. 51 932-01 
2. 3643E-01 
6 . 9742E-02 
7 . 0 3 64 E-0 1 

7.7 7 98E-0 1 
1.0000E 00 


-1 .6139 E-03 
-1. 2153 E- 03 
-9.9503E-02 
-1.1695E-01 
-1. 8798 E- 01 

3 . 1 734E-02 
2 . 546 7 E~0 2 
-3.3268E-02 
2 . 9 65 5 E- 02 

2 . 8 89 IE- 0 3 
2. 46342-02 
1. 8153 E- 02 
5 . 6 624 E-03 
-1. 1214E 00 
1.0000E 00 
-4.52982-01 


-7. 5963E-C3 
-2.9 599E- 03 
1 . 4 583 E-C 1 
1. 30862-C1 
6.6671F-C2 
-4.1431E-C2 
-6. 1825E-02 
1.5 17 1 E-0 1 
1. 5240E-01 
1.58762-02 
6.4757 E-02 
1. 3416S-01 
4 . 50 23 E-0 2 
-5.8915E-G1 
1.0C002 00 
1 .5437E-01 


-9.6448E-03 
1.2947E-03 
1 • 057 1 E-02 
4.9058E-02 
- 1 .45012-02 
1.04372-03 
4.2017E-04 
-3. 95 06 E-03 
-1.5715 E-0 1 
-1 . 9 b 96E- 02 
5.2462E-02 
1.000CE 00 
3.7240E-01 
3.4682E-02 
7.5954E-04 
-5. 024 IE-01 

1 . 000CE CC 
2.2207E-C1 
9. 8 07 OE-O 1 
6 . 5 94 6 E-0 1 
9 . 4 356 E- 0 1 
2.4CC2E-01 
2.5506E-01 
3.2162E-01 
-1.83412-01 
4 .4566E-01 
-2 . 2 4 6 1 E-0 1 
-2.4C62E-01 
-3.5 80 OE-O 1 
-2. 7674 E- 01 
2.6769E-01 
-6.Q323E-02 


-5 . 4 22 2 E-03 
-4.37 033-03 
4.3616 E— 02 
7. 0406 E-03 
1 .3393E-01 
-2. 9543E-03 
-6.0587E-03 
2 . 0596E-05 
4 . 2881 E-0 1 
5.1947E-02 
2. 2163E-02 
1.0000E 00 
3.7834E-01 
1. 5458E-01 
5.1 97 0E-03 
2. 403 1 E-C 1 

7. 9 540 E-0 1 
2. 1612E-01 
8. 9492E-01 
7 . 53 66 E-0 1 
7.4654E-01 
1 . 842 4 E— 0 1 
1 . 85 95E-0 1 
2. 9036 E- 01 
-1.7 002 E -0 1 
1.0000E 00 
-2. Obi 3E-0 1 
-2.4827E-01 
-3. 4064E-01 
-2.3649E-01 
1.9523E-01 
-7 . 9764 E-02 


1 . 00 COE 00 
-2.5169 E-0 2 
2 . 4587 E-0 1 
3.47 94 E-0 1 
9.36 15E-01 
2 . 83 56 E— 0 1 
3 . 6444 E-0 1 
9.5577E-02 
-4. 1 2 90 E-0 2 
1 . 94 52E-02 
-5. 1124E-02 
-1 .4180E-02 
-5.14 1 IE-02 
8. 1712E-02 
3. 9 108 E-0 1 
2. 1239E-02 


-1.1 306E.-03 
3.0615E-04 
-1 .9025 E-02 
-3.7 422 E-02 
-8. 5138E-02 
4.495OE-03 
4. 1578 E-02 
-2.8259E-02 
-3.4904 E-02 
-3.5608 E-03 
-2.0 958 E-02 
-3.97 52 E-02 
- 1 . 3 1 49E-02 
-9.1803 E-0 1 
1.0000E 00 
-1 .9882E-01 


5.4593E-04 
1 .8559E-03 
2. 2971 E- 01 
2.4739E-01 
3. 6252 E-0 1 
-6.2816E-02 
-1.3524 E-0 1 
8 . 8 84 3 E-0 2 
-7.74 1 7E-02 
-8.2808 E-03 
-5. 3315E-C2 
-5.0237 E-02 
-1.7179 E-02 
5.2013E-C2 
1.0000E 00 
9.31 19E-C 1 

-1.1 038E-C 1 
1.1496 E-03 
-3. 1574E-02 
-4. 2919 E-02 
- 1 . 084 1 E-0 1 

- 3. 20 85 E-02 
-4. 1316 E-02 
-1 .2849E-02 

4. 5725 E-03 
-1.4966 E-03 
6.2350E-03 
1 . 54 1 0 E-03 
1.G000E 00 
-2.5 397 S-02 
-4.5756E-02 

- 1 . 4766 E-02 
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THE MATRIX OF MODE SHAPES, IN MAG. AND .ANGLE { DtG . ) FORM 



1 

2 

3 

4 

5 

6 

7 

8 

1 

-0. 22 23D-02 

-0. 6617D-02 

-0.4936D-02 

G. 1302D-C1 

-35.09 

-0. 1199D-06 

0. 9571 D-02 

0.4124D-04 

2 

0.7576D-03 

-C. 1848D-02 

-0.4493D-02 

0. 1422D-C1 

-66.91 

C. 7 0 02D-07 

-0. 1532D-01 

-0. 1614D-01 

3 

0. 3935D-01 

0.7226 

-0.2161 

0.4211 

26. 18 

-0.7036D-05 

0.2330D-C1 

-0.3476D-01 

4 

1.000 

1.000 

-0.3335 

0.3412 

26.78 

-0.11 37D-04 

-0. 684 5 D-0 1 

-0.1005 

5 

-0.1181 

-0.4666 

-0.1149 

0.4691 

14.06 

-0.7054D-07 

-C.1239D-01 

-0.7854D-C1 

0 

0.7C 96D-Q3 

C.5806D-02 

0.7966D-02 

0.6104D-01 

-119.4 

0 . 6830D-06 

-0. 3659 D-0 2 

0. 1 04 1 D-0 1 

7 

0. 1646D-02 

0. 1422D-C1 

C. 1750 D-01 

C. 1466 

-116.3 

0.17 6 8D-07 

-C.76C3D-0 3 

0. 2649D-0 1 

3 

-0. 1107D-02 

-0. 2047 D- 01 

0.1025D-01 

0. 26 37D-01 

107,2 

0.2842D-06 

0.171 0D- 02 

-0.7839D-02 

9 

0.4521D-03 

0. 1370D-01 

0.4195 

0.8531 

101.2 

-0 . 3521 D-05 

1.000 

1.000 

10 

0. 5S30D-04 

0. 1794D-02 

C. 521 ID-01 

0.1019 

97.88 

-0.4267D-C6 

0.1197 

0.1182 

11 

-0. 5175D-02 

-G.5731D-01 

-0.1266 

0.4247 

146.0 

-0 . 3574D-07 

0. 1489 

0.2413 

12 

-0. 1S55D-01 

0. 274 8D-02 

1.0CC 

0.6766 

126.2 

0 • 2349D- 04 

-0.6826 

0.966 9D-0 1 

13 

-0.7773 D- 02 

0. 1085D-C2 

C . 38 05 

0.3246 

123.9 

0.8807D-05 

- 0. 2530 

0. 3556D-0 1 

14 

-0.6916D-03 

0. 2727D-01 

0.2153 

0.3545 

58. 16 

0.1120D-05 

-0. 1492D-01 

-0. 1 556D-0 1 

15 

0. 2324E-04 

- 0. 3209D-03 

-0.69 06 D-0 2 

0. 1303 

115.2 

0.4624D-07 

0.7415D-03 

-0 . 2637D-0 1 

16 

-0. 1462D-01 

-0. 2220D- 0 1 

-0.3566 

1.000 

0.0000 

1.000 

-0.1959 

-0.291 9 


9 10 11 12 13 14 15 16 


1 

-0.4818 D-02 

0.2175D-02 

130.5 

0. 1186 D-C 1 

-C. 1543 

1.000 

1.00C 

0.5444 

2 

-0. 2400D-02 

0. 1272D-02 

111.3 

-0.51 07 D-02 

-0.1573D-02 

0.2131 D-0 1 

0.6818 

0.4674 

3 

0.6348 D-01 

0.2162D-01 

-107.7 

0. 1544 

-0. 1417D-01 

0.7 188D- C 1 

0.9459 

0.6342 

4 

0. 5782D-01 

0. 1980D-01 

-108.8 

0. 1437 

- 0 . 2548 D-0 1 

0. 1460 

0.8324 

0.5316 

5 

-C. 1794D-C2 

0. 2522D-C1 

-178.6 

0.3676 

-0.1289 

0.7462 

0.7339 

0. 3850 

6 

-0.2065D-01 

0. 5459D-G2 

58.72 

0.1565 

-0.1161 

0.7325 

0.6897 

0.3685 

7 

-0.6907D-01 

0.5372D-01 

97. 4C 

0.3482 

-C.1195 

0.7302 

0.6909 

0.3677 

8 

0. 971SD-01 

C. 4 1 1 5D-0 1 

-31. 51 

0. 2710 

-0 . 7649 D- 0 1 

0.47 15 

0.8567 

0.5203 

9 

0 ♦ 7 296 D-0 1 

0. 4149D-01 

-66.83 

0.1295 

-0. 37 26 D-01 

0.2345 

0.7181 D-02 

-0. 2497D-0 1 

1C 

0.7695 D-02 

C. 4303C-C2 

-67. 7C 

0. 1 122D-C1 

0 • 1 059 D -0 1 

-0.1160 

-C.2252D-01 

1.00C 

1 1 

0. 3740D-01 

0.2366D-01 

-62.30 

0.81 3 5D-C 1 

-0. 3230D-01 

0.2064 

-0.6476D-01 

-0.7 07 ID-01 

12 

0 . 6 8 09 D- 0 1 

0.3890D-01 

-66.64 

0. 1374 

-0 . 452 1 D-G 1 

0.2858 

-0. 1769D-0 1 

-0.803 6D- 0 1 

13 

0.2302D-01 

0. 1 301D-01 

- 67 . 27 

0.4011 D-0 1 

1.000 

0.9434 

-0.2622D-C1 

-0. 1056 

14 

-0.3004 

0. 3055 

175.1 

0.4260 

-0.87 82D-0 1 

0.3851 

-C. 97 34 D-01 

-0.149C 

15 

1.000 

1.0CC 

C.OOCO 

1,000 

-0.1302 

0.7719 

0.7235 

0.3848 

16 

-0.20 74 D-01 

C . 6 3 04 D-0 1 

178. 1 

0.79 31 

-0. 4 8 30 D-01 

0.1346 

0.1180 

0.3666D-01 
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APPENDIX E 


COMPARISON OF NONLINEAR AND FULL-STATE LINEAR 
MODEL TRANSIENT RESPONSES 

Presented in this appendix is a comparison of transients run on the nonlinear and 
full-state linear models. Responses of states and outputs are shown in figures 9 and 10 
for a + 3-percent change in fuel flow and nozzle area. The operating point is sea-level- 
static idle. (The plots already shown in the main text are omitted.) The responses are 

for P T7 , P T4>5 » P T 7M’ T T2 . 5H* T T2.5C’ T T3’ T T5’ T T6C’ T T7M’ WFAN ’ T T4’ 
SMAF, SMHC, DP 2 5Qp> an d DP 2 55c* The offset at zero time for the nonlinear 

transient is due to the small startup transient in the simulation. Sometimes these 
startup-transient offsets seem quite large. For example, for the SMAF comparisons 
for a +3-percent change in fuel flow (fig. 9(f)), the startup transient caused an offset of 
-0.00075 at zero time. The actual surge margin is 0.10418, thus the drift is 0.7 per- 
cent of the actual. This is a very small error, and the startup transient is essentially 
negligible. The scale used makes the difference look worse than it actually is. The 
immediate response at zero time for the output responses (e.g., fig. 9(f)) are inherent 
in the linear model (ref. 8). Output variables can change instantaneously with time. In 
the linear model this is accomplished through the D matrix. The steady- state errors 
in the comparisons are due to the averaging method used to generate the B and D ma- 
trices. If these matrices had been calculated by forcing the steady-state match only in 
the positive direction, the steady-state error would be almost zero, and only off by the 
startup transient. Note from the comparisons that the linear model at sea-level -static 
idle does approximate the dynamics of the nonlinear simulation quite well. The same 
comparisons were made at the other operating points, and the results were as good or 
better than those presented. Thus, the linear models can be used to represent the en- 
gine over its full operating envelope . 

Some states of the original 16, specifically, T T4 5HI , T T4 5LO , T T4HI , and 
T T4LO* are not shown because they were not output in the nonlinear simulation. 
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APPENDIX F 


THE PROBIFM OF STARTUP TRANSIENTS 

Startup transients occur when the calculations made at zero time in a simulation, 
in which no input is provided, have a small error. For an engine simulation this means 
that a balance point is not reached (i.e. , there is an error between the compressor and 
turbine match) . In steady-state calculations this is not a problem because a balance is 
forced. However, in transient calculations a different calculation path is generally 
made through the simulation, and, depending on how tables are read and how quantities 
are calculated, small errors may result. 

In general, most digital transient simulations do have small startup transients. 
Usually, however, when gross transients (e.g., acceleration from idle to full power) 
are being studied, startup transients are negligible. 

In model validation small transients are of interest, such as a +3-percent change in 
nozzle area. Because a small nozzle-area change has a small effect on an engine, a 
question arose as to how a startup transient would effect this validation. In the F100 
simulation the startup transient was essentially zero so no problems resulted. To in- 
vestigate this problem, a small error was introduced at zero time. This was accom- 
plished by changing the pressure ratio across the compressor slightly from the steady- 
state calculation. This caused a difference between the power extracted by the com- 
pressor and that supplied by the turbine. The results are shown in figure 11 for a 
5-percent change in nozzle area input to the simulation at 0. 1 second. From zero time 
to time equal o.l second, nothing was done to the simulation. Without the small error 
in compressor pressure ratio, the response to the step (figs. 11(a) and (b)) is well be- 
haved and approaches the steady-state value calculated by the simulation (steady-state 
part). Also, the transients show no change from 0 to 0.1 second. 

When the small error in the compressor-pressure ratio is introduced the startup 
transient becomes much larger than the input. Figure 11(a) shows that at 0. 1 second 
there is a loss of 5 rpm from the operating point. The error built up is large when the 
step is applied. 

However, the speed never recovers to the predicted steady-state value. If this 
transient is used to compare with the linear model, the comparison will be very poor. 
The linear model is derived from the steady-state part of the simulation. Thus, a pos- 
itive step in nozzle area would match the curve without the startup transient. In effect, 
the comparison would be what is shown in figure 11(a). 

Figure 11(b) shows the same transients but for fan speed. At 0. 1 second there is 
little difference in the two transients (the error input was into the compressor). But at 
0 . 1 second the curves diverge since there is a large error buildup in the compressor- 
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turbine match. The fan speed must change to reflect this error, and thus when there is 
a starting transient, the fan speed does not recover to the predicted steady-state value. 
Again, when the error is eliminated, the speed behaves well. 

The results of this analysis show that great care must be used to make sure that 
startup transients have been eliminated when validating linear models that require small 
inputs . 
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APPENDIX G 


FOURTH-ORDER MODELS 

Presented in this appendix are the ABR, BBR, CBR, and DBR for the reduced- 
order models. The models are 

X = ABR X + BBR U (Gl) 

Y = CBR X + DBR U (G2) 

Also presented are the eigenvalues for the ABR matrices 

|abR-\T| =0 (G3) 

The six operating points are the same as listed in appendix D and are shown again 
in table VII for reference. 
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Operating Point 1 (sea-level-static idle) 


MATRIX 3P = 



1 

2 

3 

4 

1 

-2.149 

0. 6785 

85.74 

-662. 7 

2 

0. 7537 D- 01 

-5. <07 

12 3.8 

-145.6 

3 

0 . 86 82 D- 01 

u.924 

-14Q.1 

25.70 

4 

0. 13^9 

-0. 1935 

5.714 

-243.6 


MATE 

1 

IX 33R = 

1 

1.275 

2 

15.3 1 

3 

- 14. 06 

4 

6.183 

5 

0.5555D 

05 

2 

1.038 

-6 .352 

-1. 571 

-34.67 

0. 5512D 

05 

3 

1.099 

2.528 

2. 329 

39.43 

4298. 


4 

0 . 4224D-01 

-241.1 

1. 240 

-1.610 

0. 4163D 

05 


MATRIX CBR = 

1 

2 

3 

4 

1 

0.2275 

-0. 3059 

9. 286 

372.8 

2 

0. 1749D-01 

0. 2800D-02 

-0. 1506D-01 

-1.143 

3 

0.31 65D-02 

-0.4691 

5. 903 

1.989 

4 

-0. 24 8 4 D- 04 

0. 3969D-04 

- C . 4 8 2 3D- 03 

- 0. 3662D-C 1 

5 

0. 1650D-04 

0. 5936D-03 

-C.1473D-01 

0. 36 2 3D-0 2 

6 

0. 50 9 7 D- 04 

0. 6505D-05 

C.4075D-04 

-0.271 0D- 0 1 

7 

0.211 ID- 04 

-0. 2884D-C5 

0.3373D-04 

-0.5 188D-0 2 


MATRIX DBR = 

1 

2 

3 

4 

5 

1 

0. 5337D-01 

-0.4695 

2. 0 18 

-2.628 

0. 17 67D 05 

2 

0 . 16 79 D- 0 3 

C. 45 1 ID- 02 

0.1664 

0. 2495 D- 01 

148.0 

3 

0 .6032 

-0.7252 

-0.2431 

-4.857 

1841. 

4 

0. 7166D-05 

0. 6971D-C3 

-0. 1359D- 02 

Q.7972D-03 

4.720 

5 

-0. 1484D-04 

-C. 2731D-02 

0.3182D-03 

0 .7678D-02 

-3.855 

6 

0. 4533D-05 

0. 354 2D- 0 3 

0 .576 CD- 04 

0.7493D-04 

3.817 

7 

0. 672CD-06 

C. 3846D-04 

0 .2092D-03 

-0 . 207 ID-04 

0.7410 


Eigenvalues 


PEAL PART 

-0.87832969 

-2.6616826 

-142.80557 

-245.00836 


IMAGINARY PART 

0.00000000 
0 .00000000 
0.00000000 
0.00000000 
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Operating Point 2 


MATRIX ABR = 

1 2 

1 -3.684 -0.6320 

2 1.200 -6.137 

3 -1.659 5.819 

4 0.2584 -0.2028 


3 4 

103.0 -789.0 

79.95 -522.8 

-154.8 729.0 

5.581 -103.7 


MATRIX B3R = 

1 2 

1 0.7260 -257.7 

2 0.6320 -458.4 

3 0.8555 647.3 

4 0. 1891 D- 01 -1229. 


3 

-107.9 
13.40 
-16. 04 
10.16 


4 

-1.972 

-73.93 

59.29 

-1.938 


5 

-0. 1500D 05 
-0.1327D 05 
0.7941D 05 
-399.6 


MATRIX CBR = 


1 

2 

3 

4 

5 

6 
7 


1 

1.109 

0. 1382D-01 
0. 1093 
0.731 5D-04 
-0. 3066D-04 
0 .4003D-04 
0. 1332D-04 


2 

-0. 8876 
0. 3152D-05 
-0. 1181 
0.516 9D-05 
0. 1 294D-03 
-0. 3016D-04 
C.4075D-05 


3 

24.65 

-0.7751 D-07 
1. 136 

0 . 5422D-06 
-0. 2594D-02 
0.11 66 D-03 
-C.9470D-05 


4 

70.15 

-0. 9650D-02 
-24.47 
-0. 1508D-0 1 
0. 1449D-01 
-0. 1430D-01 
-0.3745D-02 


MATRIX DBR = 


1 

2 

3 

4 

5 

6 
7 


1 

0. 1647 
0 . 6344D-04 
0. 1329 
0.3239D-05 
-0. 8887 D- 08 
0.31 1 1D-07 
0 . 427 1 D-07 


2 

-232. 1 
0.1782 
-23.97 
-C. 3606D-0 2 
0. 1247D-01 
-0. 1163D-01 
-0. 2992D-02 


3 

41 . 22 
0.7403 

1.113 

-0.31 10D-02 
-0.291 5D- 03 
0.1737D-02 
0.7025D-03 


4 

-8.291 
0. 2073D-03 
-1.199 
0 . 6736D-04 
0.2437D-02 
0. 3235D-03 
0.5 1 75D-C4 


5 

-3459. 

0. 7989D-02 
-2464. 

0. 1239D-01 
-0. 1710D-01 
0. 5144D-02 
0.2557D-02 


Eigenvalues 


REAL' PART IMAGINARY PART 


-3. 1516863 
-5.5467604 
-60. 662508 
-198.9 1636 


0.00000000 

0.00000000 

0.00000000 

0.00000000 
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Operating Point 3 


MATRIX ABR = 
1 

2 

3 

4 

1 -7.879 

0.9396 

98. 10 

-592.0 

2 -3.031 

-8.405 

67.88 

-97.23 

3 5.137 

7.877 

-152.2 

159.8 

4 0.2454 

-0. 3006 

5.846 

-97.76 


MATRIX BBR = 
1 

2 

3 

4 

5 

1 0.6766 

-325.5 

-263.6 

26.37 

-6251. 

2 0.6269 

105.6 

- 139.0 

-67.69 

-5763. 

3 0.8802 

-98.04 

225.4 

109. 1 

0.3298D 05 

4 0. 4199 D- 01 

-1675. 

13. 53 

-5.122 

-349.4 


MATRIX CBR = 
1 

2 

3 

4 

1 

-0. 1689 

-1.501 

28.78 

154.6 

2 

0. 3305D-01 

0. 6804D-03 

-0.17 09 D-0 2 

-0. 3215D-01 

3 

-0. 7572 D- 01 

-0. 8768D-0 1 

0.8535 

-3.239 

4 

0. 3228D-04 

0. 4222D-04 

-0.1078D-03 

-0. 1977D-02 

5 

0. 7168 D- 04 

0. 8988D-04 

-0.1622D-02 

0. 2062 D- 02 

6 

0. 1824D-04 

0.1 069D-04 

-0 . 26 67D-04 

-0 .3672D-02 

7 

0 . 92 30 D- 06 

0. 9995D-05 

-0.2572D-04 

-0.91 07D-03 


MATRIX DBR * 
1 

2 

3 

4 

5 

1 

0.1879 

-183.5 

-8. 129 

-25.23 

-1559. 

2 

0 . 2079D- 03 

0.9867 

1.794 

0. 1238 D-0 1 

-2.719 

3 

0. 85 09 D-0 1 

5.206 

-4. 855 

-2.705 

-576.3 

4 

0.4685D-06 

0.456 ID-02 

-0.3128D-02 

0.7268D-03 

-0. 4405D-02 

5 

0. 15 1 9 D- 06 

-0. 8408D-03 

0.2865D-02 

0.2832D-02 

-0.7376D-02 

6 

0 . 2653D-06 

0.5091D-02 

-0. 4673 D- 03 

0. 1529D-03 

0. 4778D-02 

7 

-0. 3444D-08 

0. 9200D-03 

0.4092D-03 

0. 1885D-03 

-0.9336D-03 


Eigenvalues 


REAL PART IMAGINARY PART 


-5. 9363691 
-5.9363691 
-81.667400 
-172.74548 


1.8557230 

-1.8557230 

0.00000000 

0.00000000 
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Operating Point 4 


MATRIX ABB = 

1 

2 

3 

4 


1 -0.9575 

-0. 1766 

106.5 

- 877.5 


2 0.3233 

-1 .662 

80. 38 

-503.0 


3 -0.4455 

1.536 

- 150.4 

680.9 


4 0.7116D-01 

- C. 54S3D-0 1 

5. 52b 

-97. 28 


MATRIX BBR = 

1 

2 

3 

4 

5 

1 0.7S61 

-48.54 

-28 .97 

-0.6660 

-9199. 

2 0.6307 

1 05. C 

6.0 19 

- 19.68 

-7538. 

3 0.3F35 

-135.1 

-7. 237 

15.11 

0.4364D 05 

4 0.3E87D-01 

-296.0 

2. 321 

-0.5710 

-491 .7 


MATRIX CBF. = 
1 

2 

3 

4 

1 

0 . 2564 

-C. 2845 

28. 84 

- 1.747 

2 

0 . 4 1 73D- 02 

0. 1 307D-05 

- C .7787 D-0 6 

-0. 1404D-C 1 

3 

0.1136 

-C. 1270 

4. 087 

-88.94 

4 

0 . 155CD-03 

0. 6621 D- 05 

0.5959D-C5 

-0.6216D-01 

C 

-0.374 ID-04 

0. 146 ID-03 

-0. 1065D-01 

0.591 3D- 0 1 

b 

0.4192D-04 

-G. 3C76D-C4 

C.4211D-C3 

-0.5070D-0 1 

7 

0. 1366D-04 

0.4574D-05 

-0.3548D-04 

-0. 1314D-C1 


MATRIX DBR = 
1 

2 

3 

4 

5 

1 

0 . 1577 

74.45 

7.321 

-2.664 

-2472. 


0 . 9665D- 04 

0. 1070 

C . 173 c 

0. 17 64 D-0 3 

-3. 120 

3 

0.5122 

19.07 

1. 758 

-1.358 

-5532. 

L 

0. 1305D-04 

C .2140D-C 1 

-C.4331D-02 

0 . 706 9 D-0 4 

-0.7896 

5 

0. 1503D-05 

-C. 1 133D-01 

-0.61 18D-03 

0. 1966D-02 

-0. 1 4 87 D -0 1 

6 

-0 . 1225D-06 

C. 100 ID-01 

C . 1576D-02 

-0. 3559D-C 3 

-0 . 5224 D- 0 1 

7 

0. 21 30D-06 

0. 2575D-C2 

0.5 367D-C3 

0.3574D-C4 

-0. 1348D-01 


Eigenvalues 


RLAL PART 

- f . 6 6 6 t 00 3 1 
-1.5322735 
-56 . 951700 
-190. 5343a 


IMAGINARY PARI 

0.00000000 
0.00000000 
0 . 0C0O0"OO 
0. 00000000 
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Operating Point 5 


MATRIX ABR = 

1 

2 

3 

4 

1 -1.961 

-0.3216 

104. 1 

-839.4 

2 0.6604 

-3.232 

79.31 

-508.2 

3 -0.8983 

3.016 

-151 . 3 

696.9 

4 0.1382 

-0. 1064 

5.552 

-98.96 


MATRIX BBR = 
1 

2 

3 

4 

5 

1 0.7680 

-147.4 

-54.62 

0.3312 

-9176. 

2 0.5902 

C.7827 

10.89 

-36.39 

-9957. 

3 0.9134 

-32.71 

-13.01 

27.99 

0.5235D 05 

4 0 . 3266D-01 

-607.4 

4.467 

-0.8482 

-652.6 


M A' 

TRIX CBR = 

1 

2 

3 

4 

1 

0 . 5007 

-0.5556 

28. 88 

-6.866 

2 

0. 7768D-02 

0. 3799D-05 

-C. 1414D-04 

-0. 1636 D- 01 

3 

0.1132 

-0. 1197 

2.067 

-44.95 

n 

0. 1900D-03 

0. 6005D-05 

0.8149D-06 

-0.3027D-01 

u 

-0.3545D-C4 

0. 1350D-03 

-0. 5067 D- 02 

0. 2853D-01 

(' 

0 . 4068D-C4 

-0. 3010D-04 

0.2080D-03 

- 0 . 2589D-01 

7 

0. 131 2D- 04 

0.4535D-05 

-0. 1877D-0L 

-0.6749D-02 


MATRIX DdR = 
1 

2 

3 

4 

5 

1 

0. 18b3 

139.2 

13. 86 

-4.317 

-3131. 

2 

0. 89 1SD-04 

0. 1 345 

0.329 1 

0.274 ID-04 

-0.1415 

3 

0.25 12 

0.9063 

1.563 

-1.113 

-3236. 

4 

-0. 8940D-05 

-0. 19 8 3D- 0 1 

- 0 . 5 1 0 ID- 0 3 

0.9520D-C4 

-0. 1344D-01 

c. 

0. 2426D-05 

-0. 5765D-03 

-0.5172D-03 

0. 2271D-02 

0. 5789D-01 

6 

-0. 8103D-06 

0. 3437D-02 

0. 1529D-02 

- 0 . 2 852D- 03 

-0 . 5343D-0 1 

7 

0 . 1269D-07 

0. 1 166D-02 

C.5840D-03 

0.540 ID-04 

-0. 1107D-01 


Eigenvalues 


REAL PART IMAGINARY PART 


1.667 97 6 8 
.3 .0882016 
57. 6^3485 
19 3. 06924 


0. oooocooo 
0.00000000 
0.00000000 
o .oooocooo 



Operating Point 6 


MATRIX ABR = 

1 

2 

3 

4 

1 -1.159 

0. 2999 

91.41 

-553.4 

2 -0.4326D-01 

-2.154 

63. 36 

-138.0 

3 0.2093 

1.827 

-107.2 

114.5 

4 0.74 1 ID- 01 

-C.6340D-01 

4. 156 

-77. 34 


MATRIX BBE = 

1 

2 

3 

4 

5 

1 1.176 

-15.01 

-10. 02 

4. 136 

-8622. 

2 0.9912 

-3.966 

-4.557 

-16.10 

-4955. 

3 1.031 

14.07 

6. 113 

20.40 

0.3420D 05 

<■ 0.4048D-01 

-189.8 

0.7533 

-0.8214 

-344.0 


MA 

TRIX CBR = 





1 

2 

3 

4 

1 

-0. 2873D-02 

-0. 2996 

16.83 

238.6 

2 

0 . 1089D-01 

0. 1620D-02 

-0 . 3626 D-02 

-0.7139 

3 

-0. 1834 D-0 1 

-0. 1824 

4.257 

-16.35 

U 

0. 1275D-C3 

0. 1 1C1D-03 

-C.243CD-Q3 

-0 . 48 26 D-C 1 

0 

0. 4972D-04 

0. 2068D-03 

-0.1462D-01 

0.2112D-01 

6 

C . 901 2D- 04 

0. 1431D-C4 

-0.4249D-C4 

-0.697 CD- 01 

7 

0. 3266D-04 

0. 4030D-C5 

-0.97 17D-05 

- 0 . 15 83D- 01 


MATRIX DBR = 

1 

2 

3 

4 

5 

1 

0. 15 94 

-47.58 

-0.7485 

-3.780 

-1127. 

2 

0. 1975D-03 

-0.6159D-C1 

C .1365 

0.2133D-01 

1 . 668 

3 

0.7961 

-2.578 

-1.016 

-3.319 

-4555. 

4 

0. 3630D-05 

-0.7 057D-02 

-0.4148 D-02 

0. 1479D-02 

0.4707 

5 

0. 4405 D- 05 

0. 2 8 86 D- 02 

0.9740D-03 

0.41 98D-02 

-0.2233 

6 

0. 21 18D-06 

-0.6671 D-02 

0 . 67 20 D-0 4 

0. 1953D-03 

-0. 9091D-01 

j 

0. 2037D-06 

-0. 1510D-02 

0.308QD- 03 

0.5381D-04 

-0.9650D-02 


Eigenvalues 


REAL PART IMAGINARY PART 

-1.0627800 0.00000000 
-1.5121531 0.00000000 
-65.526587 0.00000000 
-119.70158 0.00000000 
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APPENDIX H 


COMPARISON OF FULL-STATE AND REDUCED-ORDER (4th ORDER) 

MODEL TRANSIENT RESPONSES 

Transients run with the full-state and reduced-order linear models are compared 
herein. The responses of reconstructed states and outputs are shown in figures 12 
and 13 for a +3-percent change in fuel flow and nozzle area. The operating point is sea- 
level-static idle. (Plots shown in the main text are omitted.) The responses shown are 

for Tipg 5 JJ* ^T2 5C* ^TS’ I ^'t4* ^TS’ ^TGC’ ^TTM’ ^VFAN, SMHC, 

DP g 5 Qp, and DPg 55c • I n the comparisons for the reconstructed states, there are 
immediate responses at zero time . This is due to the loss of frequency information 
when generating a reduced-order model. The discarded states are approximated by 
X = 0, and then reconstructed from linear relations with the retained states an inputs. 
The loss of high-frequency information is observed by the faster rise time for the re- 
constructed states. In general, after 0. 1 second the reconstructed states from the 
reduced-order model agree well with the full-state model . Immediate responses for the 
outputs (figs. 12 and 13(h) to (m)) are inherent in the models (ref. 8). The output var- 
iables can change instantaneously with time through the D matrix. In the comparisons 
of output responses, the offset remained the same for the reduced-order model. The 
difference in rise time between the output responses can be attributed to loss of high- 
frequency information. Again, after 0.1 second, the fourth-order responses approxi- 
mate the full-state responses quite well. 

The same procedure was used at the other five operating points, and the agreement 
between transient responses was as good or better than that shown for the sea-level- 
static-idle operating point. From this analysis the fourth-order models can be used for 
control studies since they adequately represent the engine over its whole operating 
range . 
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TABLE I. - TYPICAL ENGINE FLIGHT CONDITIONS 


Ope rating 
point 

Flight condition 

Significance 

1 

Sea-level takeoff 
(SLTO) 

Maximum nonaugmented (intermediate) 
thrust and maximum augmented thrust 
rating points (T T4 , stability, thrust 
primary importance). 

2 

Sea-level ram 

Burner case pressure limit as well as 
extremely high fan and compressor 
blade and vane loading along this 
boundary. High structural and 
vibratory stress. High fan and com- 
pressor speeds may be reached here. 

3 

Supersonic dash 

High steady-state inlet distoration and 
subsequent reduced stability margin. 
Possible inlet disturbance due to 
popping or swallowing the inlet shock 
waves during maneuvers. 

4 

i ~ 

Extreme envelope 

definition 

Control must satisfy performance, sta- 
bility, and safety limits. Also distor- 
tion and augmentor light problems. In 
addition, may reach or exceed fan and 
compressor speed limits. 

3 

Low-pressure 
limited- flight 
condition 

Ignition limits for both main burner 
and augmentor occur along this bound- 
ary. Flame out and igntion instability 
possible. 

G 

Subsonic-cruise 
and combat 
maneuvering 
condition 

Thrust specific fuel consumption 
(TSFC - thrust per unit mass of fuel) 
and therefore fuel flow and F/A ratio 
important. Transient engine response 
characteristics during combat mis- 
sion. Also severe transient inlet dis- 
tortion due to aircraft attitude changes. 



TABLE n. - EIGENVALUES FOR TWO DIFFERENT PERTURBATION SIZES AT THE SEA-LEVEL-STATIC-IDLE OPERATING POINT 


(a) Positive perturbations only (b) Positive and negative perturbations 


Change in state, AX, percent 


Change in state, AX, percent 

+0 

.5 

+1 

0 



+0 

.5 


+1.0 

Eigenvalues of A matrix 


Eigenvalues of A matrix 

Real 

Imaginary 

Real 

Imaginary 


Real 

Imaginary 

Real 

Imaginary 

-860.24312 

0 

-849.45871 


0 


-852.11287 


3 

-869.70083 


0 

-267.76119 

0 

-264.38152 


0 


-268.16126 



-422.82995 


0 

-106.45859 

0 

-118.70167 




-135.01001 



-133.54989 


3 

-52.072414 

24.775151 

-50.905953 

15.076307 


-56.317567 



-50.563185 


7.5119157 

-52.072414 

-24.775151 

-50.905953 

-15.076307 


-49.669363 



-50.563185 

-7.5119157 

-51.736044 

0 

-51.779746 


3 


-48.432318 



-52.253414 


3 

-43.674746 

0 

-43.546884 




-40.376271 



-42.606607 



-28.475954 

0 

-26.871969 




-14.119784 



-23.852079 



-4.8312091 

2.4935793 

-8.1525361 




-2.8587858 



-12.163400 



-4.8312091 

-2.4935793 

-3.9743200 

1 



-21.319590 


.83925693 

-19.973461 


.52786244 

-20.034792 

.38079509 

-20.057634 


.23031854 


-21.319590 

-.83925693 

-19.973461 


-.52786244 

-20.034792 

-.38079500 

-20.057634 


-.23031854 


-19.747053 

C 


-.74937136 


.24305228 

-2.0269044 

0 

-2.0128254 

0 


-1.9813104 

c 


-.74937136 


-.24305228 

-.68641882 

.22144251 

-.69555555 


.20482900 


-.74410687 


.20040877 

-2.4917681 

0 

-.68641882 

-.22144251 

-.69555555 


-.20482900 


-.74410687 

- 

.20040877 

-1.9956926 

0 

-21.033073 

0 

-20.914347 

( 



-20.347614 

0 


-21.620234 

C 



TABLE HI. - EIGENVALUES FOR 
SEA- LEVEL-STATIC-IDLE 


OPERATING POINT 


Vector 

Eigenvalues of 

A matrix 


Real 

In' aginary 

1 

-852.11287 


0 

2 

-268.16126 



3 

-135.01001 



4 

-56.317567 



5 

-49,669363 



6 

-48.432318 



7 

-40.376271 



8 

-14.119704 



9 

-2 . 8587858 



10 

-21.319598 


.83925693 

11 

-21,319590 


-.83925693 

12 

-19,747853 

0 

13 

-1.9813104 

0 

14 

-.74410687 


. 20040877 

15 

-.74410687 


-.20040877 

16 

-20.347614 

c 
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TABLE IV. - MATRIX OF MODE SHAPES, IN MAGNITUDE AND ANGLE PEG) FORM FOR SEA-LEVEL-STATIC-IDLE 

OPERATING POINT 


State 









or 1 

VECTOR 1 

2 

3 

4 

5 

6 

7 

8 

1 

- 0 . 50 91 D- 02 

0.111 2D-0 1 

-0.1061 D-Q 1 

-0. 1119D-01 

-C.3240D-02 

-0.1674D-01 

0 • 8529D-02 

0.3946D-01 

2 

0.6996D-03 

0.4697D-02 

-0.5217D-02 

-0. 3222D-02 

0. 6876D-04 

0 . 1756 D-02 

-0. 1572D-01 

-0. 9271D-02 

3 

0.4265D-02 

-0.6182 

1.000 

-0.4488D-01 

-0.2938D-03 

-0.4445D-02 

-0. 8495 D-02 

-0.1791 

4 

1.000 

0.4298 

0.6677 

-0. 4734D-01 

0.5837D-03 

0.5958D-02 

-0.4413D-01 

-0 . 285 1D-0 1 

5 

-0.9348D-01 

1.000 

0.2819 

0.27 4 ID- 02 

0.1060D-02 

-0.160 lD-04 

-0. 1Q31D-03 

0.683CD-0 1 

6 

0 . 4284D- 03 

-0.1C59D-01 

-0.8042D-02 

0.57 23D-03 

-C.7569D-05 

0. 1542D-02 

-0.3850D-02 

0.1023 

7 

0.5358D-03 

-0. 137 ID- 0 1 

-0.1104D-01 

0. 1607D-03 

0.11 03D-03 

0. 1434D-02 

-0.5 Q78D-02 

0.1360 

8 

- 0 . 3 148D- 03 

0. 1986D-01 

-0.3447D-01 

0.5197D-02 

- 0 . 3029D- 04 

-0. 1337D-02 

0. 7 846 D-02 

0.1324 

9 

0 . 955 1 D- 04 

-0.15b7D-01 

0 . 1 978 D-0 1 

0. 2910 

-0 . 46 84 D- 02 

-0.1246 

1.000 

0.3339 

10 

0. 1329D-04 

-0. 2056D-0 2 

0.2554D-02 

0 . 3579D- 0 1 

-0.5702D-03 

-0.1512D-01 

0.1188 

0. 302 1D-0 1 

11 

-0. 17 91 D- 02 

0.7303D-01 

-0.3301 D-0 1 

-0. 4046 D- Cl 

0.2507D-04 

-0.4099D-02 

0.2084 

0.2486 

12 

*0. 1639D-01 

-0. 9005D-0 1 

0.9364D-01 

1.000 

0.2054 

1.000 

-0. 1710 

0.3722 

13 

- 0. 6535D-02 

-0.3562D-C1 

C. 3660 D- 01 

0. 3781 

C.77C4D-01 

0.3743 

-0.6317D-01 

0.1121 

14 

-0.1167 D-02 

-0. 4395 D-G 1 

0.4630D-01 

0.1088 

0. 8835D- 02 

0.2722D-01 

0. 1 306D-0 1 

0.9822 

1 5 

0 • 9337 D- 06 

0. 5654D-03 

C.2050D-G2 

- 0 . 6855D-04 

-0.8388D-04 

-0 . 1040 D-02 

0.4891 D-02 

0.4569 

16 

-0.1 959D- 01 

-0.251 1D-0 1 

-0.1308 

-0.5305 

1.000 

0.6876 

-0.7342D-01 

1.000 


9 

10 

11 

12 

13 

14 

15 

16 

1 

1.000 

0. 1108D-01 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

2 

-0.2465D-01 

0.4686 D-02 

G.110CD 08 

0. 1200 D C8 

0.13C0D C8 

0. 1400D 08 

0. 1500D 08 

0.U0CD 08 

3 

-0. 1354 

0. 1584 

-121.3 

-0. 29 11 D- 03 

0 • 8936D-01 

1.000 

0.0000 

0.2817D-02 

4 

0.5330D-01 

0.9993D-01 

-90.69 

-0. 1233D-03 

-0. 3463 D-02 

0.2065 

7.058 

0.1992D-02 

5 

0. 1493 

C. 5802D-0 1 

29.69 

0. 1331D-G1 

-0. 1498D-01 

0.9429 

6.993 

-0.6495D-02 

6 

0.2047 

0.4624D-01 

21.40 

0.2282D-C1 

0 • 3274D- 02 

0.4792 

4.804 

0.1 178D-0 1 

7 

0.2076 

C. 468 1 D-C 1 

1.219 

0.3413D-01 

0. 1278D-01 

0.1468 

-0.1646 

0.3 1 02D-0 1 

8 

0.6082D-01 

0. 1276 

-149.6 

-0.31 97D-02 

0. 1705D-01 

0.1331 

-2.988 

0.1911 D-02 

9 

0.1172 

0.1364 

-148.5 

-0. 1164D-01 

0. 1767D-01 

0.1073 

-5.304 

0. 1674D-0 1 

10 

-0. 1521D-01 

0. 1438D-0 1 

58.19 

0.1637D-01 

0.3054D-02 

0.2734 

4.476 

-0 . 2850D-0 1 

11 

0.1146 

0.7760D-01 

101.8 

0. 1632D-0 1 

0. 9846 D- 02 

0.2582 

-171.5 

-0.4952D-01 

12 

0.1764 

0.1071 

101.2 

0 . 1 686D-02 

-0.3006D-02 

0.6820 

80.10 

-0.5149D-02 

13 

-0.1750 

0. 3622D- 0 1 

109.1 

0 • 943 1 D- 02 

0 . 9992D-02 

0.2990 

-171.6 

-0. 2933D-0 1 

14 

0.2362 

0.7247 

112.6 

0. 1905D-C1 

0. 1547D-01 

0.3887 

-168.3 

-0.3569B-0 1 

15 

0.2336 

0.5801 

112.1 

0.6321D-02 

1 .oco 

0.5185 

-174.6 

-0.11 95D-0 1 

16 

-0.1137 

1.000 

-22.76 

1.000 

0. 1178D-01 

0.4734 

-169.2 

1.000 


TABLE V. - MODAL REDUCTION AT SEA-LEVEL-STATIC-IDLE 


OPERATING POINT 


(a) States and corresponding 
eigenvalues 


States 

Eigenvalues, ra d/sec 

1 

-2.86, -0 . 744±0 . 20 

2 

-0 . 744 =t0 . 20 

3 

-135.01, -0 . 744±0 . 20 

5 

-268.16 


(b) Possible models 
(states 1, 2, 3, and 5) 


Model 


Eigenvalues, rad/sec 


1 -2.86, -0 . 744±0 . 20, -268.16 


2 -0 . 744±0 . 20, -135.01, -268.16 


TABLE VI. - OPE RATING- POINT DEFINITION 


FOR 16TH-ORDER LINEAR MODELS 


Operating 

point 

Power 

lever 

angle, 

deg 

Altitude, 

km 

Mach 

number 

Data 

(pages) 

1 

20 

0 

0 

33 - 35 

2 

83 

0 

0 

36 - 38 

3 



0 

1.2 

39 - 41 

4 



13.7 

. 9 

42 - 44 

5 



9.14 

.9 

45 - 47 

6 

20 

9.14 

,9 

48 - 50 
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TABLE VII. - OPERATING- POINT DEFINTION 


FOR FOURTH-ORDER LINEAR MODELS 


Operating 

Power 

Altitude, 

Mach 

Data 

point 

lever 

angle, 

deg 

km 

number 

(pages) 

1 

20 

0 

0 

55 

2 

83 

0 

0 

56 

3 



0 

1.2 

57 

4 



13.7 

. 9 

18 

5 

1 


9.14 

. 9 

59 

6 

20 

9.14 

. 9 

60 
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pressure 
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-Envelope 

boundary 


-Blade loading and 
burner case limits 



Figure 3. - Linear engine model. 


Mach number 

Figure 2 . - Typical engine operating envelope with constraints. 




Perturbation, 

percent 

O +0.5 

□ +1.0 

▲ Eigenvalues coincide 


-4 -2 -1 6 -20 -10 -6 -4 -2 -1 

Real part of eigenvalue, a, rad/sec 

(a) Positive perturbation. (b) Positive and negative perturbations averaged. 

Figure 4. - Eigenvalue shift (upper left-half plane). 
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Change in compressor discharge pressure, P-p, kPa Change in fan speed, Nj_, rpm 




Full-state model 

Reduced-order model 



(e) Augmentor pressure. 



L. .1 1 1 1 J 

.75 LOO 1.25 1.50 1.75 2.00 

Time, sec 


(f) Engine net thrust. 
Figure 7. - Concluded. 



Change in compressor-discharge pressure, kPa Change in fan speed, N|_, rpm 



(c) Compressor-discharge pressure. (d) Interturbine pressure. 

Figure 8. - Comparison of full-state and reduced-order linear-model responses to a +3 percent step in nozzle area; sea-level-static-idle 
operation. 
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! 


17.5i — 3.51- 



Time, sec 

(c) Augmentor pressure. (d) Fan inside-diameter discharge pressure. 

Figure 9. - Comparison of nonlinear- and linear-model responses to a +3 percent step in fuel flow; sea -level -static -idle operating point. 



i 



50 . 75 1.00 1.25 1.50 1.75 2.00 0 . 25. .50 . 75 1.00 1.25 1.50 1.75 


Time, sec 

(g) Fan-turbine-exit temperature. (h) Duct-exit temperature. 


Figure 9. - Continued. 



Change in burner-exit temperature, T-j- 4 , K Change in exhaust-nozzle-inlet temperature, Tjy^, 



(k) Burner-exit temperature. ( 1 ) Fan stall margin. 

Figure 9. - Continued. 


Nonlinear model 

Linear model 



J L I I I I I 

(m) Compressor stall margin. 



(n) Empirical fan-exit AP/P parameter. 


Figure 9. - Concluded. 


(o) Theoretical fan-exit AP/P parameter. 



Change in augmentor pressure, F^, kPa Change in compressor-discharge pressure, P^, kPa 



Figure 10. - Comparison of nonlinear- and linear-model responses to a +3 percent step in nozzle area; sea-level-static-idle operating 
point. 


Change in fan-turbine-exit temperature, T T5 , K Change in duct temperature, 




1.0 — 
.5 — 



•2 5 1 1 I I 1 1 I 

' 0 .25 .50 .75 LOO 1.25 1.50 



(g) Fan-turbine-exit temperature. 


(h) Duct-exit temperature. 


Figure 10. - Continued. 
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Change in burner-exit temperature, T T4 , K Change in exhaust-nozzle-inlet temperature, Tjy M , 
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Figure 10. - Continued. 
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(l) Fan stall margin. 
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Change in compressor-discharge temperature, T T3 , K Change in fan inside-diameter discharge temperature, T T2 5H> 




.9 


— - Full-state model 
— Reduced-order model 


.8 



.50 .75 LOO 1.25 1.50 1.75 2.00 0 .25 .50 .75 1.00 1.25 1.50 

Time, sec 


(g) Exhaust-nozzle-inlet temperature. <h) Total engine airflow. 

Figure 12. - Continued. 
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